Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,348,7490AG47.4% 16.2 / 97.3 78intergenic (‑70/+404)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (16/20);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 8.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-02
Rejected as consensus: Frequency below/above cutoff threshold.

TTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCC  >  NC_002947/4348633‑4348898
                                                                                                                      |                                                                                                                                                     
ttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                      <  2:1285738/150‑1 (MQ=17)
ttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                      <  1:178208/150‑1 (MQ=32)
ttCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCCGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTg                                                                                                                      <  1:234578/150‑1 (MQ=17)
    gCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  <  1:427966/150‑1 (MQ=17)
    gCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  <  2:1268206/150‑1 (MQ=32)
      aTTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGACTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                                >  2:669815/1‑150 (MQ=18)
       ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                               <  1:120743/150‑1 (MQ=25)
           ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                           <  1:159295/150‑1 (MQ=25)
           ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                           <  2:430776/150‑1 (MQ=35)
            gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAg                                                                                                          >  1:567147/1‑150 (MQ=35)
                 cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                     >  1:894288/1‑150 (MQ=35)
                 cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTGTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtctc                                                                                                     >  2:184908/1‑150 (MQ=25)
                  aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  <  2:2836/136‑1 (MQ=18)
                  aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  >  1:2836/1‑136 (MQ=18)
                  aaCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                                                                    <  2:1287162/150‑1 (MQ=25)
                   aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTGGCACAGTGGCaa                                                                                                                  >  1:1303079/1‑135 (MQ=18)
                   aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTGGCACAGTGGCaa                                                                                                                  <  2:1303079/135‑1 (MQ=18)
                    cGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  <  2:320312/150‑1 (MQ=25)
                    cGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  <  2:890484/150‑1 (MQ=35)
                     gagaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCTTCGCCTCTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                                                                 <  1:527843/150‑1 (MQ=255)
                      agaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  1:492003/1‑150 (MQ=35)
                           gCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAAt                                                                                           <  2:1203821/150‑1 (MQ=35)
                                gCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                                      <  2:894288/150‑1 (MQ=35)
                                 ccTGCTGGCGGCATAACGTGGTCGGCCCTTGTTGGCGCGGTTAAC‑‑CAATACTTAGTACTGCTTCTGCGGATGTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAATCGTACAGATCTCCCAGGGTAATTcgcgc                                                                                     <  2:264288/150‑1 (MQ=11)
                                  cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                               >  2:843515/1‑123 (MQ=25)
                                  cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCg                                                                                                               <  1:843515/123‑1 (MQ=25)
                                   tACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGa                                                                                   >  1:1297435/1‑150 (MQ=25)
                                      tGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  <  1:438339/116‑1 (MQ=33)
                                      tGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                                  >  2:438339/1‑116 (MQ=33)
                                      tGGCGGCATAACCTCTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCt                                                                                <  1:984787/150‑1 (MQ=25)
                                          ggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCt                                                                            <  1:1078742/150‑1 (MQ=35)
                                            cATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                          >  1:1255541/1‑150 (MQ=35)
                                             aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                                                         >  1:745283/1‑150 (MQ=25)
                                              tAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCtt                                                                        >  2:1238781/1‑150 (MQ=21)
                                               aaCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                                                                       <  2:567147/150‑1 (MQ=35)
                                               aaCCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                                                                       <  2:312594/150‑1 (MQ=34)
                                                       cGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCACCTTTCGACGCGCGGTAGCACAGTGGCAAGCCTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTc                                                               <  1:556606/150‑1 (MQ=255)
                                                          cccTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCGGTGGACGCGCGCTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGa                                                            >  1:452062/1‑150 (MQ=255)
                                                           ccTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGCCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  2:704390/1‑113 (MQ=255)
                                                           ccTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGCCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                <  1:704390/113‑1 (MQ=255)
                                                            cTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATc                                                          <  2:1255541/150‑1 (MQ=35)
                                                               ggtggCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                                       <  1:669815/150‑1 (MQ=35)
                                                                  gtggcgGTTAACCGCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACg                                                      >  1:470792/4‑150 (MQ=255)
                                                                   gcgcGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  <  2:888235/103‑1 (MQ=255)
                                                                   gcgcGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  >  1:888235/1‑103 (MQ=255)
                                                                     gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGt                                                                >  1:291541/1‑135 (MQ=37)
                                                                     gcgGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGt                                                                <  2:291541/135‑1 (MQ=21)
                                                                           caC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                             <  2:329651/115‑1 (MQ=38)
                                                                           aaC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCCCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                             >  1:329651/1‑116 (MQ=37)
                                                                            aC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcc                                          >  2:54440/1‑150 (MQ=25)
                                                                             c‑‑cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcct                                         >  2:1146834/1‑149 (MQ=25)
                                                                                cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGt                                        >  2:765788/1‑150 (MQ=25)
                                                                                 gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTg                                       >  1:766531/1‑150 (MQ=35)
                                                                                  ccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                          >  2:56141/1‑114 (MQ=255)
                                                                                  ccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                          <  1:56141/114‑1 (MQ=255)
                                                                                  ccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                                                         >  1:1338468/1‑115 (MQ=38)
                                                                                  ccACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                                                         <  2:1338468/115‑1 (MQ=38)
                                                                                    aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                    <  1:1238781/150‑1 (MQ=25)
                                                                                    aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCaa                                    <  2:766531/150‑1 (MQ=35)
                                                                                      ttACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                                                           >  2:178652/1‑125 (MQ=37)
                                                                                      ttACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                                                           <  1:178652/125‑1 (MQ=35)
                                                                                          tacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTg                              <  1:396190/150‑1 (MQ=25)
                                                                                          tacGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTg                              >  2:378292/1‑150 (MQ=35)
                                                                                           acGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTgg                             >  1:300841/1‑150 (MQ=21)
                                                                                            cGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTggg                            <  1:765788/150‑1 (MQ=25)
                                                                                            cGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  <  2:189631/80‑1 (MQ=255)
                                                                                            cGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                                  >  1:189631/1‑80 (MQ=255)
                                                                                                 tCGGCTGACTTCTGGACGCTCCTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCACGGTAATTCGTGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTg                       <  1:410329/150‑1 (MQ=255)
                                                                                                    gCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgc                                                                                       <  1:1243201/83‑1 (MQ=40)
                                                                                                    gCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgc                                                                                       >  2:1243201/1‑83 (MQ=39)
                                                                                                        aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                >  2:42765/1‑70 (MQ=255)
                                                                                                        aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                                <  1:42765/70‑1 (MQ=255)
                                                                                                        aCTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTTGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAttt                >  1:1369220/1‑150 (MQ=21)
                                                                                                          ttCTGCACGCTCCTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt              <  2:24305/150‑1 (MQ=25)
                                                                                                              gCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCa          >  1:861071/1‑150 (MQ=25)
                                                                                                               cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCAc         >  2:146890/1‑150 (MQ=25)
                                                                                                                    ttaTCCCTTCTTCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAc    <  2:78119/148‑1 (MQ=255)
                                                                                                                    tCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca                                                 <  1:1264941/105‑1 (MQ=38)
                                                                                                                    tCATCCCGTCGCCTTTCGACGAGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca                                                 >  2:1264941/1‑105 (MQ=255)
                                                                                                                      aTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                                                                    <  2:861358/84‑1 (MQ=40)
                                                                                                                      aTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                                                                    >  1:861358/1‑84 (MQ=40)
                                                                                                                      aTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCAccc  >  2:273612/1‑150 (MQ=35)
                                                                                                                      |                                                                                                                                                     
TTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAAC‑‑CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCACCTCACCC  >  NC_002947/4348633‑4348898

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: