Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,382,332 A→T 10.9% intergenic (+12/+15) dinB → / ← PP_1204 DNA polymerase IV/lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,382,3320AT10.9% 184.7 / 13.6 64intergenic (+12/+15)dinB/PP_1204DNA polymerase IV/lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (3/4);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

AGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTC  >  NC_002947/1382192‑1382475
                                                                                                                                            |                                                                                                                                               
agTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGcccc                                                                                                                                        <  1:998046/150‑1 (MQ=255)
  ttATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGCCCCtt                                                                                                                                      <  2:767770/150‑1 (MQ=255)
  ttATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGCCCCt                                                                                                                                       >  1:638318/1‑149 (MQ=255)
  ttATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGCCCCt                                                                                                                                       <  2:638318/149‑1 (MQ=255)
     tCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCGCTTCTa                                                                                                                                   >  1:909089/1‑150 (MQ=255)
     tCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTa                                                                                                                                   <  1:598122/150‑1 (MQ=255)
        gCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCTCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTa                                                                                                                                <  1:251960/150‑1 (MQ=255)
          aaTTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAAt                                                                                                                              >  1:1060792/1‑150 (MQ=255)
          aaTTGCAGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAAt                                                                                                                              >  1:443865/1‑150 (MQ=255)
                          gTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGCCCCtttca                                                                                                                                   <  2:731461/129‑4 (MQ=255)
                          gTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCATTGCAGCCCCtttca                                                                                                                                   >  1:731461/1‑126 (MQ=255)
                            tCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGcca                                                                                                            <  2:838852/150‑1 (MQ=255)
                                ccGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGc                                                                                                                                           <  1:331931/115‑1 (MQ=255)
                                ccGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGc                                                                                                                                           >  2:331931/1‑115 (MQ=255)
                                ccGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGGGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCaa                                                                                                        >  2:696252/1‑150 (MQ=255)
                                   tGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGATGTTTCCGCCCAAATGAAAGGGGCTGCAATGTAGCGCGTTCTATTAATGCACGCCTGGATCCGCCACCAACCg                                                                                                     >  1:569875/1‑150 (MQ=255)
                                       ggcAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCa                                                                                                 >  2:871795/1‑150 (MQ=255)
                                             ccGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAAt                                                                                                                              >  1:966565/1‑115 (MQ=255)
                                             ccGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAAt                                                                                                                              <  2:966565/115‑1 (MQ=255)
                                              cGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCt                                                                                          <  2:281008/150‑1 (MQ=255)
                                                gTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTg                                                                                        <  1:871795/150‑1 (MQ=255)
                                                      ttGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGc                                                                                  <  1:1085706/150‑1 (MQ=255)
                                                         cTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGAc                                                                               >  2:627080/1‑150 (MQ=255)
                                                                taggggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGa                                                                        <  1:152483/150‑1 (MQ=255)
                                                                    ggtGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTc                                                                    <  2:608167/150‑1 (MQ=255)
                                                                      tGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCct                                                                  >  2:614010/1‑150 (MQ=255)
                                                                       gAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCctg                                                                 <  2:633494/150‑1 (MQ=255)
                                                                               gTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCt                                                         >  2:1105177/1‑150 (MQ=255)
                                                                                 gACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGTATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTcc                                                       <  1:627080/150‑1 (MQ=255)
                                                                                 gACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTcc                                                       >  1:102325/1‑150 (MQ=255)
                                                                                    ttGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCa                                                                           >  1:897331/1‑127 (MQ=255)
                                                                                    ttGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGACCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCa                                                                           <  2:897331/127‑1 (MQ=255)
                                                                                           gcgcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCtgctgc                                                             <  2:998626/134‑1 (MQ=255)
                                                                                           gcgcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCtgctgc                                                             >  1:998626/1‑134 (MQ=255)
                                                                                             gcgcATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCgg                                           <  2:1060792/150‑1 (MQ=255)
                                                                                                 aTGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                   >  2:791538/1‑122 (MQ=255)
                                                                                                 aTGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                   <  1:791538/122‑1 (MQ=255)
                                                                                                 aTGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCa                                       >  2:425437/1‑150 (MQ=255)
                                                                                                  tGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAgg                                                                         <  1:190332/115‑1 (MQ=255)
                                                                                                  tGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAgg                                                                         >  2:190332/1‑115 (MQ=255)
                                                                                                    aaCAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGGCTGCAGGACCTCCTGCTGCAGCGCGGGCTCGTTGCGGGTCAGCt                                    <  2:443865/150‑1 (MQ=255)
                                                                                                     aCAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTc                                   >  1:568502/1‑150 (MQ=255)
                                                                                                         ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGa                                 >  2:445840/1‑148 (MQ=255)
                                                                                                         ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGa                                 <  1:445840/148‑1 (MQ=255)
                                                                                                         ttGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATc                               >  1:320103/1‑150 (MQ=255)
                                                                                                               cTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGctc                         >  1:579280/1‑150 (MQ=255)
                                                                                                                tGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGctct                        <  2:579280/150‑1 (MQ=255)
                                                                                                                 gTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGAt                                                   >  1:1121863/1‑122 (MQ=255)
                                                                                                                 gTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGAt                                                   <  2:1121863/122‑1 (MQ=255)
                                                                                                                 gTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTg                       <  1:425437/150‑1 (MQ=255)
                                                                                                                  tttCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGt                      >  1:790762/1‑150 (MQ=255)
                                                                                                                     ccGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTcc                   <  2:790762/150‑1 (MQ=255)
                                                                                                                        cccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                   >  2:1197658/1‑99 (MQ=255)
                                                                                                                        cccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTcc                                                                   <  1:1197658/99‑1 (MQ=255)
                                                                                                                         ccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCgg                                           >  2:711335/1‑122 (MQ=255)
                                                                                                                         ccAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCgg                                           <  1:711335/122‑1 (MQ=255)
                                                                                                                            aaTGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCAc            <  2:89384/150‑1 (MQ=255)
                                                                                                                             aTGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACt           <  2:102325/150‑1 (MQ=255)
                                                                                                                              tGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTg          <  1:1105177/150‑1 (MQ=255)
                                                                                                                              tGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCCCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTg          <  1:398880/150‑1 (MQ=255)
                                                                                                                                   ggggCTGCATTGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTc     >  1:466240/1‑150 (MQ=255)
                                                                                                                                   ggggCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTc     >  1:683189/1‑150 (MQ=255)
                                                                                                                                      gCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCctc  <  2:683189/150‑1 (MQ=255)
                                                                                                                                      gCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCctc  >  2:941281/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                               
AGTTATCGGCAATTGCTGGGGCAGGCGTTCGCCCGTGGCGGCAAGCCGGTGCGCTTGCTAGGGGTAGGGGTGAGGTTGCGTGACTTGCGGGGCGCGCATGAACAGTTGGAGCTGTTTCCGCCCAAATGAAAGGGGCTGCAATGCAGCCCCTTCTATTAATGCACGCCTGGATCCGCCACCAACCGGCCAGTGGCCTTGGTAAGCGACTGCAGGAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTC  >  NC_002947/1382192‑1382475

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