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

CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  NC_000913/606470‑606759
                                                                                                                                                    |                                                                                                                                                   
cgacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATc                                                                                                                                             <  1:1260491/151‑1 (MQ=255)
 gacgaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAAt                                                                                                                                              <  2:1126204/149‑1 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCCGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGAGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  2:226928/1‑151 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  1:55766/1‑151 (MQ=255)
    gaTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGc                                                                                                                                         >  2:716417/1‑151 (MQ=255)
     aTCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGcc                                                                                                                                        >  2:896802/1‑151 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCg                                                                                                                                          <  1:895498/145‑1 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCg                                                                                                                                          <  1:943409/145‑1 (MQ=255)
         ccAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCa                                                                                                                                       <  1:5337/148‑1 (MQ=255)
                ggggaggggTGAGGCCCCCTTCAACCCAATGGCGACGGGCGTGAAAACGGTTGAATTTATCAAGCAGGTATTTTTTTTCCTGATTTTAAACCGGGCTTTCCGTCATTACCCGGGGGGGGGTAGCATAAATTAATCCCGCCATATTTACTGc                                                                                                                             <  2:127240/144‑1 (MQ=255)
                     tCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGc                                                                                                                        >  2:1193751/1‑151 (MQ=255)
                      gggTGATGCCCCCTTCACCCCAAAAGCGACGGGCGTAAAAGGGGTTAAATTTATACAGCAGGTATTTTTTTTCGTTTTTTTAAACGGGGCTTTCCGTAAGTACCGGGGGGGGGGTAGCTTAAATAAATCCCGCCATATTTACTGCGTGGCt                                                                                                                       <  2:1027552/150‑1 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:1139802/1‑151 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:1219349/1‑151 (MQ=255)
                           atgaCGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAggg                                                                                                                  >  2:937224/1‑151 (MQ=255)
                               ccccTACAACCCCATAGCGAAGGGCCTGAAAAGGGTTAATTTTATACAGCAGGTAATTTTTTCCCTGATTTTTAAGCGGGCGTTCCGTCAGTACCCGGGGGGGGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAgg                                                                                                                   <  1:1094054/144‑1 (MQ=255)
                                 ccTTCAACCCCATAGCGACGGCCGGGGAAGGTTTTAATTTTATACAGCAGGTATTTTTTTTCCTTTTTTTAAAGGGGGCTTTCCTTCAGTACCCGGGGGGGGGTACCCTAAATAAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAAc                                                                                                             <  2:928257/150‑1 (MQ=255)
                                   ttCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACggg                                                                                                          >  1:1007474/1‑151 (MQ=255)
                                   ttCAACCCCAAAGCGACAGGCTGGAAAGGGGTTAATTTTATACAGAAGGTATTTTTTTTCCTGATTTTAAAGCGGGCTTTCCGTCAGAACCCGGGGGGGGGTGGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAAc                                                                                                             <  1:720375/148‑1 (MQ=255)
                                      aaccccATAGCGACAGGCCGGAAAACGTTTAAATTTATACAGGAGGTATTTTTTTTCCTGATTTTAAAACGGGCGTTCCGTCAGTACCCGGGGGGCGGTAGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTa                                                                                                               <  1:1077876/140‑1 (MQ=255)
                                           cATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCAt                                                                                                  >  1:327178/1‑151 (MQ=255)
                                            aTAGCGACAGGCTTGAAAACGTTTAATTTTATACAGCAGTTAATTTTTTCCCTGTTGTTTAAGCGGGCGTTCCGTCAGTACCGGGGGGGGGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATc                                                                                                 <  1:812651/151‑1 (MQ=255)
                                            aTAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATc                                                                                                 >  1:944115/1‑151 (MQ=255)
                                                  aCAGGCTGGGAAGCGTTTGAATTTATACGGCAGTAATTTTTTTTCGTTTTTTTTAAGCGGGCGTCCCGTAATTACCCGGGGGGCGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCAc                                                                                               <  2:932834/147‑1 (MQ=255)
                                                   cgggCCTGGAAACGTTTGAATTTATACGGCAGGTAATTTTTTCCTTGTTTTTAAAGGGGGCGTTCCGAAATTAACCGGTGGGGGGTGGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGa                                                                                           <  1:690434/148‑1 (MQ=255)
                                                    gggCTTGGAAGCGTTTGAATTTATCAAGATGTTATTTTTTTCCGTGGTTTTAAACCGGGCGTCCCTCCAGAACCTAGGGGGGGGTAGCTAAAATTATTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAc                                                                                         <  1:1003103/150‑1 (MQ=255)
                                                    aGGCATGGAAACGTTTAAATTTATACAGCAGTTATTTTTTTTCCTGATTTTAAAGCGGGCTTTCCTTCAGTAACGGGTGGGGGGTAGCGTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCg                                                                                            <  1:730842/148‑1 (MQ=255)
                                                     ggCGTGAAAAGGTTTGATTTTATACAGCAGTTATTTTTTTTCCTGATTTTTAACCGGGCTTCCCGTCAGAACCGGGGGGGGGGTGGCCGAAATAACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAAcc                                                                                        <  2:960691/151‑1 (MQ=255)
                                                     ggCCTGAAAAGGTTTAATTTTATACAGGAGTTAATTTTTTTCCTGATTTTAAAGGGGGCTTTCCTTAAGTACCCGGGGGGGGTTAGCCTAATTTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGaa                                                                                          <  1:1185455/149‑1 (MQ=255)
                                                      gCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCt                                                                                       >  2:796348/1‑151 (MQ=255)
                                                            aaaCGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCg                                                                                 >  2:109135/1‑151 (MQ=255)
                                                               agggttttattttATACAGAGGTTATTTTTTTCCAGTTTTTTTAAGCGGGCCTTCCTAAAGAACCCGGTGGGGGGCTCCAGAAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGt                                                                                <  1:38146/141‑1 (MQ=255)
                                                                cGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGc                                                                              >  2:184984/1‑151 (MQ=255)
                                                                    tGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCaa                                                                           >  1:468499/1‑150 (MQ=255)
                                                                      aaTTTATACAGCAGTTATTTTTTTTCCTGTTGTTAAAAGGGGCGTCCCGTCAGAAACCAGTGGGCGGTACCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAAt                                                                        <  2:632578/151‑1 (MQ=255)
                                                                          ttATACAGCAGGTATTTTTTTTCTTGATTTTAAAGCGGGCGTCCCGTCAGAACCCAGGGGGGGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGa                                                                      <  1:202555/150‑1 (MQ=255)
                                                                              aCAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACt                                                                 >  1:29146/1‑151 (MQ=255)
                                                                                   agGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTgggg                                                             >  1:487993/1‑150 (MQ=255)
                                                                                    gttattttttttCTTGATGTTAAAGCGGGCGTCCCGTCAGAAACCGGGGGGGGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  2:506098/144‑1 (MQ=255)
                                                                                    ggTAATTTTTTTCCTGATGTTTAAGGGGGCTTTCGGTAAGTACCAGGTGGGCGGTCGCTTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  1:1095635/151‑1 (MQ=255)
                                                                                      tAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGa                                                          >  1:1235974/1‑150 (MQ=255)
                                                                                           tttttCCTGTTGTTAAACCGGGCGTCCCGTAAGTACCCGGGGGGGGGTGGCTTAAATTACTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCg                                                           <  1:103735/144‑1 (MQ=255)
                                                                                           tttttCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATcc                                                    >  1:740720/1‑151 (MQ=255)
                                                                                              ttCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATccc                                                   >  2:684422/1‑149 (MQ=255)
                                                                                                ccTGTTGTTAAAGGGGGCGTCCCGTCAGTACCCGGGGGGCGTTACCTTAAATAAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGt                                                <  1:964002/150‑1 (MQ=255)
                                                                                                ccTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCCGTg                                               >  1:1145283/1‑151 (MQ=255)
                                                                                                   tGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              >  2:610177/1‑150 (MQ=255)
                                                                                                       tttttAAGCGGGCGTTCCGTAAGTAACGGGTGGGGGGTACCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGc                                              <  1:1259723/145‑1 (MQ=255)
                                                                                                        ttttAAACCGGGCTTCCCGTCAGTAACCGGGGGGCGGTAGCCTAAATTAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAg                                                 <  1:65046/143‑1 (MQ=255)
                                                                                                        ttttAAACCGGGCGTTCCGTAAGTAACCGGGGGGGGGGAGCCAGAATTACCCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAg                                            <  2:360338/148‑1 (MQ=255)
                                                                                                        ttttAAACGGGGCGTCCCGTAAGAACCCGGGGGGGGGTAGCCTAAATAAACCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTg                                               <  2:468497/145‑1 (MQ=255)
                                                                                                           tAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCttt                                      >  1:1059203/1‑151 (MQ=255)
                                                                                                             aagcGGGCGTCCCGTCAGAACCGGGGGGGCGGTAGCCTAAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCt                                        <  1:159217/145‑1 (MQ=255)
                                                                                                             aaCGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTa                                     >  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)
                                                                                                                                                    |                                                                                                                                                   
CGACGATCTCCAGTGTGCGGATCGGTGATGACGCCTTCAAGCCCATAGCGACAGGCCTGGAAACGGTTGAATTTATACAGCAGGTAATCTTTTTCCTGATGTTTAAACGGGCGTTCCGTCAGTAACCAGTGGGCGGTAGCCTGAATTAATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCG  >  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: