Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091312,4860GA38.9% 112.5 / 28.6 72Q108Q (CAG→CAA) dnaKchaperone Hsp70, with co‑chaperone DnaJ
Reads supporting (aligned to +/- strand):  ref base G (37/7);  new base A (0/28);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 8.90e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

ACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGC  >  NC_000913/12340‑12632
                                                                                                                                                      |                                                                                                                                                  
aCGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGaaa                                                                                                                                                 >  2:872882/1‑150 (MQ=255)
 cGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAAt                                                                                                                                               <  2:55305/151‑1 (MQ=255)
   aaCCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATgg                                                                                                                                             <  2:291406/151‑1 (MQ=255)
         cAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCAccgc                                                                                                                                       >  2:1211522/1‑151 (MQ=255)
          aaaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCAccgcc                                                                                                                                      >  2:1095747/1‑151 (MQ=255)
          aaaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGATTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCAccgcc                                                                                                                                      >  2:894738/1‑151 (MQ=255)
           aaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCAccgccg                                                                                                                                     >  2:187289/1‑151 (MQ=255)
                      ttGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTctgc                                                                                                                          >  1:1172611/1‑151 (MQ=255)
                        gCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTctgctg                                                                                                                        >  1:795140/1‑151 (MQ=255)
                           aTTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGa                                                                                                                       >  1:626628/1‑149 (MQ=255)
                             tAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAg                                                                                                                     >  1:534049/1‑149 (MQ=255)
                               accaccttttggttCGCCGCTCCCAGGACAAAGAAGTTCAGGGTGATGTTACCAAC‑ATGCTGTAAAAAATTATTGTTGCGGATAACGGCGACGCATGGCT‑CGAAT‑TAAAGGGCC‑AGAAAGTGCCACCGCCGCAGATTTCTGCTGAAGTGCt                                                                                                                 <  1:445029/138‑1 (MQ=255)
                               aaCGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCt                                                                                                                 >  1:275306/1‑151 (MQ=255)
                                   ccTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTAAAAATTATGGCTGCTGAAAACGGCGACGCATGGGT‑CGAAT‑TAAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaa                                                                                                              <  2:859743/150‑1 (MQ=255)
                                   ccTGATTGGTCGCCGCTTCCAGGACGAAAAAGTACAGCGTGATGTTTCCACC‑ATGCCGTTAAAAATTTTTGCTGCTGATACCGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTg                                                                                                                <  1:90671/148‑1 (MQ=255)
                                       aTTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          >  1:90166/1‑150 (MQ=255)
                                       aTTGGTCGCCGCTTCCAGGACAAAGAAGTACAGGGTGATGTTCCCATC‑ATGCCTTTAAAAATTATTGCTGCTAATACCGGCGCCGCTTGGGT‑CGAAT‑TTAAAGGCC‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          <  1:809808/150‑1 (MQ=255)
                                        ttGTTCGCCGTTCCCAGGCCAAATAAGCACGGGGGAATGTTCCCCCC‑ATCCGTTTAAAAATTTTTGCGGCTAATACCGGCAACGCATGGGG‑AGAAT‑TAAAAGGAA‑AAAAAATGCCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          <  2:1233586/149‑1 (MQ=255)
                                          ggTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATg                                                                                                       >  2:1219898/1‑150 (MQ=255)
                                            tCGCCGCTCCAAGGACAAAGAATTACAGCGTGTTGTTTCCATC‑ATGCCGTTAAAAATTTTTGCTGCTGATAACGGCGACGCATGGGT‑GAAAG‑TAAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaa                                                                                                     <  2:369716/150‑1 (MQ=255)
                                             cgccccTTCCAGGACGAAGAAGTAGAGCGTGATTTTTCCATC‑ATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑GGAAT‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaa                                                                                                     <  2:1153006/149‑1 (MQ=255)
                                               ccgcTCCCAGGACGAAGAAGTAGAGCGGGATGTTCACCTC‑CTGCCGTTCAAAATTATGGTTGCGGAAACCGGCGACGAATGGGT‑CGAAG‑AAAAAGGCA‑AAAAAATGGCTCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGaaa                                                                                                 <  2:1172547/151‑1 (MQ=255)
                                                cgcTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGaaa                                                                                                 >  2:598266/1‑150 (MQ=255)
                                                cgcTTCCAGAACAAAGAAGTACAGCGTGATTTTTCCATC‑ATGCCGTAAAAAATTATTGCTGCTGAAAACGACGACGAATGGGT‑CGAAG‑TAAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaaga                                                                                                   <  2:1035541/148‑1 (MQ=255)
                                                   ttCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCg                                                                                             >  1:1179681/1‑151 (MQ=255)
                                                     caaGGACGAGGAAGTACAGCTTGATTTTTCCTTC‑TTGCCGTTCAAATTTTTTGCTGCTGATAACGCCGCCGCATGGGT‑CGAAG‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCg                                                                                             <  1:501454/147‑1 (MQ=255)
                                                             aagaagTACAGCGTGTTGTTTCCATC‑ATGCCTTTAAAAATTATTGCTGCTGAAAACGGCGACGCATGGGT‑GGAAG‑TAAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                    <  2:758665/150‑1 (MQ=255)
                                                             aagaagTACAGCGTGATGTTGCCATC‑ATGCCGTTCAAAATTATTGCTGCTGAAACCGGCGACGCATGGGT‑CGAAG‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                    <  2:605776/150‑1 (MQ=255)
                                                                 agTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTg                                                                               >  1:700767/1‑151 (MQ=255)
                                                                     cAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTg                                                                           >  1:1042636/1‑151 (MQ=255)
                                                                       gCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGAATGGGT‑CGAAG‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGa                                                                          <  2:789234/150‑1 (MQ=255)
                                                                          tgttttttCCATC‑ATCCCGTTCAAAATGATGGCTGCTAAAAGCGGCG‑CGCCTTGGTCCGAAT‑TTAAAGCCC‑AGAAAAGGGAACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                        <  1:1245346/144‑1 (MQ=255)
                                                                            aTGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  1:175095/1‑151 (MQ=255)
                                                                            aTGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  2:546183/1‑151 (MQ=255)
                                                                             tGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  2:434880/1‑150 (MQ=255)
                                                                               tttCCATC‑ATGCCTTTAAAAATTATTGCTGCTAATAACGGGGAGGCAGGGGT‑GGAAT‑TTAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCgg                                                                     <  2:413219/147‑1 (MQ=255)
                                                                               tttCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAAc                                                                 >  2:1207873/1‑151 (MQ=255)
                                                                               tttCCATC‑ATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAcc                                                                       <  2:1034025/145‑1 (MQ=255)
                                                                               tttCCAACAATGCCGTAAAAAATTATTGCGGCAGAAACCCGCGGCGAATGGGT‑‑GAAA‑TTAAAGGCCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGaa                                                                         <  2:1268708/144‑1 (MQ=255)
                                                                                      c‑atGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TAAAAGGCC‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGc                                                           <  1:780733/150‑1 (MQ=255)
                                                                                         tGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        >  1:1165796/1‑151 (MQ=255)
                                                                                         tGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        >  2:48179/1‑151 (MQ=255)
                                                                                              ttCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑CTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATc                                                    >  1:13979/1‑150 (MQ=255)
                                                                                              ttCAAAATTATGGCTGCTGATAACGGCGACGAATGGGT‑CGAAT‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        <  1:60833/146‑1 (MQ=255)
                                                                                               taaaaaTTATTGCTGTAGAAACCGGCGACGGATGGGT‑CGAAT‑TAAAAGGAC‑AAAAAAGGGAACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTAt                                                     <  2:770971/146‑1 (MQ=255)
                                                                                                   aaTTATTGCTGCTGATACCGGCGACGAATGGGT‑CGAAG‑TAAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTa                                              <  2:766858/151‑1 (MQ=255)
                                                                                                   aaTTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAGGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGt                                               >  1:1276607/1‑150 (MQ=255)
                                                                                                     ttattGCTGCTAATAACGGCGACGCATGGGT‑CGAAT‑TAAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTAc                                             <  2:530069/150‑1 (MQ=255)
                                                                                                       attGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          >  2:4437/1‑151 (MQ=255)
                                                                                                        ttGCTGCTAATACCGGCGACGAATGGGT‑GGAAG‑TTAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          <  2:881192/150‑1 (MQ=255)
                                                                                                         tgctgcTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          >  2:351469/1‑149 (MQ=255)
                                                                                                         tgctgcTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGc                                         >  1:1059856/1‑150 (MQ=255)
                                                                                                          gctgctGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAt                                       >  1:1031311/1‑151 (MQ=255)
                                                                                                              ctgATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACtt                                   >  1:1138036/1‑151 (MQ=255)
                                                                                                                gATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  >  1:1016034/1‑150 (MQ=255)
                                                                                                                gATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  >  1:458135/1‑150 (MQ=255)
                                                                                                                aaTAGCGGCGACTGATGGGT‑GGAGT‑TTAAAGGCC‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACt                                    <  1:240874/147‑1 (MQ=255)
                                                                                                                 aaagcGGCGAGGCATGGGT‑AGAAG‑TAAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTaa                                <  2:543730/147‑1 (MQ=255)
                                                                                                                 aTAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTa                                 >  2:218108/1‑150 (MQ=255)
                                                                                                                 aTAACGGCGACGCAGGGGT‑CGAAT‑TTAAAGCCC‑AGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  <  1:21234/149‑1 (MQ=255)
                                                                                                                  tAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAc                               >  1:670/1‑151 (MQ=255)
                                                                                                                     cGGCGACGCATGGGT‑CGAAG‑TTAAAGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACg                              <  2:99282/149‑1 (MQ=255)
                                                                                                                      ggCGACGCATGGGT‑CAAAG‑TTAAGGGCC‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATg                           <  1:1095807/151‑1 (MQ=255)
                                                                                                                      ggCGACGCATGGGT‑‑GAAGTTTAAAGGCC‑AGAAAAGGACACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGAt                            <  1:759590/150‑1 (MQ=255)
                                                                                                                         gACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTc                        >  2:1147570/1‑151 (MQ=255)
                                                                                                                            gaaTGGGT‑CGAAG‑TTAAAGGCA‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCa                       <  2:93435/147‑1 (MQ=255)
                                                                                                                                ggggacaaaG‑TAAAAGGCC‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTc                  <  2:1021601/144‑1 (MQ=255)
                                                                                                                                gggT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTc                  >  1:1227427/1‑150 (MQ=255)
                                                                                                                                      gAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCa             >  1:1333633/1‑150 (MQ=255)
                                                                                                                                             aaaGGCC‑AAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCa         <  2:795113/148‑1 (MQ=255)
                                                                                                                                             aaaGGCA‑AAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaaa       <  1:881946/150‑1 (MQ=255)
                                                                                                                                                   c‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGc  >  2:775895/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATC‑ATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGT‑CGAAG‑TTAAAGGCC‑AGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGC  >  NC_000913/12340‑12632

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑