Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,156,185 A→T 8.3% L217* (TTA→TAA)  PP_3659 ← transcriptional regulator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,156,1850AT 8.3% 287.9 / 11.6 93L217* (TTA→TAA) PP_3659transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base A (44/41);  new base T (5/3);  total (49/44)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGC  >  NC_002947/4156037‑4156334
                                                                                                                                                    |                                                                                                                                                     
ggTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaa                                                                                                                                                      >  2:2029160/1‑150 (MQ=255)
 gTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaac                                                                                                                                                     >  2:696122/1‑149 (MQ=255)
 gTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaac                                                                                                                                                     <  2:2746314/150‑2 (MQ=255)
      aCTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGc                                                                                                                                                <  1:998010/150‑1 (MQ=255)
      aCTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGAGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGc                                                                                                                                                >  1:1309543/1‑150 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTATTTTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCa                                                                                                                                               <  2:2720475/150‑1 (MQ=255)
       cTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAACCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCa                                                                                                                                               <  1:1062144/150‑1 (MQ=255)
           aCATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGGGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGaa                                                                                                                                           <  1:1540222/150‑1 (MQ=255)
            cATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAag                                                                                                                                          <  2:833585/150‑1 (MQ=255)
             aTGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAaga                                                                                                                                         <  1:2029160/150‑1 (MQ=255)
               gCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAgaga                                                                                                                                       >  2:814156/1‑150 (MQ=255)
                 tCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGaaa                                                                                                                                     >  1:1205369/1‑150 (MQ=255)
                    aGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACgg                                                                                                                                  <  1:1602389/150‑1 (MQ=255)
                      tCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGc                                                                                                                                >  2:2559025/1‑150 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            <  1:761078/150‑1 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            <  1:1116146/150‑1 (MQ=255)
                          gggCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            <  1:2816654/150‑1 (MQ=255)
                            gCGAGATCATCGAGTTCAGGCCACGTGTGTGAAGAGAGCAGTTCTAAGCGTCGCCGCATCATGATTGTCCTACTTGCGGCGTTACGGTATTTAACAAGGAAGTTACGTGGCGTACCATTTAAAAAGGCAGGAAGAGAAACGGGCGTTTgg                                                                                                                          <  2:2867867/150‑1 (MQ=255)
                                gCTCATCGAGTGCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGTCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCaaa                                                                                                                      <  2:1309543/150‑1 (MQ=255)
                                 cTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAAc                                                                                                                     <  2:414532/150‑1 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                  >  2:2852323/1‑150 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCtt                                                                                                                  >  1:1212458/1‑150 (MQ=255)
                                    aTCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACCGGCGTTTGGCAAACCtt                                                                                                                  >  2:2725555/1‑150 (MQ=255)
                                     tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCg                                                                                                                               >  2:2029868/1‑136 (MQ=255)
                                     tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCg                                                                                                                               <  1:2029868/136‑1 (MQ=255)
                                     tCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTg                                                                                                                 >  2:1726008/1‑150 (MQ=255)
                                            cAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCa                                                                                                          <  2:2794100/150‑1 (MQ=255)
                                               gCCACGTGTGTGGAGAGAGTAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGaa                                                                                                       >  1:2358684/1‑150 (MQ=255)
                                               gCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            >  1:916726/1‑129 (MQ=255)
                                                  aCGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            <  2:1067549/126‑1 (MQ=255)
                                                  aCGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGttt                                                                                                                            >  1:1067549/1‑126 (MQ=255)
                                                  aCGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGt                                                                                                                              >  2:1019013/1‑124 (MQ=255)
                                                  aCGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGt                                                                                                                              <  1:1019013/124‑1 (MQ=255)
                                                       tgtgGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTAATGTCTGAGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCtcttc                                                                                               >  1:2382583/1‑150 (MQ=255)
                                                         tgGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTtctc                                                                                             <  2:1205369/150‑1 (MQ=255)
                                                           gagagaGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGa                                                                                           >  1:1408177/1‑150 (MQ=255)
                                                           gagagaGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGa                                                                                           >  1:976064/1‑150 (MQ=255)
                                                           gagagaGCAGTTTTAATAGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGa                                                                                           >  2:2932846/1‑150 (MQ=255)
                                                                agCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAgaga                                                                                                                                       <  1:1020737/101‑1 (MQ=255)
                                                                agCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAgaga                                                                                                                                       >  2:1020737/1‑101 (MQ=255)
                                                                   aGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCg                                                                                   >  1:732714/1‑150 (MQ=255)
                                                                     ttttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAAGCCAGAAAGTGCTTCTCGAACACAGCGGa                                                                                 >  2:1478729/1‑150 (MQ=255)
                                                                       ttAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAgg                                                                               <  1:2868924/150‑1 (MQ=255)
                                                                         aaTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCacg                                                                                                                                             >  1:238036/1‑84 (MQ=255)
                                                                         aaTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTTAAAAGGCacg                                                                                                                                             <  2:238036/86‑3 (MQ=255)
                                                                          aTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGc                                                                                                            >  2:2855971/1‑118 (MQ=255)
                                                                          aTCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGc                                                                                                            <  1:2855971/118‑1 (MQ=255)
                                                                              ccgccgGATCATAATGCTCGTACTGTCCTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAc                                                                                                              <  1:2466665/112‑1 (MQ=255)
                                                                              ccgccgGATCATAATGCTCGTACTGTCCTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTAc                                                                                                              >  2:2466665/1‑112 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                   >  2:1319852/1‑64 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                   >  2:1319796/1‑64 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                   <  1:1319852/64‑1 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTaaaaa                                                                                                                                                   <  1:1319796/64‑1 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTa                                                                                                                                                       >  2:1172451/1‑60 (MQ=255)
                                                                                         tAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTa                                                                                                                                                       <  1:1172451/60‑1 (MQ=255)
                                                                                            tGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGa                                                          >  2:2221244/1‑150 (MQ=255)
                                                                                             gCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCCGGCGAAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAt                                                         <  1:1726008/150‑1 (MQ=255)
                                                                                             gCTCGTACTGTCGTCGTTACGGTATTTAACAAGGAGATTGACTGCAGTTGCTCTTTACAGGACAGGAAGAGAAAGGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGAt                                                         >  1:837698/1‑150 (MQ=255)
                                                                                                 gTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTAGCtttac                                                     >  1:517764/1‑148 (MQ=255)
                                                                                                   aCTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTc                                                   >  1:2274349/1‑150 (MQ=255)
                                                                                                    cTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                         >  2:438323/1‑127 (MQ=255)
                                                                                                    cTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTAt                                                                         <  1:438323/127‑1 (MQ=255)
                                                                                                    cTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTca                                                  >  2:2766454/1‑150 (MQ=255)
                                                                                                         gtcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACggg                                                                                                                                 >  2:1666591/1‑66 (MQ=255)
                                                                                                         gtcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACggg                                                                                                                                 <  1:1666591/66‑1 (MQ=255)
                                                                                                          tcgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAg                                            >  1:230437/1‑150 (MQ=255)
                                                                                                           cgtTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgc                                           <  1:814156/150‑1 (MQ=255)
                                                                                                             ttACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAgcgc                                         >  2:125520/1‑150 (MQ=255)
                                                                                                                 ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAaga                                                                                                                                         >  1:811033/1‑50 (MQ=255)
                                                                                                                 ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAaga                                                                                                                                         <  2:811033/50‑1 (MQ=255)
                                                                                                                 ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGc                                                                                                                                <  1:1020446/59‑1 (MQ=255)
                                                                                                                 ggTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGc                                                                                                                                >  2:1020446/1‑59 (MQ=255)
                                                                                                                     tttCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATa                                 >  1:888188/1‑150 (MQ=255)
                                                                                                                         acaAGGAAGTTACTTTGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTc                                                   >  2:1416649/1‑128 (MQ=255)
                                                                                                                         acaAGGAAGTTACTTTGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTc                                                   <  1:1416649/128‑1 (MQ=255)
                                                                                                                           aaGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCt                           >  2:2980644/1‑150 (MQ=255)
                                                                                                                           aaGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCt                           >  2:2265849/1‑150 (MQ=255)
                                                                                                                              gAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTgtc                        <  2:1212458/150‑1 (MQ=255)
                                                                                                                                 gTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAg                                                                                >  2:2476407/1‑91 (MQ=255)
                                                                                                                                  ttACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCt                                                    <  1:1180941/118‑1 (MQ=255)
                                                                                                                                  ttACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCt                                                    >  2:1180941/1‑118 (MQ=255)
                                                                                                                                        ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGc              >  1:293753/5‑150 (MQ=255)
                                                                                                                                        ttcctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCACTTCTCGAACACAGCGGAGGTTTTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCTATAGTCGCTGTCGTCCAccaga              >  2:646026/5‑147 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCACCACTTGGTTGCCAGAGAGCTCTGTTCGAACCCAGAGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga           <  1:691493/150‑1 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAAcaca                                                                                      <  2:2073100/75‑1 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAAcaca                                                                                      >  1:2073100/1‑75 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga           <  1:1478729/150‑1 (MQ=255)
                                                                                                                                           gTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACAATGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga           <  1:2637066/150‑1 (MQ=255)
                                                                                                                                           ctGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTga           <  2:293753/149‑1 (MQ=255)
                                                                                                                                            tGCCTTTTTAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTgag          >  1:1613094/1‑150 (MQ=255)
                                                                                                                                             gCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGc                                                                     >  2:1136186/1‑90 (MQ=255)
                                                                                                                                             accTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGc                                                                     <  1:1136186/89‑1 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCtcttc                                                                                               <  1:1504815/61‑1 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCtcttc                                                                                               >  2:1504815/1‑61 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg      >  1:2472076/1‑150 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg      >  1:2603901/1‑150 (MQ=255)
                                                                                                                                                ttttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGaggagg      <  1:696122/150‑1 (MQ=255)
                                                                                                                                                 tttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAggagga     >  1:1886045/1‑150 (MQ=255)
                                                                                                                                                 tttAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAggagga     <  1:2852323/150‑1 (MQ=255)
                                                                                                                                                    aaaaaGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGc  >  2:583851/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAGGAAGTTACTTGGCGTGCCTTTTAAAAAGGCAGGAAGAGAAACGGGCGTTTGGCAAACCTTGTACGCCAGAAAGCTCTTCTCGAACACAGCGGAGGTTCTATCGGCCCTCCACTCGATTTCCTCACAAAAGCGCATTCGATAGTCGCTGTCGTCCACCGGCTGAGGAGGACGC  >  NC_002947/4156037‑4156334

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: