Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,5910A.35.9% 218.8 / 126.1 92intergenic (+18/‑16)metW/metVtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (55/4);  new base . (30/3);  total (85/7)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.20e-01

CCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947453‑2947740
                                                                                                                                            |                                                                                                                                                     
ccgcggaaCCAATTACATTTTGTTGAATTCTCGGATTCCTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                                           <  5:607190/145‑1 (MQ=14)
          aaTTAAAATTTGTTGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATTACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCGGTACGGTGACGCgggg                                                                                                                                 <  7:1378446/151‑1 (MQ=255)
                                                   tGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag                                                                                                                            >  8:284830/1‑115 (MQ=255)
                                                   tGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGag                                                                                                                            <  7:285281/115‑1 (MQ=255)
                                                      agcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                                                                                     >  6:620934/1‑151 (MQ=39)
                                                                     cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTTATGAAGTaaaa                                                                                                                                                    >  4:351021/1‑73 (MQ=255)
                                                                     cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                                                                    <  3:351047/73‑1 (MQ=255)
                                                                                                     aaaTCCGGCCCCCGCAACCAATGAATTATCATCAAGTAACAGCTGTAGGTTGAAGCGGGCTCGAGTAGCCACGAAGCACGTCGCGCACTTATAACGAAGGGCTGAGGGGCACCTCCGGCCCACGCAAACAATCATTGAACAcccaacaggc                                      >  7:752107/1‑144 (MQ=255)
                                                                                                          cGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGGGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTGGTCGGTGCAACTCCGGCCCCCGCAACCAATTTTTGAACACCCTAACGGGAGttt                                 >  1:354373/1‑151 (MQ=17)
                                                                                                            gCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttt                               >  2:37725/1‑151 (MQ=255)
                                                                                                            gCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttt                               >  7:508197/1‑151 (MQ=255)
                                                                                                            gCCCCCGCAACCAATCAAATTTCATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGAAGCTCGTCGGGCTCATAACGCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGCTGtttgt                               >  3:58370/1‑149 (MQ=255)
                                                                                                             cccccGCAACCAATCAAATTTGATGAAGTAAACGCAGTACGGTGACGAGGGGTGGAGCAGACTGGTAGCTCGTCGGGCTCATAACCCGAAGGTAGTCGGTTCAAATCCCGCCCCCGCAACCAATTATTGAACACCCTACCGGGTGtttttt                              >  1:1179331/1‑151 (MQ=18)
                                                                                                              ccccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  4:1319340/1‑151 (MQ=255)
                                                                                                               cccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTg                            >  2:1615366/1‑151 (MQ=255)
                                                                                                                ccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGt                           >  8:814599/1‑151 (MQ=255)
                                                                                                                ccGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACCGGTGTTTTTTTGt                           >  6:343414/1‑151 (MQ=255)
                                                                                                                 cGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACACGCAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAATGAACACGCTTTCTGGTGATTATTTggt                          >  7:1474665/1‑149 (MQ=11)
                                                                                                                  ggctctgccgatctaTTTCATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGtt                          <  4:1225084/137‑1 (MQ=17)
                                                                                                                  gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             <  5:842222/147‑1 (MQ=255)
                                                                                                                  gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  6:843030/1‑147 (MQ=255)
                                                                                                                  gCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGttt                         >  5:211253/1‑151 (MQ=255)
                                                                                                                   cAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTATAGGCGGTGTTTTTGTTTc                        >  1:1508017/1‑151 (MQ=11)
                                                                                                                   cAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTTATCACCCAATAGCGTGatttttgggatc                        >  6:1356382/1‑145 (MQ=11)
                                                                                                                    aaCCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCATCAAGGGGGTTTTTTTGtttgt                       >  6:824056/1‑149 (MQ=11)
                                                                                                                    aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCGGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  1:288770/1‑151 (MQ=255)
                                                                                                                    aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGTACACCCTAACGGGTGTTTTTTTGTTTCt                       >  3:595918/1‑151 (MQ=255)
                                                                                                                    aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  1:972538/1‑151 (MQ=255)
                                                                                                                    aaCCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCt                       >  5:242486/1‑151 (MQ=255)
                                                                                                                      ccAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATACCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg                     >  6:948128/1‑151 (MQ=255)
                                                                                                                      ccAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg                     >  1:1533932/1‑151 (MQ=255)
                                                                                                                      ccAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTgg                     >  8:1695263/1‑151 (MQ=255)
                                                                                                                       cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGCCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCAGAAGCTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAACCTAACGGGTGTTTTTTTGTTTCTGGt                    >  6:234045/1‑151 (MQ=255)
                                                                                                                       cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  1:756217/1‑151 (MQ=255)
                                                                                                                       cAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGt                    >  5:1069624/1‑151 (MQ=255)
                                                                                                                        caattaaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTGATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACGCTATATGGTGTTTTTTTGGTGCTggg                    >  1:205251/7‑150 (MQ=11)
                                                                                                                        caattaaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACAGCCTGACGGGTGTTTTTTTGTTTCTGGt                    >  7:1150751/7‑151 (MQ=14)
                                                                                                                        aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGGTTTTTTGTTTGTGGtc                   >  4:1772397/1‑151 (MQ=255)
                                                                                                                        aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  6:112556/1‑151 (MQ=255)
                                                                                                                        aaTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtc                   >  6:913777/1‑151 (MQ=255)
                                                                                                                         aattaaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAAGCTAACGGGTGATTTAGTGTTTCTGGtc                   >  8:1746609/6‑151 (MQ=14)
                                                                                                                         aTCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtct                  >  1:1513471/1‑151 (MQ=255)
                                                                                                                          attaaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAACTCCGGCCCCCGCAACCAATTATTGAACAACCAAACCGGTGTTTTTTTGGATCTGGtcg                  >  7:599789/5‑150 (MQ=14)
                                                                                                                           ttaaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTATAGCGGGTTTTATTGTTTCTGGtgtg                 >  7:1435443/4‑148 (MQ=11)
                                                                                                                           tCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGtctc                 >  1:1033296/1‑151 (MQ=255)
                                                                                                                            taaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCCTAAACGGTGATTTTTTGGTACTGGTCTcc                >  3:144528/3‑151 (MQ=14)
                                                                                                                            taaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGATTGAACACCCAAACCGGTGTTTTTTGGTTTCgggggtcc                >  4:1286252/3‑146 (MQ=14)
                                                                                                                            cAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTcc                >  4:1110226/1‑151 (MQ=255)
                                                                                                                             aaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACAGGTGATTAGTTTTTTCTGGTCTccc               >  3:470800/2‑151 (MQ=14)
                                                                                                                             aaaaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATAACACCCTAAAGCTGGTTTAATTGTTACAGGTCTcc                >  3:1779001/2‑150 (MQ=1)
                                                                                                                              aaaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  4:30091/1‑151 (MQ=255)
                                                                                                                               aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  1:569854/1‑150 (MQ=255)
                                                                                                                               aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  1:1815951/1‑150 (MQ=255)
                                                                                                                               aaTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTccc               >  6:1022334/1‑150 (MQ=255)
                                                                                                                               aaTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTTAACACCATAAAGGGTGTGTTTTTGTTTCTGGTGTCCCa              >  7:33139/1‑150 (MQ=17)
                                                                                                                                aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  3:523845/1‑151 (MQ=255)
                                                                                                                                aTTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  3:398257/1‑151 (MQ=255)
                                                                                                                                aTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             <  2:1844278/134‑1 (MQ=255)
                                                                                                                                aTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttttttt                             >  1:1846555/1‑134 (MQ=255)
                                                                                                                                aTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAAAGGGTGATTTTTTGTTTCTGGTCTCCCAt             >  8:1329638/1‑150 (MQ=255)
                                                                                                                                aTTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCATAACGGGTGTTTTTTTGtt                          >  3:1224917/1‑137 (MQ=255)
                                                                                                                                 tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  1:405456/1‑150 (MQ=255)
                                                                                                                                 tttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCAt             >  8:757677/1‑150 (MQ=255)
                                                                                                                                 tttGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGttt                                 >  6:1396394/1‑129 (MQ=255)
                                                                                                                                 tttGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTTTTTGTTTGGGGtctggaata            >  1:1799078/1‑144 (MQ=17)
                                                                                                                                 tttGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAACACCCAAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa            >  2:1153024/1‑150 (MQ=255)
                                                                                                                                  ttGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATa            >  7:1177549/1‑150 (MQ=255)
                                                                                                                                  ttGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTTAACACCAAAAAAGGTGTTTTATTTTTTCTGGTCTCCCATaa           >  3:576298/1‑150 (MQ=14)
                                                                                                                                   tgatgaAGTAAAAGCAGTACGGTGACGGGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCGCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  3:185339/1‑150 (MQ=255)
                                                                                                                                   tgatgaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaa           >  3:432459/1‑150 (MQ=255)
                                                                                                                                   tgatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAAAGGGTGTTTTGTGTTTTCTGGGCTCCCATaaa          >  1:378442/1‑150 (MQ=17)
                                                                                                                                   tgatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAATGAACACGCACAAGGCGGTTTAGTGGTTTGAGGTCTCCCATaac          >  1:387317/1‑149 (MQ=11)
                                                                                                                                    gatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGaagagc                                              >  4:969613/1‑109 (MQ=9)
                                                                                                                                    gatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACAcc                                              <  3:968670/113‑1 (MQ=25)
                                                                                                                                    gatgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAATGATCAACATAAAGGCTCTGTTTTTTTTTCGGGTGTCCCaaaaaa         >  3:123591/1‑150 (MQ=11)
                                                                                                                                     atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTTAACACCCTAACGGCTCTTTTTTTGTATCGGGTCTCCCATaaaaac       >  2:1064141/1‑150 (MQ=17)
                                                                                                                                     atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACAGGTGATTTTTTGTTTCAGGTCTCCCATaaaaa        >  4:1040719/1‑150 (MQ=255)
                                                                                                                                     atgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  3:1344999/1‑150 (MQ=255)
                                                                                                                                      tgaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGCGGGTTTTTTTGGTTCTGGTCTCCCATaaaaacg      >  3:469985/1‑149 (MQ=255)
                                                                                                                                       gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  8:1625128/1‑149 (MQ=255)
                                                                                                                                       gaAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATaaaaa        >  3:970828/1‑149 (MQ=255)
                                                                                                                                       gaAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTTTTTGTTGCGGGTGTCCCATAAAAAAgc     >  7:911309/1‑151 (MQ=18)
                                                                                                                                        aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCACCCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  4:1191812/1‑151 (MQ=255)
                                                                                                                                        aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  5:1380484/1‑151 (MQ=255)
                                                                                                                                        aaGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgc     >  8:1095409/1‑151 (MQ=255)
                                                                                                                                         aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTACCGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  7:978247/1‑151 (MQ=255)
                                                                                                                                         aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  2:416813/1‑151 (MQ=255)
                                                                                                                                         aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAgcg    >  7:470027/1‑151 (MQ=255)
                                                                                                                                         aGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTGTTTTTGTGTCTGGTCTCCAATAAAAAAgcg    >  3:1006829/1‑151 (MQ=255)
                                                                                                                                           tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  4:1346434/1‑151 (MQ=255)
                                                                                                                                           tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  3:740737/1‑151 (MQ=255)
                                                                                                                                           tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  8:97565/1‑151 (MQ=255)
                                                                                                                                           tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAACCGAAGGTCGTCGGTTCAAATCAGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  5:223534/1‑151 (MQ=255)
                                                                                                                                           tAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGCTAGCTCGTCGGGCTCATAACCCCAAGCTAGTCGCTTCAAATCCGGCACCCGCAACCAATTATTGAACACCCTAACGCGTGTTTTTTTGTTTCTAGTCTCCCATAAAAAAGCGcc  >  2:1150596/1‑151 (MQ=21)
                                                                                                                                            aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  1:921471/1‑150 (MQ=255)
                                                                                                                                            aaaaGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGGGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGcc  >  2:52338/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGTGTTTTTTTGTTTCTGGTCTCCCATAAAAAAGCGCC  >  NC_000913/2947453‑2947740

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