Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913606,6120GA46.1% 94.7 / 34.8 76Q258* (CAG→TAG) ybdKweak gamma‑glutamyl:cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (31/10);  new base A (0/35);  total (31/45)
Fisher's exact test for biased strand distribution p-value = 5.92e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  NC_000913/606470‑606759
                                                                                                                                                    |                                                                                                                                                   
cgacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATc                                                                                                                                             <  1:1260491/151‑1 (MQ=255)
 gacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAAt                                                                                                                                              <  2:1126204/149‑1 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCCGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGAGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  2:226928/1‑151 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  1:55766/1‑151 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  2:716417/1‑151 (MQ=255)
     aTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGcc                                                                                                                                        >  2:896802/1‑151 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCg                                                                                                                                          <  1:895498/145‑1 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCg                                                                                                                                          <  1:943409/145‑1 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCa                                                                                                                                       <  1:5337/148‑1 (MQ=255)
                ggggaggggTGAGGCCCCCTTCAACCCAATGGCGACGGGCGTGAAAA‑CGGTTGAA‑TTTATCAAGCAGGTATTTTTTTTCC‑TG‑ATTTTAAAC‑C‑GGGCTTTCCGTCATTACCCGGGGGGGGGTAGCATAAATTAATCCCGCCATATTTACTGc                                                                                                                             <  2:127240/144‑1 (MQ=255)
                     tCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGc                                                                                                                        >  2:1193751/1‑151 (MQ=255)
                      gggTGATGCCCCCTTCACCCCAAAAGCGACGGGCGTAAAAG‑GGGTTAAA‑TTTATACAGCAGGTATTTTTTTTCG‑T‑‑TTTTTTAAA‑CGGGGCTTTCCGTAAGTACCGGGGGGGGGGTAGCTTAAATAAATCCCGCCATATTTACTGCGTGGCt                                                                                                                       <  2:1027552/150‑1 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:1139802/1‑151 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:1219349/1‑151 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:937224/1‑151 (MQ=255)
                               ccccTACAACCCCATAGCGAAGGGCCTGAAAA‑GGGTTAAT‑TTTATACAGCAGGTAATTTTTTCCC‑TG‑ATTTTTAAG‑C‑GGGCGTTCCGTCAGTACCCGGGGGGGGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAgg                                                                                                                   <  1:1094054/144‑1 (MQ=255)
                                 ccTTCAACCCCATAGCGAC‑GGCCGGGGAA‑GGTTTTAATTTTATACAGCAGGTATTTTTTTTCC‑TT‑TTTTTAAAG‑G‑GGGCTTTCCTTCAGTACCCGGGGGGGGGTACCCTAAATAAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAAc                                                                                                             <  2:928257/150‑1 (MQ=255)
                                   ttCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACggg                                                                                                          >  1:1007474/1‑151 (MQ=255)
                                   ttCAACCCCAAAGCGACAGG‑CTGGAAAGGGGTTAAT‑TTTATACAGAAGGTATTTTTTTTCC‑TG‑ATTTTAAAG‑C‑GGGCTTTCCGTCAGAACCCGGGGGGGGGTGGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAAc                                                                                                             <  1:720375/148‑1 (MQ=255)
                                      aaccccATAGCGACAGGCCGGAAAA‑CGTTTAAA‑TTTATACAGGAGGTATTTTTTTTCC‑TG‑ATTTTAAAA‑C‑GGGCGTTCCGTCAGTACCCGGGGGGCGGTAGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTa                                                                                                               <  1:1077876/140‑1 (MQ=255)
                                           cATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCAt                                                                                                  >  1:327178/1‑151 (MQ=255)
                                            aTAGCGACAGGCTTGAAAA‑CGTTTAAT‑TTTATACAGCAGTTAATTTTTTCCC‑TG‑TTGTTTAAG‑C‑GGGCGTTCCGTCAGTACCGGGGGGGGGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATc                                                                                                 <  1:812651/151‑1 (MQ=255)
                                            aTAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATc                                                                                                 >  1:944115/1‑151 (MQ=255)
                                                  aCAGGCTGGGAAG‑CGTTTGAA‑TTTATACGGCAGTAATTTTTTTTCG‑TT‑TTTTTTAAG‑C‑GGGCGTCCCGTAATTACCCGGGGGGCGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAc                                                                                               <  2:932834/147‑1 (MQ=255)
                                                   cgggCCTGGAAA‑CGTTTGAA‑TTTATACGGCAGGTAAT‑TTTTTCC‑TT‑GTTTTTAAAGG‑GGGCGTTCCGAAATTAACCGGTGGGGGGTGGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGa                                                                                           <  1:690434/148‑1 (MQ=255)
                                                    gggCTTGGAAG‑CGTTTGAA‑TTTAT‑CAAGATGTTATTTTTTTCCGTG‑GTTTTAAAC‑C‑GGGCGTCCCTCCAGAACCTAGGGGGGGGTAGCTAAAATTATTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAc                                                                                         <  1:1003103/150‑1 (MQ=255)
                                                    aGGCATGGAAA‑CGTTTAAA‑TTTATACAGCAGTTATTTTTTTTCC‑TG‑ATTTTAAAG‑C‑GGGCTTTCCTTCAGTAACGGGTGGGGGGTAGCGTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCg                                                                                            <  1:730842/148‑1 (MQ=255)
                                                     ggCGTGAAAA‑GGTTTGAT‑TTTATACAGCAGTTATTTTTTTTCC‑TG‑ATTTTTAAC‑C‑GGGCTTCCCGTCAGAACCGGGGGGGGGGTGGCCGAAATAACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAcc                                                                                        <  2:960691/151‑1 (MQ=255)
                                                     ggCCTGAAAA‑GGTTTAAT‑TTTATACAGGAGTTAATTTTTTTCC‑TG‑ATTTTAAAG‑G‑GGGCTTTCCTTAAGTACCCGGGGGGGGTTAGCCTAATTTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGaa                                                                                          <  1:1185455/149‑1 (MQ=255)
                                                      gCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCt                                                                                       >  2:796348/1‑151 (MQ=255)
                                                            aaa‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCg                                                                                 >  2:109135/1‑151 (MQ=255)
                                                               agggttttattttATACAG‑AGGTTATTTTTTTCC‑AGTTTTTTTAAG‑C‑GGGCCTTCCTAAAGAACCCGGTGGGGGGCTCCAGAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGt                                                                                <  1:38146/141‑1 (MQ=255)
                                                                cGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGc                                                                              >  2:184984/1‑151 (MQ=255)
                                                                    tGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaa                                                                           >  1:468499/1‑150 (MQ=255)
                                                                      aa‑TTTATACAGCAGTTATTTTTTTTCC‑TG‑TTGTTAAAA‑G‑GGGCGTCCCGTCAGAAACCAGTGGGCGGTACCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAAt                                                                        <  2:632578/151‑1 (MQ=255)
                                                                          ttATACAGCAGGTATTTTTTTTCT‑TG‑ATTTTAAAG‑C‑GGGCGTCCCGTCAGAACCCAGGGGGGGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGa                                                                      <  1:202555/150‑1 (MQ=255)
                                                                              aCAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACt                                                                 >  1:29146/1‑151 (MQ=255)
                                                                                   agGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgggg                                                             >  1:487993/1‑150 (MQ=255)
                                                                                    gttattttttttCT‑TG‑ATGTTAAAG‑C‑GGGCGTCCCGTCAGAAACCGGGGGGGGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  2:506098/144‑1 (MQ=255)
                                                                                    ggTAATTTTTTTCC‑TG‑ATGTTTAAG‑G‑GGGCTTTCGGTAAGTACCAGGTGGGCGGTCGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  1:1095635/151‑1 (MQ=255)
                                                                                      tAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGa                                                          >  1:1235974/1‑150 (MQ=255)
                                                                                           tttttCC‑TG‑TTGTTAAAC‑C‑GGGCGTCCCGTAAGTACCCGGGGGGGGGTGGCTTAAATTACTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  1:103735/144‑1 (MQ=255)
                                                                                           tttttCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATcc                                                    >  1:740720/1‑151 (MQ=255)
                                                                                              ttCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc                                                   >  2:684422/1‑149 (MQ=255)
                                                                                                cc‑TG‑TTGTTAAAG‑G‑GGGCGTCCCGTCAGTACCCGGGGGGCGTTACCTTAAATAAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGt                                                <  1:964002/150‑1 (MQ=255)
                                                                                                cc‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCCGTg                                               >  1:1145283/1‑151 (MQ=255)
                                                                                                   tG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  2:610177/1‑150 (MQ=255)
                                                                                                       tttttAAG‑C‑GGGCGTTCCGTAAGTAACGGGTGGGGGGTACCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              <  1:1259723/145‑1 (MQ=255)
                                                                                                        ttttAAACC‑GGGCTTCCCGTCAGTAACCGGGGGGCGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAg                                                 <  1:65046/143‑1 (MQ=255)
                                                                                                        ttttAAACC‑GGGCGTTCCGTAAGTAACCGGGGGGGGGGAGCCAGAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAg                                            <  2:360338/148‑1 (MQ=255)
                                                                                                        ttttAAA‑CGGGGCGTCCCGTAAGAACCCGGGGGGGGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTg                                               <  2:468497/145‑1 (MQ=255)
                                                                                                           tAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCttt                                      >  1:1059203/1‑151 (MQ=255)
                                                                                                             aagc‑GGGCGTCCCGTCAGAACCGGGGGGGCGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCt                                        <  1:159217/145‑1 (MQ=255)
                                                                                                             aa‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTa                                     >  2:40559/1‑150 (MQ=255)
                                                                                                              aacggggCTTTCCGTAAGTAACCGGGGGGGGGTAGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTa                                     <  1:1182017/147‑1 (MQ=255)
                                                                                                                  gggCTTCCCGTCAATAACCGGGGGGCGGTCGCTTAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTaa                                    <  2:1257330/148‑1 (MQ=255)
                                                                                                                   ggCTTTCCGTAAGAAACCGGGGGGCGGTAGCAAAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGc                                 <  1:850236/150‑1 (MQ=255)
                                                                                                                    gCGTTCCGTCAGTAACCAGTGGGCGGAAGCCTGAATTAATCCCGCCTTATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCt                                >  2:128102/1‑150 (MQ=255)
                                                                                                                       ttCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCg                            >  1:1130751/1‑151 (MQ=255)
                                                                                                                       ttCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCg                            >  1:208122/1‑151 (MQ=255)
                                                                                                                       ttCCGTCAAAACCAAGGGGGGGGTGGCGTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGt                              <  1:1139817/149‑1 (MQ=255)
                                                                                                                         ccGTAAGAAACCGGGGGGCGGTAGCGTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTg                               <  1:131769/146‑1 (MQ=255)
                                                                                                                            taaggaacaggggggCGGTGACTTGATTTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGatca                        <  1:555530/138‑1 (MQ=255)
                                                                                                                              aGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCg                     >  1:208829/1‑151 (MQ=255)
                                                                                                                               gTAACCAGGGGGCGGTGCCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCg                     <  1:443304/150‑1 (MQ=255)
                                                                                                                                 accttgTGGGGGGTACCCAAATTCAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgt                    <  2:1006992/144‑1 (MQ=255)
                                                                                                                                 aaCCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtgg                  >  1:46584/1‑151 (MQ=255)
                                                                                                                                      gTGGGGGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTa              <  1:857380/150‑1 (MQ=255)
                                                                                                                                        gggggttAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCgtggtg                <  2:45518/140‑1 (MQ=255)
                                                                                                                                        ggggggTCGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACt           <  1:1086989/151‑1 (MQ=255)
                                                                                                                                          ggggTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCa         <  1:716421/149‑1 (MQ=255)
                                                                                                                                             gaaGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGAc      <  1:86093/149‑1 (MQ=255)
                                                                                                                                             gTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCACCCGTCCCAAAATGAGGACTGGGGCGAATACCCCAGTCCGGATCTTTAATCCTGTCGATAATGGTGGGGTaa             >  1:703177/1‑144 (MQ=255)
                                                                                                                                             gTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGAc      <  1:786368/151‑1 (MQ=255)
                                                                                                                                                 ccTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCg  <  1:266305/151‑1 (MQ=255)
                                                                                                                                                  cTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCg  >  1:257944/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAA‑CGGTTGAA‑TTTATACAGCAGGTAATCTTTTTCC‑TG‑ATGTTTAAA‑C‑GGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  NC_000913/606470‑606759

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