Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,421 T→G 5.5% intergenic (+21/+223) intA → / ← yfjH CP4‑57 prophage; integrase/CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,4210TG 5.5% 323.5 / 8.7 109intergenic (+21/+223)intA/yfjHCP4‑57 prophage; integrase/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (51/52);  new base G (2/4);  total (53/56)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

Only 100 of 114 total aligned reads displayed.

GCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAA  >  NC_000913/2757272‑2757570
                                                                                                                                                     |                                                                                                                                                     
gcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctct                                                                                                                                                       >  1:703694/1‑150 (MQ=255)
   tataACCGCAGTGTTTATCTTGAGCAACGTGGTCCGATTATGCAATGGTGGGCTGATATCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTTAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCTGCTCTCTg                                                                                                                                                    <  2:131836/150‑1 (MQ=255)
         cGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGa                                                                                                                                              <  1:1459433/150‑1 (MQ=255)
         cGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGa                                                                                                                                              <  1:1366233/150‑1 (MQ=255)
         cGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGa                                                                                                                                              <  2:1394799/150‑1 (MQ=255)
          gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTATCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGat                                                                                                                                             <  2:294178/150‑1 (MQ=255)
            aGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGatat                                                                                                                                           >  2:618844/1‑150 (MQ=255)
               gATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATAttta                                                                                                                                        >  1:1016776/1‑149 (MQ=255)
                 ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                      <  2:1702511/150‑1 (MQ=255)
                   aTCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaaa                                                                                                                                    >  2:762361/1‑150 (MQ=255)
                    tCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaaaa                                                                                                                                   >  2:106702/1‑150 (MQ=255)
                      ttGAGCAAAGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGa                                                                                                                                 >  1:1340158/1‑150 (MQ=255)
                       tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                                <  2:416831/150‑1 (MQ=255)
                       tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                                <  1:676429/150‑1 (MQ=255)
                         aGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCg                                                                                                                              <  1:153525/150‑1 (MQ=255)
                          gCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGTATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                             <  1:751552/150‑1 (MQ=255)
                          gCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                             <  2:1044903/150‑1 (MQ=255)
                            aaCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                           >  1:183174/1‑150 (MQ=255)
                              cgtcgtCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTa                                                                                                                         <  2:541577/150‑1 (MQ=255)
                                  gtcCTATGATGCAATGGTGGTCTGATCTCGTAAAAGCAGCAGATAGTTGTAGCATCGTTTTTACTCATTTGAGCTAAATTCGTCTTGTCGGATAAAAAATTTATAAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATa                                                                                                                     <  1:414666/150‑1 (MQ=255)
                                  gtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATa                                                                                                                     >  1:1625850/1‑150 (MQ=255)
                                   tcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAg                                                                                                                    >  2:1088772/1‑150 (MQ=255)
                                     cGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTATCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCt                                                                                                                  >  2:1051959/1‑150 (MQ=255)
                                      gatgatGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTa                                                                                                                 >  2:481369/1‑150 (MQ=255)
                                          atgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcgcg                                                                                                             <  2:400386/150‑1 (MQ=255)
                                          atgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcgcg                                                                                                             <  1:1706815/150‑1 (MQ=255)
                                           tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                            <  2:919054/150‑1 (MQ=255)
                                             cAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCATCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATATATACGCGTaa                                                                                                          >  1:350851/1‑150 (MQ=255)
                                                 ggtggGCTGATCTCGTAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGc                                                                                                                                                 >  2:1682431/1‑106 (MQ=255)
                                                 ggtggGCTGATCTCGTAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGc                                                                                                                                                 <  1:1682431/106‑1 (MQ=255)
                                                  gtggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTg                                                                                                     >  1:1075593/1‑150 (MQ=255)
                                                   tggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGc                                                                                                    >  1:69483/1‑150 (MQ=255)
                                                   ggggCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGc                                                                                                    <  2:493316/149‑1 (MQ=255)
                                                        tGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTAc                                                                                               <  2:1340158/150‑1 (MQ=255)
                                                        tGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTAc                                                                                               <  1:133941/150‑1 (MQ=255)
                                                         gATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACt                                                                                              <  2:369650/150‑1 (MQ=255)
                                                         gATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACt                                                                                              >  1:1403150/1‑150 (MQ=255)
                                                             tcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCg                                                                                          <  1:594400/150‑1 (MQ=255)
                                                             tcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATAATTCTAAAAGATCGCAATAAAAAGCTACGCTGAATAGTGCTCTACTTTct                                                                                          >  2:848158/1‑149 (MQ=255)
                                                              cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCATCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                         >  1:1496641/1‑150 (MQ=255)
                                                               gTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGtt                                                                                        <  1:1486141/150‑1 (MQ=255)
                                                                tAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGttt                                                                                       <  2:714938/150‑1 (MQ=255)
                                                                   aaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGc                                                                                    <  2:1403150/150‑1 (MQ=255)
                                                                   aaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGc                                                                                    <  1:106702/150‑1 (MQ=255)
                                                                   aaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGc                                                                                    <  1:660402/150‑1 (MQ=255)
                                                                    agcagcAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCTGATAAAATATTTATCAGCCGAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGcc                                                                                   <  1:210190/150‑1 (MQ=255)
                                                                    agcagcAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGcc                                                                                   <  1:1442315/150‑1 (MQ=255)
                                                                    agcagcAGATAGTGGTAGCATCGTTTTAACTCAATTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATGTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGcc                                                                                   >  1:1055946/1‑150 (MQ=255)
                                                                       agcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCtt                                                                                >  1:1664106/1‑150 (MQ=255)
                                                                       agcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCtt                                                                                >  2:817559/1‑150 (MQ=255)
                                                                        gcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATAATTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCttt                                                                               >  2:571957/1‑150 (MQ=255)
                                                                          agATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCa                                                                             >  2:213126/1‑150 (MQ=255)
                                                                            aTAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCaga                                                                           <  2:165167/150‑1 (MQ=255)
                                                                            aTAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCaga                                                                           <  2:1316711/150‑1 (MQ=255)
                                                                               gTGGTTGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGt                                                                        <  2:223571/150‑1 (MQ=255)
                                                                                     gCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACa                                                                  >  1:44773/1‑150 (MQ=255)
                                                                                     gCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACa                                                                  >  2:765207/1‑150 (MQ=255)
                                                                                     gCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGATACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACa                                                                  >  1:556640/1‑150 (MQ=255)
                                                                                      cATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACAt                                                                 >  2:886295/1‑150 (MQ=255)
                                                                                          gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAATATATTTATCAGCCCAGCTCTCTGTGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                             <  1:74951/150‑1 (MQ=255)
                                                                                          gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                             <  2:556669/150‑1 (MQ=255)
                                                                                          gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                             >  1:1072724/1‑150 (MQ=255)
                                                                                            tttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGc                                                           <  2:703694/150‑1 (MQ=255)
                                                                                              tAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCaa                                                         >  1:801310/1‑150 (MQ=255)
                                                                                                 cTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAAAGGTAAGACATTAGGGCAAttt                                                      >  2:925515/1‑150 (MQ=255)
                                                                                                    aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGATTAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                   >  2:508783/1‑150 (MQ=255)
                                                                                                    aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                   <  2:1589291/150‑1 (MQ=255)
                                                                                                      ttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                                 <  2:357825/150‑1 (MQ=255)
                                                                                                      ttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                                 <  2:1664106/150‑1 (MQ=255)
                                                                                                       tGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGt                                                >  1:473464/1‑150 (MQ=255)
                                                                                                          gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                             >  1:1000984/1‑150 (MQ=255)
                                                                                                          gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                             >  1:911157/1‑150 (MQ=255)
                                                                                                          gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                             <  1:213126/150‑1 (MQ=255)
                                                                                                           cAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGa                                            >  1:729398/1‑150 (MQ=255)
                                                                                                             aaaTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACt                                          >  2:1464020/1‑150 (MQ=255)
                                                                                                                ttCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTgcg                                       <  1:896526/150‑1 (MQ=255)
                                                                                                                 tCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAAAAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTgcgc                                      <  2:172274/150‑1 (MQ=255)
                                                                                                                  cGTCTTTTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCa                                     <  1:426452/150‑1 (MQ=255)
                                                                                                                  cGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCa                                     >  1:104771/1‑150 (MQ=255)
                                                                                                                  cGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCAGGGCTGATATTTCTAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAg                                    >  2:326377/1‑150 (MQ=255)
                                                                                                                      ttGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCa                                 >  1:1278596/1‑150 (MQ=255)
                                                                                                                         tCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCg                              <  1:762361/150‑1 (MQ=255)
                                                                                                                           ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacag                            <  1:817559/150‑5 (MQ=255)
                                                                                                                           ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacag                            <  2:1552286/150‑5 (MQ=255)
                                                                                                                            gATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATTGAGt                           <  1:1034836/150‑1 (MQ=255)
                                                                                                                             aTAAAATATTTATCAGCCCATCTCTCTGGGCTTATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTAAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacagtg                          <  2:686415/150‑7 (MQ=255)
                                                                                                                                  agaaatATCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgg                     >  2:1231718/6‑150 (MQ=255)
                                                                                                                                    ttttATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGATTAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgggg                   <  1:508783/149‑1 (MQ=255)
                                                                                                                                        aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt               >  2:927470/1‑150 (MQ=255)
                                                                                                                                        aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt               >  1:890915/1‑150 (MQ=255)
                                                                                                                                        aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAAAAGTGCTCTACTTTCGTTTAGCCCTTTCCGAGGTAAGACATTAGGGCAAATTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt               <  1:202196/150‑1 (MQ=255)
                                                                                                                                         tCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTg              >  1:313840/1‑150 (MQ=255)
                                                                                                                                          cAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGc             >  2:679371/1‑150 (MQ=255)
                                                                                                                                           aGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCt            >  2:867105/1‑150 (MQ=255)
                                                                                                                                           aGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCt            >  2:1475692/1‑150 (MQ=255)
                                                                                                                                            gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTc           <  2:1072724/150‑1 (MQ=255)
                                                                                                                                            gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTAACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCta           >  2:77545/1‑149 (MQ=255)
                                                                                                                                                  ctctctGGGATGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg     <  2:352257/150‑1 (MQ=255)
                                                                                                                                                    ctctGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCa   <  1:775699/150‑1 (MQ=255)
                                                                                                                                                     tctGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCaa  <  1:177903/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTGCAA  >  NC_000913/2757272‑2757570

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