Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC29.7% 261.7 / 138.3 139intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (13/84);  new base C (2/39);  total (16/123)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 180 total aligned reads displayed.

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGC  >  NC_000913/2947317‑2947596
                                                                                                                                                      |                                                                                                                                  
 cAAATCGGTGATCGGGTTAAAAACAAGGGTGCATGAAAACGCGAGCGTAGTATAGTCCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa                                                                                                                                   <  5:1404312/151‑1 (MQ=18)
 cAAATCCGTGATGGTGTTAAAAAAAAGGCTGCATGAAAAATAGATAGTAGTATAGTGCGCATACACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaa                                                                                                                                   <  8:1373907/151‑1 (MQ=17)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                   >  2:220636/1‑150 (MQ=255)
 cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                   <  8:1412772/151‑1 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatt                                                                                                                                  >  6:1466776/1‑149 (MQ=255)
  aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGAAACCAATTatt                                                                                                                                  >  5:1223021/1‑149 (MQ=255)
   aaggtgTGATGGGGGTAAAATTAACGTTGCACGAAAAAGAGAGTGTTGTATAGTGCGCCTACAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                 <  2:1252299/146‑2 (MQ=11)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                                                                                                                                 <  6:1414978/151‑1 (MQ=255)
   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                                                                                                                                 <  1:499067/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGCGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCATGGTAGGTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAAGCCGGCCCCCGCACCCAATAAAAtt                                                                                                                               >  2:1546079/1‑151 (MQ=25)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                 >  2:924112/1‑147 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                <  3:1209788/151‑1 (MQ=255)
    aTCCGTGATCGGGGTAAAAAAAACGTTCCACGAAAACGAGAGCGGAGTAGACTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                <  6:1735317/151‑3 (MQ=18)
     tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaatga                                                                                                                                >  3:1787224/1‑147 (MQ=255)
      gggTGATCGGGGTAAAAAAAAGGTTGCATGAAAACTCGAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                              <  7:923732/149‑7 (MQ=255)
      gcttGATCGTGTTAAAAATAAGGTTGCATGAATTCTCGTTTGGAGTATATTGCGCTTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                               <  1:137638/147‑1 (MQ=14)
      cggggATCGGCGTAAAAAAAAGCGTCCATGAAAAAGCGAGCGGAGTATACTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                               <  2:1812404/146‑1 (MQ=17)
      ccGTGATCGGGGTAAAAAAACGGTTGAATGAAAACGCTAGCGGTGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                                                                                                                               <  1:67490/150‑1 (MQ=18)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                                               <  4:190124/150‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                              <  1:592022/151‑1 (MQ=255)
      ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                              <  8:260680/151‑1 (MQ=255)
       cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCCTCAACGGAAGGTCGTCGGTTCACATCCGGCCCCCGCACCGAATTAAAAttgg                                                                                                                             >  6:1044086/1‑149 (MQ=25)
        gttaTCTGGTTATAAAAAAGGGTGCATGAAAAATCGAGGGCAGTATTCTGCGATTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTtgat                                                                                                                           <  3:206179/148‑1 (MQ=14)
          gATCGGTGTAAAAAAAAGGTTTAATGAAAACTCGATCTGATTATAGTTCGCATCCACGGAAGAGGGGTGGAGCAGCCTGTTAGCTCGTCGGTCTCATAACACGAAGTTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                           <  4:1499713/150‑1 (MQ=18)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                          >  8:296959/1‑141 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgagg                                                                                                                          >  8:843842/1‑149 (MQ=255)
          gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAATTtgtgg                                                                                                                          >  6:262519/1‑148 (MQ=255)
           atgggcttataaatcaGGTTGCATGAAATCTAGAGCGTTGTATTCTGCGCTTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa                                                                                                                        <  5:739432/136‑1 (MQ=11)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                         <  4:1361409/151‑1 (MQ=255)
           aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatga                                                                                                                         <  6:64083/151‑1 (MQ=255)
            tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                         >  3:1575490/1‑139 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatttaacacccta                                                                                                                        >  6:492519/1‑138 (MQ=255)
             cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCCGCCCCCGCAACCAATTattgaag                                                                                                                              >  4:1736678/1‑138 (MQ=255)
              ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGt                                                                                                                      <  1:1421117/151‑1 (MQ=255)
               gggTAAAAATAAGGTTGCATGAAAACGCTAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccctaacg                                                                                                                     <  7:139823/151‑16 (MQ=17)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                                   <  5:1533/151‑1 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                                   <  6:886830/151‑1 (MQ=255)
                ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                                   <  7:1778705/151‑1 (MQ=255)
                 gTAAAAAAAAGGGTCCATCAAAAAGACATCGTAGTAGAGAGCGCATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                 <  3:1902577/151‑2 (MQ=11)
                   ataactaaGGATGCAGGAAAAAGCGTTCGGTGTAGAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaag                                                                                                                <  8:1755479/144‑3 (MQ=11)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                               <  1:1562696/151‑1 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                               <  4:1261898/151‑1 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                               <  6:193149/151‑1 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                               <  3:30693/151‑1 (MQ=255)
                    aaaaaaaGGTTGCATGAAAACGCGAGCGGAGAATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                               <  1:1247342/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                              >  5:855133/1‑130 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                              <  6:973238/151‑1 (MQ=255)
                     aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                              <  6:473738/151‑1 (MQ=255)
                     aaaaaaGGTTGAATGAAAAACATATCGGTGTATATTAGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagcag                                                                                                             <  2:382529/151‑6 (MQ=11)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                             <  8:1463258/151‑1 (MQ=255)
                      aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                             <  3:238575/151‑1 (MQ=255)
                        cccgccgccATGAAATCGCTTGCGCAGTATAGAGCCCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                           <  1:1010816/143‑1 (MQ=11)
                         aaggctcaaggAAAAAGCGAGCGGAGTAGAGTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                          <  7:822381/141‑1 (MQ=11)
                         aaGGTTGCATGAAAACGCTAGCGGAGTATAGTGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacccta                                                                                                                        <  1:1538526/138‑13 (MQ=17)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatagaacacccta                                                                                                                        >  6:1539049/1‑126 (MQ=255)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                          >  3:113423/1‑151 (MQ=255)
                         aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTGAATAAATAAAAGCAGTACg                                                                                                          >  1:1884958/1‑151 (MQ=18)
                         aaGGTTGAATGAAAACGAGATCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                          <  6:1418475/151‑1 (MQ=14)
                           ggtgcaaTGAAAACCCCATCGCAGTATACTCGGCCTCCCCTGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                        <  5:1945808/145‑1 (MQ=11)
                            gtggCATGAAAACGCCTGCGCTGTATTCTCCGCTGCCCCGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                       <  3:843318/148‑1 (MQ=11)
                             gtgaatcaaaaCTAGATCGCAGTATAGTCGGAATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                      <  6:495590/144‑1 (MQ=11)
                              tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                      <  3:590918/150‑1 (MQ=255)
                               gCATGAAAACGCGATCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGATCAAATCCGGCGCCCGCAACCAATTAAAATTTGATGAAGTACAGCAGTACGGTGAcg                                                                                                    <  3:1748495/151‑1 (MQ=255)
                               gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcg                                                                                                    <  5:1163363/151‑1 (MQ=255)
                               ccaggAAAACGCGATGGGAGTATTGTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                     <  6:384058/146‑1 (MQ=11)
                                aaTGAAATCTCTAGCGTAGTATATTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc                                                                                                   <  6:1885740/150‑1 (MQ=1)
                                 aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgcg                                                                                                  <  3:1338398/151‑1 (MQ=255)
                                 aTGAAAACGAGAGCGGAGTAGAGTCGGAATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                  <  8:1902060/151‑1 (MQ=11)
                                   gAAAACGCGAGCGGAGTATAGTGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                          <  5:1349210/126‑11 (MQ=14)
                                   gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGAGGAGGTAAAGCAGTACGGTGAcgcgcg                                                                                                >  2:479657/1‑149 (MQ=255)
                                    aattagagttggtaGTATTGTGGGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                                               <  8:768622/138‑1 (MQ=11)
                                    aaatagcgtttgttgTATAGTTGGCATCCCGGGACGCGGGGTGGAGGAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                                               <  3:800989/137‑1 (MQ=11)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                          <  1:1066829/125‑11 (MQ=255)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaca                                                                                                                             >  6:765482/1‑115 (MQ=255)
                                    aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                 >  1:1644709/1‑115 (MQ=255)
                                    aaaaCGCGAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                 <  6:1645601/118‑4 (MQ=255)
                                     acaCGCGCGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                              <  1:785212/149‑1 (MQ=255)
                                     aatcgagatcGGAGTATAGTGCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt                                                                                                                      <  6:1487916/118‑1 (MQ=14)
                                     aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                              <  8:603271/151‑1 (MQ=255)
                                      aaCGCGAGCGGAGTATAGTGCGAATCCACGGACGCGTGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAATGTCGTCGGTGCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                                                             <  5:40018/151‑1 (MQ=255)
                                         gcgAGCGGAGTATAGTGCGCATCCACGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  7:753468/151‑1 (MQ=255)
                                         gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  2:838368/151‑1 (MQ=255)
                                         gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGag                                                                                          <  6:1345542/151‑1 (MQ=255)
                                           gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                                        <  4:1697201/151‑1 (MQ=255)
                                               gtaggatagggatcATCCCCGGAAGCTGGGTGGAGCAGCCTGGTAGGTCGGCGGTCTCATAACCCGAAGGGCGGCGGTTCAAATCAGTCCCCCGCAACCAATTATATTTTGATGAAGTATCGGAGTACGGGGACGCGGGGTGGAGCAGCCt                                                                                    <  3:902323/138‑1 (MQ=9)
                                               ggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaca                                                                                                                             <  5:1705101/111‑8 (MQ=255)
                                                 aGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                  <  2:12718/151‑1 (MQ=255)
                                                 aGTATAGTGAGCATCCACGGACGCGGGGTGCAGCAGGCTGGTAGCTCGTCGGTCTCATAACCCTAATGTCGGCGGTTCAAATCATCCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACTCTGACGCGGGGTGGAGCAGCCTg                                                                                   <  6:739710/150‑1 (MQ=17)
                                                    atacggataaTCCACGGACTCGGGTGGGATCAGCCTTTTAGCTCTTCGGGATCAGAACCCTAAGGGCGTAGGTTCAAATCCGGCCCCCGCAACCAATAAAACTTTGATGAAGAAACCCAGTACCGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                               <  4:279973/142‑1 (MQ=9)
                                                    ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAGATTTCGGAAGAAAAGCAGTAGGCTGTCGCGtggtggccccgtctgtgcg                                                                               >  6:624958/1‑138 (MQ=17)
                                                    ataGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAg                                                                               <  1:955376/151‑1 (MQ=255)
                                                         gcgcATGGACCGACTCGTGGTGGAGTAGCCTTGTAGCTCTTCGGGATCATATCCCGAATCGCGGAGCTTCAAATACCGCCGGCGCAAGCAAATACAAATAGTAGAAGAAACGGACTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt                                                                          <  4:839812/151‑1 (MQ=2)
                                                         gcgcATCCACGGAAGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTGATATCCGTCCGCCGCAACGAATTAAAATTTGATGAAGTTAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt                                                                          <  1:1598324/151‑1 (MQ=17)
                                                          cgcATCCACGGAAGCGGGTGGGAGCAGCCTGGTAGCTCGTCGGGCTCAAAACCCGAAGGTCGTCGGTGCAAATCAGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtc                                                                         <  8:1539715/151‑1 (MQ=17)
                                                                cACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATaaaat                                                                                                                                 >  7:1260276/1‑89 (MQ=18)
                                                                                           cgatctgcgctataaTAAGCCTAATGGCGTCGGGTCCAATCATCCCCCCGCAACGAATAAAACATTGCTGAATACAAGCAGTACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAAt                                         <  2:643733/136‑1 (MQ=2)
                                                                                                gtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgtcg                                                                        >  7:447576/1‑114 (MQ=255)
                                                                                                             ccggAAGGTCGTAGGTGGATATCCGACCGCTGCAACGAATTATATTTTGATGTATTCTAGCAGTACGGTGACGCGGGGTGGAGCAGCGTGGTAGCTGGTCGGGCTCATAAGCCGAAGGTCGTCGGTTCAAATCCGGCCCGCGCAACGAAt                       <  7:320598/147‑1 (MQ=2)
                                                                                                                   gggcggcgGTTCAACTCCGGCCCCCGCAAACAATTAAACATGGATGAAGTAACGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAttt                <  8:113695/145‑1 (MQ=14)
                                                                                                                                 atccttcccccGCAATCAATTAAAATTTGATGAATTAACGCATTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGc  <  3:1406510/145‑1 (MQ=18)
                                                                                                                                                      |                                                                                                                                  
CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGC  >  NC_000913/2947317‑2947596

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