Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,348,7490AG53.3% 6.4 / 169.5 135intergenic (‑70/+404)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (33/30);  new base G (34/38);  total (67/68)
Fisher's exact test for biased strand distribution p-value = 6.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 135 total aligned reads displayed.

TTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAC  >  NC_002947/4348632‑4348896
                                                                                                                      |                                                                                                                                                   
tttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAg                                                                                                                      >  2:444950/1‑149 (MQ=11)
tttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                     <  2:482898/150‑1 (MQ=14)
 ttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGCTCGGACCTAAATGGGGGTGTTTCCCGCCCGTACGTATGGCATCGGCTGTCTTCTGCTCGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                    <  2:1968834/150‑1 (MQ=2)
    tGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCa                                                                                                                 >  2:1420638/1‑150 (MQ=32)
     gCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                <  2:287530/150‑1 (MQ=17)
     gCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                <  2:1693319/150‑1 (MQ=17)
     gCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCCCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGGAGCACAGTGGCaa                                                                                                                >  2:438865/1‑150 (MQ=12)
        ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                             >  2:1659705/1‑150 (MQ=35)
          aCGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                           <  2:2097634/150‑1 (MQ=25)
          aCGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                           <  2:620671/150‑1 (MQ=35)
                aCCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtc                                                                                                     >  2:2059527/1‑150 (MQ=25)
                  cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                   >  1:635924/1‑150 (MQ=25)
                  cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                   >  2:1221511/1‑150 (MQ=35)
                  cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                   >  1:948793/1‑150 (MQ=35)
                     cGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                <  2:1727976/150‑1 (MQ=25)
                       agaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAGCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                              <  2:990037/150‑1 (MQ=255)
                         aCGGCGACGCCCCCTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGt                                                                                            <  2:417431/150‑1 (MQ=255)
                                 gCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                          >  2:2150702/1‑128 (MQ=37)
                                 gCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                          <  1:2150702/128‑1 (MQ=37)
                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                          <  2:1516697/126‑1 (MQ=35)
                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                          >  1:1516697/1‑126 (MQ=35)
                                   cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgcg                                                                                  >  1:1133161/1‑150 (MQ=35)
                                    tACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGa                                                                                 >  2:2038910/1‑150 (MQ=35)
                                     aCTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                                                                <  1:1060962/150‑1 (MQ=35)
                                       tGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                   >  1:640261/1‑145 (MQ=25)
                                       tGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                   <  2:640261/145‑1 (MQ=25)
                                         gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTc                                                                            >  2:868056/1‑150 (MQ=25)
                                         gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTc                                                                            >  1:2155232/1‑150 (MQ=25)
                                         gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTc                                                                            >  1:556305/1‑150 (MQ=35)
                                           ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCt                                                                          <  1:298095/150‑1 (MQ=25)
                                           ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCACGCGTa                                                                                                           <  2:215925/117‑1 (MQ=25)
                                           ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                           >  1:215925/1‑117 (MQ=38)
                                            gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                              >  1:2065936/1‑129 (MQ=37)
                                             cATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                        <  1:71152/150‑1 (MQ=25)
                                              aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                                                       >  2:890116/1‑150 (MQ=25)
                                               tAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                        <  2:2147298/148‑1 (MQ=35)
                                                 aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                    <  2:635924/150‑1 (MQ=25)
                                                 aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                                    >  2:383020/1‑150 (MQ=35)
                                                    tGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                 <  1:2088399/150‑1 (MQ=25)
                                                    tGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                                 >  2:2010342/1‑150 (MQ=34)
                                                     gTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCt                                                                <  2:556305/150‑1 (MQ=35)
                                                        cGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGta                                                             <  1:2038910/150‑2 (MQ=25)
                                                         ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                           <  2:477508/135‑1 (MQ=37)
                                                         ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                           >  1:477508/1‑135 (MQ=34)
                                                          gCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCgg                                                           >  2:1102293/1‑150 (MQ=25)
                                                          gCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCgg                                                           >  1:1459783/1‑150 (MQ=25)
                                                          gCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCgg                                                           >  1:1411943/1‑150 (MQ=34)
                                                           cccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                                                                <  1:749897/128‑1 (MQ=37)
                                                           cccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                                                                >  2:749897/1‑128 (MQ=37)
                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATg                                                                   >  1:2039348/1‑140 (MQ=37)
                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGATTATg                                                                   <  2:2039348/140‑1 (MQ=255)
                                                            ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                                                         >  2:669094/1‑150 (MQ=34)
                                                             cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATc                                                        <  1:1221511/150‑1 (MQ=35)
                                                             cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATc                                                        >  2:1183909/1‑150 (MQ=35)
                                                             cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATc                                                        <  1:445554/150‑1 (MQ=35)
                                                             cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGACCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGGCGGATc                                                        >  1:153660/1‑150 (MQ=21)
                                                                ggtggCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                                     <  2:1133161/150‑1 (MQ=35)
                                                                    gcgcGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca                                                 >  1:1927315/1‑150 (MQ=35)
                                                                    gcgcGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca                                                 >  1:3869/1‑150 (MQ=35)
                                                                      gcgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                              >  1:629030/1‑87 (MQ=255)
                                                                      gcgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                         >  2:90580/1‑92 (MQ=255)
                                                                      gcgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                         <  1:90580/92‑1 (MQ=255)
                                                                          ttAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGc                                                                                                                  >  2:1876467/1‑79 (MQ=255)
                                                                          ttAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGc                                                                                                                  <  1:1876467/79‑1 (MQ=255)
                                                                          ttAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGt                                           >  2:1834065/1‑150 (MQ=25)
                                                                            aaCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTc                                         <  2:948793/150‑1 (MQ=35)
                                                                             aCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcc                                        <  1:1659705/150‑1 (MQ=35)
                                                                               cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCTCGACCTTCCTGCTTATGc                                                                  >  2:1694430/1‑122 (MQ=255)
                                                                                gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTg                                     >  1:306144/1‑150 (MQ=25)
                                                                                gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTg                                     <  1:1834065/150‑1 (MQ=25)
                                                                                   aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                  <  1:354641/150‑1 (MQ=25)
                                                                                   aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                  <  2:122796/150‑1 (MQ=35)
                                                                                               cTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                  >  2:775499/1‑139 (MQ=25)
                                                                                               cTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                  <  1:775499/139‑1 (MQ=37)
                                                                                               cTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTCTTGGGCtt                       >  2:623059/1‑150 (MQ=25)
                                                                                               cTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCtt                       >  2:623218/1‑150 (MQ=25)
                                                                                                ttCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                    >  1:1922425/1‑56 (MQ=35)
                                                                                                ttCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                    <  2:1922425/56‑1 (MQ=35)
                                                                                                ttCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                         >  1:1700302/1‑67 (MQ=255)
                                                                                                ttCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCttt                      >  2:1724150/1‑150 (MQ=35)
                                                                                                ttCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCttt                      >  2:1025978/1‑150 (MQ=35)
                                                                                                     cTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTACTTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGa                 >  1:2138240/1‑150 (MQ=25)
                                                                                                     cTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGa                 >  1:847436/1‑150 (MQ=35)
                                                                                                       gACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAtt               <  2:368772/150‑1 (MQ=25)
                                                                                                        aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTg                            <  2:1643842/136‑1 (MQ=25)
                                                                                                        aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTg                            >  1:1643842/1‑136 (MQ=25)
                                                                                                        aCTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                                                               >  2:1813327/1‑69 (MQ=255)
                                                                                                        aCTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                                                               <  1:1813327/69‑1 (MQ=255)
                                                                                                         cTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg             <  2:1724283/150‑1 (MQ=35)
                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                     >  1:966888/1‑45 (MQ=255)
                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGt                                                                                                                     <  2:966888/45‑1 (MQ=255)
                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt            >  1:394459/1‑150 (MQ=25)
                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt            <  1:868056/150‑1 (MQ=25)
                                                                                                          ttCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt            <  2:1459783/150‑1 (MQ=25)
                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt            >  2:1657749/1‑150 (MQ=35)
                                                                                                          ttCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt            <  1:1687101/150‑1 (MQ=35)
                                                                                                               cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACTTATATGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAc       <  1:1302002/150‑1 (MQ=11)
                                                                                                                aCGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGCGCTTTGGCGATTTGTGCCAcc      >  1:1257988/1‑150 (MQ=255)
                                                                                                                  gCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTc    >  2:216022/1‑150 (MQ=25)
                                                                                                                    tCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAc  <  1:384313/150‑1 (MQ=35)
                                                                                                                      |                                                                                                                                                   
TTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAC  >  NC_002947/4348632‑4348896

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