Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 8.8% pseudogene (128/129 nt) ykgQ → pseudogene, putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5660GA 8.8% 160.5 / 5.5 57pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (2/3);  total (27/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.02e-01

AGTCCTACATTATTAATGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAGAGTATCTTTCAATACCAGGCGACCAGCACCGG  >  NC_000913/317273‑317859
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     
aGTCCTACATTATTAATGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                                                                                                                                                                                                                                >  2:162096/1‑301 (MQ=255)
         ttattaATGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                                                                                                                                                                                                                       >  2:409892/1‑301 (MQ=255)
           tttaATGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttt                                                                                                                                                                                                                                                                                     <  2:259152/300‑1 (MQ=255)
              aaTGTATTTTTTGTATAAGTGATTTCAAAATCCACATACTACTTTGAGGTTACGGTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                                                                                                                                                                                                                  <  1:99336/301‑1 (MQ=255)
               tttaaTTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                                                                                                                                                                                                                                 <  2:110821/298‑5 (MQ=255)
                tGAATTTTTTGCATAAGTGTTATCTATTGCCTCATACTACTTTTAGGTTACGTTTTAACTTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAAGCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacag                                                                                                                                                                                                                                                                                <  1:572977/301‑6 (MQ=255)
                tGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                                                >  2:200356/1‑301 (MQ=255)
                ataaTTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                                                                                                                                                                                                                                <  1:618369/299‑1 (MQ=255)
                  atctttttGCATTAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACGGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAg                                                                                                                                                                                                                                                                              <  1:575004/298‑1 (MQ=255)
                     tttttGCTGTAGTGATATCAAAATATACATACTACTTTTAGGTGTCGTTTTTACGTTGACTCATTGTTCATGCCTTATGGTGGGACTGTCAGTTTAATTGTCCGTATAGTCTGCTTGTGAGTATTGCAGATGCCCATCAAAAAGCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAAGCCTGAATTTGATCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTGACCTAAgacg                                                                                                                                                                                                                                                                           <  1:57284/301‑1 (MQ=255)
                      ttttGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgagat                                                                                                                                                                                                                                                                          >  2:575719/1‑298 (MQ=255)
                                  aatatCAAATTCCACATTCTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                                                                                                                                                                                                                              <  1:310471/300‑1 (MQ=255)
                                       cAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACaa                                                                                                                                                                                                                                                         >  2:574090/1‑300 (MQ=255)
                                              cacaTACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACAACAAAAACCTTTAAACAGGGAGttt                                                                                                                                                                                                                                                  >  1:631121/1‑299 (MQ=255)
                                                      actTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCGTCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                                                                                                                                                                                                                          >  2:539644/1‑301 (MQ=255)
                                                       ctTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTAGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAGAGGGAGTGATATTGATCt                                                                                                                                                                                                                                         >  1:565528/1‑301 (MQ=255)
                                                                 tcGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                                                                                                                                                               <  2:59356/300‑1 (MQ=255)
                                                                 aCGTTTTAACTTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                                                                                                                                                                                                                               <  2:570737/301‑1 (MQ=255)
                                                                      ttctcGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATCGATCTTCACTCTGTCATAtc                                                                                                                                                                                                                          <  2:148161/297‑1 (MQ=255)
                                                                          cGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAATGATATCGATCTTCACTATGTCATAtcagat                                                                                                                                                                                                                      >  1:148035/1‑297 (MQ=255)
                                                                                   tttGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAGAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                                                                                                                                                                             <  2:566242/300‑1 (MQ=255)
                                                                                   aTTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg                                                                                                                                                                                                             >  2:71405/1‑301 (MQ=255)
                                                                                            gCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTACTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt                                                                                                                                                                                                    >  2:477684/1‑301 (MQ=255)
                                                                                              cTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTtc                                                                                                                                                                                                  >  2:34109/1‑301 (MQ=255)
                                                                                                            gACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAgag                                                                                                                                                                                    >  2:222717/1‑299 (MQ=255)
                                                                                                                gTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGcc                                                                                                                                                                                >  2:277483/1‑301 (MQ=255)
                                                                                                                   gAATTGTCCGTATAGGCTTCTTCAGAGTATTGCAGATTCCCATCTAATATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGa                                                                                                                                                                             <  1:200113/301‑1 (MQ=255)
                                                                                                                                     gCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGACCTGACGCCTGAAAAAGTGaaca                                                                                                                                                           >  2:94576/1‑301 (MQ=255)
                                                                                                                                           gagTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac                                                                                                                                                     >  1:156959/1‑301 (MQ=255)
                                                                                                                                           aagTATTGCAGTTGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagac                                                                                                                                                     <  1:71315/300‑1 (MQ=255)
                                                                                                                                             gTATTGCAGTTGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAacagacag                                                                                                                                                   <  1:161945/301‑1 (MQ=255)
                                                                                                                                                  gCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTc                                                                                                                                              <  1:409329/301‑1 (MQ=255)
                                                                                                                                                      aTGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGaacactgtctgttgttcactt                                                                                                                                          <  1:34058/301‑18 (MQ=255)
                                                                                                                                                                      gATTTCTTACTGACAGTTAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGAGTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGttt                                                                                                                          <  1:538971/301‑1 (MQ=255)
                                                                                                                                                                                 gACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGaa                                                                                                               >  2:61732/1‑301 (MQ=255)
                                                                                                                                                                                  aCAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAAc                                                                                                              <  2:432049/301‑1 (MQ=255)
                                                                                                                                                                                       tAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGaa                                                                                                         >  2:502531/1‑300 (MQ=255)
                                                                                                                                                                                          aGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCa                                                                                                      >  1:547211/1‑301 (MQ=255)
                                                                                                                                                                                                gggacgtaTTTGTTTGTTTGATTTGGGTTTTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTTTCTTGCTGTAAAAATACATGAGTGATATTAATCATGTATTGTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTg                                                                                                <  1:474185/294‑1 (MQ=255)
                                                                                                                                                                                                           gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAg                                                                                     <  1:156360/301‑1 (MQ=255)
                                                                                                                                                                                                            ttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGa                                                                                    >  2:343850/1‑301 (MQ=255)
                                                                                                                                                                                                              cGCTTGTTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTTTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAAt                                                                                  <  1:61664/301‑1 (MQ=255)
                                                                                                                                                                                                                  tGATTTGGGTTATGTTGGGGAGTCGGTAATACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTcc                                                                              <  1:343387/301‑1 (MQ=255)
                                                                                                                                                                                                                                tttgggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGc                                                                <  1:222474/298‑1 (MQ=255)
                                                                                                                                                                                                                                    tagacgtgtgctcttccgatttCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCAtt                                                            <  2:616995/280‑1 (MQ=255)
                                                                                                                                                                                                                                        tCGGTTATACATTTTAATCCTGTATTTGTTCGTTGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCa                                                        <  2:547884/301‑1 (MQ=255)
                                                                                                                                                                                                                                           ttttgctcttccgaTCTTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGa                                                     <  2:478093/288‑1 (MQ=255)
                                                                                                                                                                                                                                                      tAATCCTTAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGcc                                          <  2:631998/301‑1 (MQ=255)
                                                                                                                                                                                                                                                         tCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAACAATTTCCCGCACCGAATCTGCCATTAACATGACTGAACCTGCCAGc                                       >  1:572884/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTaaatcggaagagcgacgtgtag                                      >  1:616137/1‑280 (MQ=255)
                                                                                                                                                                                                                                                            tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGaagatcgaaagagcgtcg                                    >  1:477492/1‑284 (MQ=255)
                                                                                                                                                                                                                                                                   gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGACAAAGTGAACAACAGACAGAGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATACATTACTGCGTTGTTTAAGAAATTCCCGATCCGAATCTGCCATTATCCTAACTTCACCTACCAGCATCAGAGTaa                             >  2:323042/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         ggtttttttgatagaaaaaTTCATGATTGATATTAATTATGTATTTTACCTAAGACGAAAAAAACCTTTTAACAGGGAGTTATATTGATCTTCATTCTGTCATATCTCCGGTAATATTGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAGAGTATCTTTCa                       <  1:322638/287‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    aggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCAAGCATCAGAATATCTTTCAATACCAGTCgc            >  1:316984/5‑300 (MQ=255)
                                                                                                                                                                                                                                                                                      tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAGAGTATCTTTCAATACCAGGCGAcc          <  2:157110/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                        aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCAAGCAACAGAATATCTTTCAAAACCAGGcaccccc        >  1:474669/1‑295 (MQ=255)
                                                                                                                                                                                                                                                                                              cATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAGAGTATCTTTCAATACCAGGCGACCAACACCga  >  2:80535/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                     
AGTCCTACATTATTAATGAATTTTTTGCATAAGTGATATCAAAATCCACATACTACTTTGAGGTTACGTTTTAACGTAGACTCATTGTTCATGCCTAATGGAGGGACTGACAGTGGAATTTTCCGTATAGTCTGCTTCAGAGTATTGCAGATGCCCATCAAAAATCGATTTCTCACTGACAGATAAAGAAATGGCACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCTTCATCATCAGAAAGGGGGCCTGACGCCTGAAAAAGTGAACAACAGACAGTGTTCGGATTATCAATATTCAGTTTTTAACGTTGAACTGGATTCATTACTGCGTTGTTTAAGAATTTCCCGCGCCGAATCTGCCATTATCATGACTGCACCTGCCAGCATCAGAGTATCTTTCAATACCAGGCGACCAGCACCGG  >  NC_000913/317273‑317859

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: