Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,658,953 T→C 6.7% I295I (ATT→ATC ynfE → putative selenate reductase YnfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,658,9530TC 6.7% 227.5 / 5.2 78I295I (ATT→ATCynfEputative selenate reductase YnfE
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base C (2/3);  total (38/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 = 7.84e-01

CTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTG  >  NC_000913/1658804‑1659102
                                                                                                                                                     |                                                                                                                                                     
cttGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGAtt                                                                                                                                                       >  1:895798/1‑150 (MQ=255)
   gAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAAc                                                                                                                                                    >  2:875688/1‑150 (MQ=255)
   gAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAAc                                                                                                                                                    >  2:756422/1‑150 (MQ=255)
   gAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAAc                                                                                                                                                    <  1:411495/150‑1 (MQ=255)
   gAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAAc                                                                                                                                                    <  2:632998/150‑1 (MQ=255)
     aaaaGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGa                                                                                                                                                  <  1:52243/150‑1 (MQ=255)
        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                                                                           >  1:5328/1‑150 (MQ=255)
                                                                            gTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATg                                                                           >  2:254987/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                                      >  2:13803/1‑150 (MQ=255)
                                                                                                                 tGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTgg                                      <  1:1219281/150‑1 (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                                                           >  2:1636268/1‑108 (MQ=255)
                                                                                                                                      tGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCa                                                           <  1:1636268/108‑1 (MQ=255)
                                                                                                                                        ccTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGa               >  1:237843/1‑150 (MQ=255)
                                                                                                                                        ccTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGa               >  1:1065303/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)
                                                                                                                                                     |                                                                                                                                                     
CTTGAAAAAGCGCGCGAGAAATCGAACGCCAAAATGATTGTTATCGATCCGCGATATACCGATACGGCTGCCGGTCGTGAAGACGAATGGCTCCCTATTCGCCCGGGCACCGATGCCGCGCTGGTTGCGGGTATTGCCTGGGTATTGATTAACGAAAATCTCGTTGATCAACCTTTTCTCGATAAATACTGCGTCGGCTATGACGAAAAAACCTTACCCGCAGATGCACCCAAAAATGGTCACTATAAAGCCTATATTCTTGGTGAAGGTGACGATAAAACAGCGAAAACGCCGCAGTG  >  NC_000913/1658804‑1659102

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