Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 596,677 C→T 5.9% Q360* (CAG→TAG)  cusC → copper/silver efflux system, outer membrane component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913596,6770CT 5.9% 215.3 / 2.6 68Q360* (CAG→TAG) cusCcopper/silver efflux system, outer membrane component
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base T (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

AGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGG  >  NC_000913/596533‑596820
                                                                                                                                                |                                                                                                                                               
aGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCcagcag                                                                                                                                               >  3:299259/1‑147 (MQ=255)
       tatGGACCGCAAGTAGCGATCTATCGTCATTATTTAACGCCAGAAGCGGGATGTGGAATTTTATACCCAAAATTGAGCTCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCgggg                                                                                                                                      <  3:35412/149‑3 (MQ=255)
       tatCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCggtg                                                                                                                                      >  7:60010/1‑149 (MQ=255)
         tCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCggtggt                                                                                                                                    <  3:352825/149‑1 (MQ=255)
                ccAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATATTATACCCAAAATTGAGATAACCAATTTTAATGACGGACGCAACCAGTCCAATCTGGATATCGCAGAAATTCGCGAGCAGCCTTCGGTGGTGAATTAt                                                                                                                             >  4:62228/1‑149 (MQ=255)
                 cAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATAACAAAAATTGATATACCCATTTTAAATGCCGTACGCACCCATGCCAATCTGGATATCGCCGAAATTAGAAAGCAGCCGTCGGTGTTGAATTAt                                                                                                                             >  3:67634/1‑148 (MQ=255)
                   gTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGaa                                                                                                                          <  8:60010/149‑1 (MQ=255)
                   gTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGaa                                                                                                                          <  3:236198/149‑1 (MQ=255)
                          tCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGaaa                                                                                                                    >  8:313695/1‑148 (MQ=255)
                              tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGaaa                                                                                                                    <  1:362149/144‑1 (MQ=255)
                              tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATc                                                                                                                 >  3:88824/1‑147 (MQ=255)
                              tCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCa                                                                                                               <  8:282139/149‑1 (MQ=255)
                                    ttattTAACGCCAGCAGCGGGATGTGGAATTTTATACACAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAAATCGACAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCAAGAACGc                                                                                                         >  3:3467/1‑149 (MQ=255)
                                       tttAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCtt                                                                                                      >  3:329136/1‑149 (MQ=255)
                                           aCGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTaaa                                                                                                  <  4:299259/149‑1 (MQ=255)
                                           aCGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTaaa                                                                                                  <  7:51575/149‑1 (MQ=255)
                                               cagcagCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTa                                                                                                                              >  5:366737/1‑117 (MQ=255)
                                               cagcagCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTa                                                                                                                              <  6:366737/117‑1 (MQ=255)
                                                agcagcGGGATGTAGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGt                                                                                             <  1:233520/149‑1 (MQ=255)
                                                 gcagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTg                                                                                            <  7:166815/149‑1 (MQ=255)
                                                  cagcGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTgg                                                                                           <  7:313695/149‑1 (MQ=255)
                                                       ggATGTGGAATTTTATACCCAAAATTGAGATCCACATTTTTAATGCCGGACGCAACCAGGCCAATCTGTATATCTCCTAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGAt                                                                                      <  3:62228/149‑1 (MQ=255)
                                                             ggAATTTTATACCCAAAATTGAGATCCCCATTTTTGATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGa                                                                                       <  1:159831/142‑1 (MQ=255)
                                                             ggAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGa                                                                                       >  2:159831/1‑142 (MQ=255)
                                                              gAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTgg                                                                               >  1:252920/1‑149 (MQ=255)
                                                                     ataCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTAc                                                                         <  6:23383/148‑1 (MQ=255)
                                                                      taCCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGt                                                                       >  3:192177/1‑149 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTaaa                                                                                                  >  1:133519/1‑119 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTaaa                                                                                                  <  2:133519/119‑1 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCa                                                                     >  8:292123/1‑148 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCa                                                                     >  3:267049/1‑148 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCa                                                                     >  6:184417/1‑148 (MQ=255)
                                                                         ccAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGACAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCa                                                                     >  2:173720/1‑148 (MQ=255)
                                                                          cAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGt                                                                                             <  1:36663/123‑1 (MQ=255)
                                                                          cAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGt                                                                                             >  2:36663/1‑123 (MQ=255)
                                                                           aaaaTTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAg                                                                  >  7:254799/1‑149 (MQ=255)
                                                                                    aTCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGa                                                         >  6:23059/1‑149 (MQ=255)
                                                                                     tCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                        >  3:4308/1‑149 (MQ=255)
                                                                                      ccccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                        >  3:307518/1‑148 (MQ=255)
                                                                                      ccccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                        >  3:312690/1‑148 (MQ=255)
                                                                                       cccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACgt                                                         <  7:270195/146‑2 (MQ=255)
                                                                                       cccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGa                                                         >  8:270195/1‑146 (MQ=255)
                                                                                       cccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATc                                                       >  7:155922/1‑148 (MQ=255)
                                                                                        ccATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCa                                                      >  4:294654/1‑148 (MQ=255)
                                                                                          aTTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAAt                                                   >  7:15128/1‑149 (MQ=255)
                                                                                              ttAATGCCGTACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGTAGCAGGCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCacaacgg                                                >  3:69391/1‑142 (MQ=255)
                                                                                              ttAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGTAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                        <  4:74191/140‑1 (MQ=255)
                                                                                              ttAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGTAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGAt                                                        >  3:74191/1‑140 (MQ=255)
                                                                                              ttAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGTAGCAGTCGGGGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGGCTGAACGAt                                                        <  4:69391/140‑1 (MQ=255)
                                                                                                        aCGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTATCGCCCagcagc                                     <  4:329136/149‑1 (MQ=255)
                                                                                                         cGCAACCAGTCAAATCTGGATATCGCAGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGTAAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCACAGCAGCg                                    <  4:3467/149‑1 (MQ=255)
                                                                                                         cGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCg                                    <  4:192177/149‑1 (MQ=255)
                                                                                                              ccAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATc                               <  4:88824/149‑1 (MQ=255)
                                                                                                                 ggCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTgg                                                                                           <  3:348027/86‑1 (MQ=255)
                                                                                                                 ggCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTgg                                                                                           >  4:348027/1‑86 (MQ=255)
                                                                                                                 ggCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTgg                            >  4:62330/1‑149 (MQ=255)
                                                                                                                    cAATCTGGATATCGCCGAAATTCGCCTGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAACGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATc                                                       <  1:235013/119‑1 (MQ=255)
                                                                                                                    cAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATc                                                       >  2:235013/1‑119 (MQ=255)
                                                                                                                    cAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCAt                         <  2:252920/149‑1 (MQ=255)
                                                                                                                       tCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGc                      >  8:168874/1‑149 (MQ=255)
                                                                                                                              atCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCACCCAGCAGCGTTATCTGGCATCGCt                     <  2:298114/143‑1 (MQ=255)
                                                                                                                              atCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCACCCAGCAGCGTTATCTGGCATCGCt                     >  1:298114/1‑143 (MQ=255)
                                                                                                                                 gCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTAc            >  4:331091/1‑149 (MQ=255)
                                                                                                                                  ccGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACt           >  2:39520/1‑149 (MQ=255)
                                                                                                                                   cGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACtt          <  5:23059/149‑1 (MQ=255)
                                                                                                                                    gAAATTCGCCATCATCAGTCGGTGGTCAACTATTAACAGAAACTCCAGAACGCCTTTAAAGAAGCTCCAGTTTCGCTGGAATTACTTAACAGACTTAACGAACAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACttt         <  4:67634/149‑1 (MQ=255)
                                                                                                                                     aaaTTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTg        <  3:294654/149‑1 (MQ=255)
                                                                                                                                       aTTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGc       >  2:286515/1‑148 (MQ=255)
                                                                                                                                           gCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACgg  <  8:15128/149‑1 (MQ=255)
                                                                                                                                             cagcagcagTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAgc                                               <  8:167242/102‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
AGCGGAATATCGACCGCCAGTAGCGATCTATCGTCATTATTTAACGCCAGCAGCGGGATGTGGAATTTTATACCCAAAATTGAGATCCCCATTTTTAATGCCGGACGCAACCAGGCCAATCTGGATATCGCCGAAATTCGCCAGCAGCAGTCGGTGGTGAATTATGAACAGAAAATCCAGAACGCCTTTAAAGAAGTGGCAGATGCGCTGGCATTACGTCAAAGCCTGAACGATCAAATTAGCGCCCAGCAGCGTTATCTGGCATCGCTGCAAATTACTTTGCAACGG  >  NC_000913/596533‑596820

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: