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→CAAdnaKchaperone 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

ACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGC  >  NC_000913/12340‑12632
                                                                                                                                                      |                                                                                                                                                  
aCGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGaaa                                                                                                                                                 >  2:872882/1‑150 (MQ=255)
 cGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAt                                                                                                                                               <  2:55305/151‑1 (MQ=255)
   aaCCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATgg                                                                                                                                             <  2:291406/151‑1 (MQ=255)
         cAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAccgc                                                                                                                                       >  2:1211522/1‑151 (MQ=255)
          aaaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAccgcc                                                                                                                                      >  2:1095747/1‑151 (MQ=255)
          aaaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGATTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAccgcc                                                                                                                                      >  2:894738/1‑151 (MQ=255)
           aaaCACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCAccgccg                                                                                                                                     >  2:187289/1‑151 (MQ=255)
                      ttGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTctgc                                                                                                                          >  1:1172611/1‑151 (MQ=255)
                        gCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTctgctg                                                                                                                        >  1:795140/1‑151 (MQ=255)
                           aTTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGa                                                                                                                       >  1:626628/1‑149 (MQ=255)
                             tAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAg                                                                                                                     >  1:534049/1‑149 (MQ=255)
                               accaccttttggttCGCCGCTCCCAGGACAAAGAAGTTCAGGGTGATGTTACCAACATGCTGTAAAAAATTATTGTTGCGGATAACGGCGACGCATGGCTCGAATTAAAGGGCCAGAAAGTGCCACCGCCGCAGATTTCTGCTGAAGTGCt                                                                                                                 <  1:445029/138‑1 (MQ=255)
                               aaCGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCt                                                                                                                 >  1:275306/1‑151 (MQ=255)
                                   ccTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTAAAAATTATGGCTGCTGAAAACGGCGACGCATGGGTCGAATTAAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaa                                                                                                              <  2:859743/150‑1 (MQ=255)
                                   ccTGATTGGTCGCCGCTTCCAGGACGAAAAAGTACAGCGTGATGTTTCCACCATGCCGTTAAAAATTTTTGCTGCTGATACCGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTg                                                                                                                <  1:90671/148‑1 (MQ=255)
                                       aTTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          >  1:90166/1‑150 (MQ=255)
                                       aTTGGTCGCCGCTTCCAGGACAAAGAAGTACAGGGTGATGTTCCCATCATGCCTTTAAAAATTATTGCTGCTAATACCGGCGCCGCTTGGGTCGAATTTAAAGGCCAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          <  1:809808/150‑1 (MQ=255)
                                        ttGTTCGCCGTTCCCAGGCCAAATAAGCACGGGGGAATGTTCCCCCCATCCGTTTAAAAATTTTTGCGGCTAATACCGGCAACGCATGGGGAGAATTAAAAGGAAAAAAAATGCCACCGCCGCAGATTTCTGCTGAAGTGCTGaaaaaa                                                                                                          <  2:1233586/149‑1 (MQ=255)
                                          ggTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATg                                                                                                       >  2:1219898/1‑150 (MQ=255)
                                            tCGCCGCTCCAAGGACAAAGAATTACAGCGTGTTGTTTCCATCATGCCGTTAAAAATTTTTGCTGCTGATAACGGCGACGCATGGGTGAAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaa                                                                                                     <  2:369716/150‑1 (MQ=255)
                                             cgccccTTCCAGGACGAAGAAGTAGAGCGTGATTTTTCCATCATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGTGGAATTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaa                                                                                                     <  2:1153006/149‑1 (MQ=255)
                                               ccgcTCCCAGGACGAAGAAGTAGAGCGGGATGTTCACCTCCTGCCGTTCAAAATTATGGTTGCGGAAACCGGCGACGAATGGGTCGAAGAAAAAGGCAAAAAAATGGCTCCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGaaa                                                                                                 <  2:1172547/151‑1 (MQ=255)
                                                cgcTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGaaa                                                                                                 >  2:598266/1‑150 (MQ=255)
                                                cgcTTCCAGAACAAAGAAGTACAGCGTGATTTTTCCATCATGCCGTAAAAAATTATTGCTGCTGAAAACGACGACGAATGGGTCGAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATgaaga                                                                                                   <  2:1035541/148‑1 (MQ=255)
                                                   ttCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCg                                                                                             >  1:1179681/1‑151 (MQ=255)
                                                     caaGGACGAGGAAGTACAGCTTGATTTTTCCTTCTTGCCGTTCAAATTTTTTGCTGCTGATAACGCCGCCGCATGGGTCGAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCg                                                                                             <  1:501454/147‑1 (MQ=255)
                                                             aagaagTACAGCGTGTTGTTTCCATCATGCCTTTAAAAATTATTGCTGCTGAAAACGGCGACGCATGGGTGGAAGTAAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                    <  2:758665/150‑1 (MQ=255)
                                                             aagaagTACAGCGTGATGTTGCCATCATGCCGTTCAAAATTATTGCTGCTGAAACCGGCGACGCATGGGTCGAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGAtt                                                                                    <  2:605776/150‑1 (MQ=255)
                                                                 agTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTg                                                                               >  1:700767/1‑151 (MQ=255)
                                                                     cAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTg                                                                           >  1:1042636/1‑151 (MQ=255)
                                                                       gCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGAATGGGTCGAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGa                                                                          <  2:789234/150‑1 (MQ=255)
                                                                          tgttttttCCATCATCCCGTTCAAAATGATGGCTGCTAAAAGCGGCGCGCCTTGGTCCGAATTTAAAGCCCAGAAAAGGGAACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAc                                                                        <  1:1245346/144‑1 (MQ=255)
                                                                            aTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  1:175095/1‑151 (MQ=255)
                                                                            aTGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  2:546183/1‑151 (MQ=255)
                                                                             tGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGt                                                                    >  2:434880/1‑150 (MQ=255)
                                                                               tttCCATCATGCCTTTAAAAATTATTGCTGCTAATAACGGGGAGGCAGGGGTGGAATTTAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCgg                                                                     <  2:413219/147‑1 (MQ=255)
                                                                               tttCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAAc                                                                 >  2:1207873/1‑151 (MQ=255)
                                                                               tttCCATCATGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAAcc                                                                       <  2:1034025/145‑1 (MQ=255)
                                                                               tttCCAACAATGCCGTAAAAAATTATTGCGGCAGAAACCCGCGGCGAATGGGTGAAATTAAAGGCCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGaa                                                                         <  2:1268708/144‑1 (MQ=255)
                                                                                      catGCCGTTAAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTAAAAGGCCAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGc                                                           <  1:780733/150‑1 (MQ=255)
                                                                                         tGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        >  1:1165796/1‑151 (MQ=255)
                                                                                         tGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        >  2:48179/1‑151 (MQ=255)
                                                                                              ttCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGCTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATc                                                    >  1:13979/1‑150 (MQ=255)
                                                                                              ttCAAAATTATGGCTGCTGATAACGGCGACGAATGGGTCGAATTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGt                                                        <  1:60833/146‑1 (MQ=255)
                                                                                               taaaaaTTATTGCTGTAGAAACCGGCGACGGATGGGTCGAATTAAAAGGACAAAAAAGGGAACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTAt                                                     <  2:770971/146‑1 (MQ=255)
                                                                                                   aaTTATTGCTGCTGATACCGGCGACGAATGGGTCGAAGTAAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTa                                              <  2:766858/151‑1 (MQ=255)
                                                                                                   aaTTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAGGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGt                                               >  1:1276607/1‑150 (MQ=255)
                                                                                                     ttattGCTGCTAATAACGGCGACGCATGGGTCGAATTAAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTAc                                             <  2:530069/150‑1 (MQ=255)
                                                                                                       attGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          >  2:4437/1‑151 (MQ=255)
                                                                                                        ttGCTGCTAATACCGGCGACGAATGGGTGGAAGTTAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          <  2:881192/150‑1 (MQ=255)
                                                                                                         tgctgcTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCgg                                          >  2:351469/1‑149 (MQ=255)
                                                                                                         tgctgcTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGc                                         >  1:1059856/1‑150 (MQ=255)
                                                                                                          gctgctGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCAt                                       >  1:1031311/1‑151 (MQ=255)
                                                                                                              ctgATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACtt                                   >  1:1138036/1‑151 (MQ=255)
                                                                                                                gATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  >  1:1016034/1‑150 (MQ=255)
                                                                                                                gATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  >  1:458135/1‑150 (MQ=255)
                                                                                                                aaTAGCGGCGACTGATGGGTGGAGTTTAAAGGCCAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACt                                    <  1:240874/147‑1 (MQ=255)
                                                                                                                 aaagcGGCGAGGCATGGGTAGAAGTAAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTaa                                <  2:543730/147‑1 (MQ=255)
                                                                                                                 aTAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTa                                 >  2:218108/1‑150 (MQ=255)
                                                                                                                 aTAACGGCGACGCAGGGGTCGAATTTAAAGCCCAGAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACttt                                  <  1:21234/149‑1 (MQ=255)
                                                                                                                  tAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAAc                               >  1:670/1‑151 (MQ=255)
                                                                                                                     cGGCGACGCATGGGTCGAAGTTAAAGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACg                              <  2:99282/149‑1 (MQ=255)
                                                                                                                      ggCGACGCATGGGTCAAAGTTAAGGGCCAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATg                           <  1:1095807/151‑1 (MQ=255)
                                                                                                                      ggCGACGCATGGGTGAAGTTTAAAGGCCAGAAAAGGACACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGAt                            <  1:759590/150‑1 (MQ=255)
                                                                                                                         gACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTc                        >  2:1147570/1‑151 (MQ=255)
                                                                                                                            gaaTGGGTCGAAGTTAAAGGCAAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCa                       <  2:93435/147‑1 (MQ=255)
                                                                                                                                ggggacaaaGTAAAAGGCCAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTc                  <  2:1021601/144‑1 (MQ=255)
                                                                                                                                gggTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTc                  >  1:1227427/1‑150 (MQ=255)
                                                                                                                                      gAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCa             >  1:1333633/1‑150 (MQ=255)
                                                                                                                                             aaaGGCCAAAAAAGGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCa         <  2:795113/148‑1 (MQ=255)
                                                                                                                                             aaaGGCAAAAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCaaa       <  1:881946/150‑1 (MQ=255)
                                                                                                                                                   cAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGc  >  2:775895/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACGAACCCGCAAAACACTCTGTTTGCGATTAAACGCCTGATTGGTCGCCGCTTCCAGGACGAAGAAGTACAGCGTGATGTTTCCATCATGCCGTTCAAAATTATTGCTGCTGATAACGGCGACGCATGGGTCGAAGTTAAAGGCCAGAAAATGGCACCGCCGCAGATTTCTGCTGAAGTGCTGAAAAAAATGAAGAAAACCGCTGAAGATTACCTGGGTGAACCGGTAACTGAAGCTGTTATCACCGTACCGGCATACTTTAACGATGCTCAGCGTCAGGCAACCAAAGACGC  >  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: