Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,963 C→T 6.6% L299F (CTC→TTC)  ynfE → putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,9630CT 6.6% 233.5 / 4.2 76L299F (CTC→TTC) ynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (2/3);  total (36/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

GCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAA  >  NC_000913/1658815‑1659112
                                                                                                                                                    |                                                                                                                                                     
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:896836/150‑1 (MQ=255)
                                                                                                                gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGAc                                      <  1:1260493/150‑1 (MQ=255)
                                                                                                                gTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGAc                                      >  1:759995/1‑150 (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:62066/1‑150 (MQ=255)
                                                                                                                                      tGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAAcgccgc                <  1:756422/150‑1 (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)
                                                                                                                                                    ctcGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCaa  >  2:410990/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTGGGCTTCGCAA  >  NC_000913/1658815‑1659112

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: