Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,961 A→T 6.1% N298I (AAT→ATT)  ynfE → putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,9610AT 6.1% 256.7 / 2.8 79N298I (AAT→ATT) ynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base A (37/37);  new base T (2/3);  total (39/40)
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.00e+00

AGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCT  >  NC_000913/1658812‑1659106
                                                                                                                                                     |                                                                                                                                                 
aGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGaaaa                                                                                                                                                   >  2:897838/1‑150 (MQ=255)
   gcgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtct                                                                                                                                                >  2:1400419/1‑150 (MQ=255)
    cgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                                                                                                                                               >  2:1063233/1‑150 (MQ=255)
    cgcgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAAtctc                                                                                                                                               >  1:170363/1‑150 (MQ=255)
      cgAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGt                                                                                                                                             >  1:1133642/1‑150 (MQ=255)
       gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCTATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGtt                                                                                                                                            <  1:863127/150‑1 (MQ=255)
       gagaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATCAACGAGATTTTCGtt                                                                                                                                            >  1:908132/1‑150 (MQ=255)
        agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTATCTAAAATATCTTTg                                                                                                                                           <  2:895798/150‑1 (MQ=255)
        agaAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTg                                                                                                                                           <  2:1004679/150‑1 (MQ=255)
           aaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATc                                                                                                                                        >  2:1260493/1‑150 (MQ=255)
           aaTCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATCAACGAGATTTTCGTTAATc                                                                                                                                        <  1:1633905/150‑1 (MQ=255)
                 aCGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATCAACGAGATTTTCGTTAATCaatacc                                                                                                                                  <  2:908132/150‑5 (MQ=255)
                   gCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATCAACGAGATTTTCGTTAATCaataccca                                                                                                                                <  2:1120883/150‑7 (MQ=255)
                    ccAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtc                                                                                                                               >  1:486986/1‑150 (MQ=255)
                     cAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtct                                                                                                                              >  2:490650/1‑150 (MQ=255)
                      aaaaTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTtctc                                                                                                                             <  2:215172/150‑1 (MQ=255)
                         aTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGAt                                                                                                                          <  2:1308572/150‑1 (MQ=255)
                         aTGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGAt                                                                                                                          <  1:454502/150‑1 (MQ=255)
                          tGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTATCGAAAATCTCGTTGATCAACCTTTTCTCGATa                                                                                                                         <  2:1427108/150‑1 (MQ=255)
                            aTTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATaaa                                                                                                                       >  1:1320909/1‑150 (MQ=255)
                                 tATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTg                                                                                                                  >  2:997401/1‑150 (MQ=255)
                                 tATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTg                                                                                                                  >  1:1602244/1‑150 (MQ=255)
                                     gATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTc                                                                                                              >  1:1329630/1‑150 (MQ=255)
                                           cgATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTAt                                                                                                        <  1:1198848/150‑1 (MQ=255)
                                            gATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATg                                                                                                       <  1:897838/150‑1 (MQ=255)
                                              tataCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGAc                                                                                                     <  1:589859/150‑1 (MQ=255)
                                              aataCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTTTTGATAAACGAAAATCTCGTTGATCAAGCTTTTCTCGTTAAATACTGCGTCGGCTATGAc                                                                                                     <  2:1068562/149‑1 (MQ=255)
                                                taCCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGa                                                                                                   >  1:777369/1‑150 (MQ=255)
                                                  ccGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGaaa                                                                                                 >  1:853538/1‑150 (MQ=255)
                                                         ggCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCtt                                                                                          >  1:885095/1‑150 (MQ=255)
                                                            tGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTAcc                                                                                       >  1:6436/1‑150 (MQ=255)
                                                               cGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGc                                                                                    <  2:1133642/150‑1 (MQ=255)
                                                                ggTCGTGAAGACGAATGGCTCCCTATTCGACCGGGCACCGATGCCGCGCTGGTTGCGGGTATTCCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTCCGTCGCCTATGACGAAAAAACCTTACCCgaa                                                                                   >  2:64867/1‑148 (MQ=255)
                                                                  tCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGa                                                                                 <  1:490650/150‑1 (MQ=255)
                                                                   cGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGAt                                                                                >  2:696916/1‑150 (MQ=255)
                                                                    gTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATg                                                                               >  2:254987/1‑150 (MQ=255)
                                                                    gTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATg                                                                               >  1:5328/1‑150 (MQ=255)
                                                                     tGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGc                                                                              >  2:471548/1‑150 (MQ=255)
                                                                                ggCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCAACCAAAAATg                                                                   >  2:1219281/1‑150 (MQ=255)
                                                                                 gCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATgg                                                                  >  1:514677/1‑150 (MQ=255)
                                                                                   tCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTc                                                                <  1:64867/150‑1 (MQ=255)
                                                                                    cccTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCa                                                               >  2:1390466/1‑150 (MQ=255)
                                                                                     ccTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCAc                                                              <  1:301899/150‑1 (MQ=255)
                                                                                      cTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACt                                                             <  2:77220/150‑1 (MQ=255)
                                                                                       tATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACta                                                            <  1:1390466/150‑1 (MQ=255)
                                                                                           cGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATaaa                                                        >  2:896836/1‑150 (MQ=255)
                                                                                            gCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAg                                                       <  1:1176902/150‑1 (MQ=255)
                                                                                             cccGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGc                                                      <  1:58742/150‑1 (MQ=255)
                                                                                              ccGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGcc                                                     <  1:1400419/150‑1 (MQ=255)
                                                                                                gggCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCta                                                   <  2:1320909/150‑1 (MQ=255)
                                                                                                 ggCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCtat                                                  >  2:648062/1‑150 (MQ=255)
                                                                                                     ccGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTc                                              >  1:514104/1‑150 (MQ=255)
                                                                                                         tGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTgg                                          <  1:1219281/150‑1 (MQ=255)
                                                                                                         tGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTgg                                          >  2:13803/1‑150 (MQ=255)
                                                                                                          gCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTATTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACACAAAAATGGTCACTATAAACCCTATATTCTTGGt                                         >  2:447140/1‑150 (MQ=255)
                                                                                                            cgcgCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGa                                       <  1:1173645/150‑1 (MQ=255)
                                                                                                               gcTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAgg                                    <  1:1517524/150‑1 (MQ=255)
                                                                                                                   gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGAc                                >  1:759995/1‑150 (MQ=255)
                                                                                                                   gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGAc                                <  1:1260493/150‑1 (MQ=255)
                                                                                                                   gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGAc                                <  1:896836/150‑1 (MQ=255)
                                                                                                                      gCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGAt                             <  2:885095/150‑1 (MQ=255)
                                                                                                                        gggTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATaa                           <  1:63261/150‑1 (MQ=255)
                                                                                                                         ggTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATaaa                          <  1:648062/150‑1 (MQ=255)
                                                                                                                           tATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAAc                        <  2:1252159/150‑1 (MQ=255)
                                                                                                                              tGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCa                                                               <  1:1636268/108‑1 (MQ=255)
                                                                                                                              tGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCa                                                               >  2:1636268/1‑108 (MQ=255)
                                                                                                                                ccTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGa                   >  1:1065303/1‑150 (MQ=255)
                                                                                                                                ccTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGa                   >  1:237843/1‑150 (MQ=255)
                                                                                                                                 caaaaggttGATCAACGAGATTTTCGTTAATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGaa                  >  1:1347138/8‑150 (MQ=255)
                                                                                                                                      tATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcgc             <  2:170363/150‑1 (MQ=255)
                                                                                                                                       aTTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcgcc            <  2:634526/150‑1 (MQ=255)
                                                                                                                                         tGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcgccgc          <  1:756422/150‑1 (MQ=255)
                                                                                                                                         tGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcgccgc          >  1:62066/1‑150 (MQ=255)
                                                                                                                                             tAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTg      >  1:755788/1‑150 (MQ=255)
                                                                                                                                              aaCGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTgg     >  2:1095811/1‑150 (MQ=255)
                                                                                                                                               aCGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAACGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTggg    >  1:1665755/1‑150 (MQ=255)
                                                                                                                                               aCGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTggg    >  1:144805/1‑150 (MQ=255)
                                                                                                                                                cGAAAATCTCGTTGATGAACCTTTTCTCTATATATTCTGCGTCGGGTATGTCGTTATTTCCTTACCCGCTGTTGCTCCCATAAATGGTCTCTATTAAGCCTATATTCTTGGTGAAGTTGACGATAATTCAGCGAATAGGCCTCAGTGGGc   <  2:853538/150‑1 (MQ=255)
                                                                                                                                                cGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGc   <  1:997401/150‑1 (MQ=255)
                                                                                                                                                 gAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCt  <  1:471548/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
AGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCT  >  NC_000913/1658812‑1659106

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