Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,481,791 C→T 5.7% P392P (CCC→CCT ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,481,7910CT 5.7% 496.6 / 11.5 158P392P (CCC→CCTybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (77/72);  new base T (7/2);  total (84/74)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 163 total aligned reads displayed.

GCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGA  >  NC_002947/4481644‑4481940
                                                                                                                                                   |                                                                                                                                                     
  tgGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCa                                                                                                                                                   >  2:1466231/1‑150 (MQ=255)
  tgGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCa                                                                                                                                                   >  2:4733403/1‑150 (MQ=255)
    gAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCAcc                                                                                                                                                 >  2:697645/1‑150 (MQ=255)
      gCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCtt                                                                                                                                               >  2:1990661/1‑150 (MQ=255)
      gCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCtt                                                                                                                                               >  2:3096105/1‑150 (MQ=255)
      gCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCtt                                                                                                                                               >  1:1872418/1‑150 (MQ=255)
          aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGc                                                                                                                                           >  2:2092194/1‑150 (MQ=255)
          aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGc                                                                                                                                           >  2:4665740/1‑150 (MQ=255)
          aTCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGc                                                                                                                                           >  2:3181344/1‑150 (MQ=255)
              aaGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTg                                                                                                                                       <  2:1689087/150‑1 (MQ=255)
               aGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGa                                                                                                                                      <  1:3617923/150‑1 (MQ=255)
                      cgcCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACctc                                                                                                                               <  2:990250/150‑1 (MQ=255)
                       gcCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACctct                                                                                                                              >  1:1362091/1‑150 (MQ=255)
                          tGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTaga                                                                                                                           <  1:3502358/150‑1 (MQ=255)
                             aaCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGTAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgagga                                                                                                                        <  1:3024810/150‑1 (MQ=255)
                             aaCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgagga                                                                                                                        >  1:4897546/1‑150 (MQ=255)
                             aaCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgagga                                                                                                                        <  1:2226498/150‑1 (MQ=255)
                             aaCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgagga                                                                                                                        <  1:3332456/150‑1 (MQ=255)
                              aCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAgaggag                                                                                                                       >  2:4296090/1‑150 (MQ=255)
                              aCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTAACCACCTCTAgagga                                                                                                                        >  1:4865162/1‑149 (MQ=255)
                              aCCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTAACCACCTCTAgagga                                                                                                                        <  2:4865162/149‑1 (MQ=255)
                                 ggAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCCAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGaa                                                                                                                    <  1:3443537/150‑1 (MQ=255)
                                     cTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGAt                                                                                                                <  2:909832/150‑1 (MQ=255)
                                        aTCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCAt                                                                                                             <  2:4421186/150‑1 (MQ=255)
                                            ttACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAAc                                                                                                         <  2:1138420/150‑1 (MQ=255)
                                                tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtc                                                                                                     >  1:2100497/1‑150 (MQ=255)
                                                tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtc                                                                                                     >  1:347486/1‑150 (MQ=255)
                                                tCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtc                                                                                                     >  1:2736572/1‑150 (MQ=255)
                                                  aCCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtctc                                                                                                   >  1:3069386/1‑150 (MQ=255)
                                                  aCCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGtctc                                                                                                   >  2:3854417/1‑150 (MQ=255)
                                                      cTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGAACAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCa                                                                                               <  2:1980164/150‑1 (MQ=255)
                                                      cTTAGCCCGAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCa                                                                                               <  2:4529118/150‑1 (MQ=255)
                                                       ttAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAg                                                                                              >  2:1846920/1‑150 (MQ=255)
                                                            ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGCCCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACt                                                                                         >  2:1584694/1‑150 (MQ=255)
                                                            ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCCGTTACt                                                                                         >  1:4539715/1‑150 (MQ=255)
                                                            ccTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACt                                                                                         >  2:1584771/1‑150 (MQ=255)
                                                             cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                        >  1:1007692/1‑150 (MQ=255)
                                                             cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                        >  2:4286882/1‑150 (MQ=255)
                                                             cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                        >  1:1007682/1‑150 (MQ=255)
                                                             cTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACtt                                                                                        >  1:4065102/1‑150 (MQ=255)
                                                              tAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTc                                                                                       >  2:4381300/1‑150 (MQ=255)
                                                              tAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCGTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTc                                                                                       >  1:595556/1‑150 (MQ=255)
                                                                gaagaGGCGGCCATGGGCTACAAAAGGTTCGATCAACTGCAGGATACCTGCCGCACGGTCATCCTGGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTg                                                                                     <  2:4746429/150‑1 (MQ=255)
                                                                      gcggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGc                                                                               >  2:3949647/1‑150 (MQ=255)
                                                                      gcggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGc                                                                               >  1:2154027/1‑150 (MQ=255)
                                                                        ggcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCg                                                                             >  1:956013/1‑150 (MQ=255)
                                                                         gcCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGACTGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGt                                                                            >  1:4178447/1‑150 (MQ=255)
                                                                           cATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAg                                                                          >  2:2425946/1‑150 (MQ=255)
                                                                              gggCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCgg                                                                       <  2:1123654/150‑1 (MQ=255)
                                                                                 cTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGACTGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCaa                                                                    >  1:4864245/1‑150 (MQ=255)
                                                                                 cTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCaa                                                                    <  1:904547/150‑1 (MQ=255)
                                                                                 cTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCaa                                                                    <  2:1660519/150‑1 (MQ=255)
                                                                                  tACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAt                                                                   >  2:228769/1‑150 (MQ=255)
                                                                                     aaaaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttct                                                                >  1:1820442/1‑150 (MQ=255)
                                                                                     aaaaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGCTACGGCAAtcca                                                                <  1:3983039/150‑4 (MQ=255)
                                                                                      aaaTGGTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttctt                                                               <  1:657385/150‑1 (MQ=255)
                                                                                       aaTGTTCGATCAACAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCGAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                              <  2:566924/150‑1 (MQ=255)
                                                                                       aaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                              <  1:2194350/150‑1 (MQ=255)
                                                                                       aaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAAttcttc                                                              <  1:4285126/150‑1 (MQ=255)
                                                                                       aaTGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGACGTAGCGGCAAttctta                                                              >  1:258731/1‑149 (MQ=255)
                                                                                          gTTCGCTCAAACGCAGGATATCTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGTGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGc                                                           <  1:958849/150‑1 (MQ=255)
                                                                                              gATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCACTTCTTCTGCTCAt                                                       >  2:4139564/1‑150 (MQ=255)
                                                                                              gATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCAt                                                       >  2:4140087/1‑150 (MQ=255)
                                                                                               aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATc                                                      <  2:2289504/150‑1 (MQ=255)
                                                                                               aTCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATc                                                      >  1:213558/1‑150 (MQ=255)
                                                                                                  aaaaGCAGGATAACTGCCGCCAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                   <  2:1371539/150‑1 (MQ=255)
                                                                                                  aaaaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGAc                                                   <  2:481129/150‑1 (MQ=255)
                                                                                                    aaGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCt                                                 <  2:387487/150‑1 (MQ=255)
                                                                                                        aGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTACTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTgg                                             >  1:783579/1‑150 (MQ=255)
                                                                                                             aaCTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCAt                                        <  2:3864813/150‑1 (MQ=255)
                                                                                                              aCTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATc                                       <  1:915851/150‑1 (MQ=255)
                                                                                                                    gCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGc                                 <  1:4944740/150‑1 (MQ=255)
                                                                                                                       cggTCATCCTCGTGCCTGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTTGGCCATCGTCAGCgtg                              <  1:3965921/149‑1 (MQ=255)
                                                                                                                          tCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTc                           <  1:1059945/150‑1 (MQ=255)
                                                                                                                            aTCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACATCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                         >  1:4073190/1‑150 (MQ=255)
                                                                                                                            aTCCTCGTGCCAGGTGCTACACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCg                         <  2:4182535/150‑1 (MQ=255)
                                                                                                                             tCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCCTCGTCAGCGTGTTCCgg                        >  2:1274147/1‑150 (MQ=255)
                                                                                                                             ccaaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTcc                          >  2:1106230/6‑148 (MQ=255)
                                                                                                                             ccaaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTcc                          <  1:1106230/143‑1 (MQ=255)
                                                                                                                              ccTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCCTCGACCTGTGGGCCATCGTCAGCGTGTTCCGGa                       <  2:879097/150‑1 (MQ=255)
                                                                                                                              ccTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGa                       <  1:633152/150‑1 (MQ=255)
                                                                                                                              caaggTGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGa                       >  2:4436765/5‑150 (MQ=255)
                                                                                                                                tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCCCCTCGCGAGGAGGAATGATCAGGAACGGTCTCGCCAGTGCCTGCTGGATTGCCGTAGCGGCACTTCTTCTGCTCATCGACCTGTGGGCCATCGGCAGCGTGTTCCGGAGc                     >  2:2820696/1‑150 (MQ=255)
                                                                                                                                tCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGc                     >  2:3207278/1‑150 (MQ=255)
                                                                                                                                   tGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGAc                  >  2:4431305/1‑150 (MQ=255)
                                                                                                                                   tGCCGGGTGCAGCACCTGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGAc                  >  1:371665/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACaa                >  2:2862982/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACaa                >  2:3312363/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACaa                >  2:501792/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACaa                >  2:487359/1‑150 (MQ=255)
                                                                                                                                     ccAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCGAGAGGAGGACTGATCATGAACGGTCTCGCCCGGTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCAGCGTCAGCGTGTTCCGGCGCGACaa                >  2:3312601/1‑150 (MQ=255)
                                                                                                                                       aGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGa              <  1:143146/150‑1 (MQ=255)
                                                                                                                                           gctgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGa          <  1:3831444/150‑1 (MQ=255)
                                                                                                                                            ctgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGACCATCGTCAGCGTGTTCCGGAGCGACAAGACCGAt         >  2:370332/1‑150 (MQ=255)
                                                                                                                                             tgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATg        >  2:1178564/1‑150 (MQ=255)
                                                                                                                                             tgcACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATg        <  2:3618360/150‑1 (MQ=255)
                                                                                                                                                aCCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCg     <  1:1198532/150‑1 (MQ=255)
                                                                                                                                                  ccGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGACTGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTg   >  1:1587012/1‑150 (MQ=255)
                                                                                                                                                  ccGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTg   >  1:3171691/1‑150 (MQ=255)
                                                                                                                                                   cGGCACCTTGGGCCCTGACCACCTGTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGa  <  1:4267987/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGGAGCACATCGAAGCAGGGCGCCTGCAACCGGAACTGATCGTTACTCACCGCTTAGCCCTAGAAGAGGCGGCCATGGGCTACAAAATGTTCGATCAAAAGCAGGATAACTGCCGCAAGGTCATCCTCGTGCCAGGTGCTGCACCCGGCACCTTGGGCCCTGACCACCTCTAGAGGAGGAATGATCATGAACGGTCTCGCCAGTTACTTCTGGATTGCCGTAGCGGCAATTCTTCTGCTCATCGACCTGTGGGCCATCGTCAGCGTGTTCCGGAGCGACAAGACCGATGCCGTGA  >  NC_002947/4481644‑4481940

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: