Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,418 C→A 7.1% intergenic (+18/+226) 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,4180CA 7.1% 363.0 / 15.5 127intergenic (+18/+226)intA/yfjHCP4‑57 prophage; integrase/CP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (58/60);  new base A (4/5);  total (62/65)
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.20e-01

Only 100 of 132 total aligned reads displayed.

GCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757270‑2757567
                                                                                                                                                    |                                                                                                                                                     
gcgcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGct                                                                                                                                                      >  1:1036878/1‑150 (MQ=255)
 cgcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGATCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctc                                                                                                                                                     >  2:1178059/1‑150 (MQ=255)
 cgcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctc                                                                                                                                                     <  1:993346/150‑1 (MQ=255)
       taACCGCAGTGATTATCTTGAGCAACGTCGTCCGTTGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTggg                                                                                                                                               <  1:1403843/150‑1 (MQ=255)
        aaCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGc                                                                                                                                              >  2:1132462/1‑150 (MQ=255)
        aaCCGCAGAGATTATCTTGAGCAACGTCGTCCGATGATGCAATCGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTAGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTcggg                                                                                                                                              >  2:255922/1‑149 (MQ=255)
         aCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCt                                                                                                                                             >  1:1049939/1‑150 (MQ=255)
         aCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCt                                                                                                                                             >  2:1561344/1‑150 (MQ=255)
          ccGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTg                                                                                                                                            >  1:777281/1‑150 (MQ=255)
          ccGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTg                                                                                                                                            <  2:1176544/150‑1 (MQ=255)
            gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTTGGCTGATCTCGTAAAAGCAGCAGATAGTGGTATCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGat                                                                                                                                          <  1:807390/150‑1 (MQ=255)
            gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGat                                                                                                                                          <  1:243656/150‑1 (MQ=255)
             cAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGata                                                                                                                                         >  1:432752/1‑150 (MQ=255)
               gTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGTCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGataaa                                                                                                                                       <  1:1237698/150‑3 (MQ=255)
                tGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGTAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGataaat                                                                                                                                      <  1:278616/150‑4 (MQ=255)
                    tATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaa                                                                                                                                  <  2:1451146/150‑1 (MQ=255)
                     aTCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaaa                                                                                                                                 >  2:1423486/1‑150 (MQ=255)
                        ttGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCCGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGa                                                                                                                              >  2:852821/1‑150 (MQ=255)
                        ttGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGa                                                                                                                              >  1:679257/1‑150 (MQ=255)
                         tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             >  2:1388878/1‑150 (MQ=255)
                         tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             >  2:627952/1‑150 (MQ=255)
                            gCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          <  2:1675257/150‑1 (MQ=255)
                             cAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCa                                                                                                                         >  2:1402359/1‑150 (MQ=255)
                                cgtcgtCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTa                                                                                                                      <  2:807880/150‑1 (MQ=255)
                                          tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  1:385096/150‑1 (MQ=255)
                                          tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  2:154522/150‑1 (MQ=255)
                                             tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTCATATTTCTAAAATATCTCACTAAATAGCTACGCgg                                                                                                         >  2:1612078/1‑150 (MQ=255)
                                                  tggtggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAg                                                                                                    <  1:806816/150‑1 (MQ=255)
                                                     tggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCTCTAAATAGCTACGCGGAATAGTGc                                                                                                 >  1:637440/1‑150 (MQ=255)
                                                     tggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCTCTAAATAGCTACGCGGAATAGTGc                                                                                                 <  2:637440/150‑1 (MQ=255)
                                                       ggCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGAAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTATAAAAGATAGAACTAAATAGCTACGCGGAATAGTGctc                                                                                               >  1:1302396/1‑150 (MQ=255)
                                                         cTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCATCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTa                                                                                             >  1:1228297/1‑150 (MQ=255)
                                                         cTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTATAAAAGATCGCAATAAATAGCTACGCGGAATAGTGCtata                                                                                             >  1:167026/1‑150 (MQ=255)
                                                         cTGATCTCGAAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTTAAAGATCGCACTAAATAGATAAGCGGAATAGTGCTCTa                                                                                             >  1:105176/1‑150 (MQ=255)
                                                            aTCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGATCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACtt                                                                                          <  2:234813/150‑1 (MQ=255)
                                                            aTCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGAAAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACtt                                                                                          <  1:257771/150‑1 (MQ=255)
                                                               tcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCg                                                                                       >  1:1230592/1‑150 (MQ=255)
                                                               tcGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCg                                                                                       <  1:505666/150‑1 (MQ=255)
                                                                cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGCTAAAAAATATATCAGCCCAGCTCTCTCGGCTGATATTTCTAAACGATCGCACTAAAAAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                      <  2:726726/150‑1 (MQ=255)
                                                                cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                      >  1:1273917/1‑150 (MQ=255)
                                                                   aaaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTa                                                                                   >  2:1569301/1‑150 (MQ=255)
                                                                   aaaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTATCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTa                                                                                   >  2:1449963/1‑150 (MQ=255)
                                                                    aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  >  2:1587397/1‑150 (MQ=255)
                                                                    aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  <  2:1708351/150‑1 (MQ=255)
                                                                    aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  <  1:490102/150‑1 (MQ=255)
                                                                         agcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCtt                                                                             <  2:777281/150‑1 (MQ=255)
                                                                           cagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTc                                                                           <  2:1231966/150‑1 (MQ=255)
                                                                            agATAGTGGTAGCATCGTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCag                                                                         <  2:490602/150‑1 (MQ=255)
                                                                             gATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCag                                                                         >  1:974436/1‑150 (MQ=255)
                                                                             gATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGGCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCag                                                                         <  2:39163/150‑1 (MQ=255)
                                                                               tAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCagag                                                                       <  1:1402359/150‑1 (MQ=255)
                                                                               tAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCagag                                                                       <  2:510669/150‑1 (MQ=255)
                                                                                    gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  1:612538/150‑1 (MQ=255)
                                                                                    gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  2:1324200/150‑1 (MQ=255)
                                                                                    gTAGCATAGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  2:791159/150‑1 (MQ=255)
                                                                                        cATCGTTATAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGATCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACAt                                                              >  2:619977/1‑150 (MQ=255)
                                                                                          tCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTa                                                            <  2:847199/150‑1 (MQ=255)
                                                                                           cGTTTTAACTCATTTGATCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTCGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGCAATACTGCTCTACTTTCGATTAGCCCTTTCACAGGTAAGACATTAg                                                           >  2:1560035/1‑150 (MQ=255)
                                                                                           cGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAg                                                           <  1:1326046/150‑1 (MQ=255)
                                                                                            gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          >  1:1437119/1‑150 (MQ=255)
                                                                                             ttttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAggg                                                         <  1:1423486/150‑1 (MQ=255)
                                                                                                tAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCaa                                                      >  1:26418/1‑150 (MQ=255)
                                                                                                  aCTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAAtt                                                    >  1:1719176/1‑150 (MQ=255)
                                                                                                    tCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTc                                                  <  1:1476503/150‑1 (MQ=255)
                                                                                                     cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCa                                                 <  2:1049939/150‑1 (MQ=255)
                                                                                                      aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                <  1:1647239/150‑1 (MQ=255)
                                                                                                      aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                <  1:852821/150‑1 (MQ=255)
                                                                                                      aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGAACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                >  1:277962/1‑150 (MQ=255)
                                                                                                        ttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                              <  1:1186276/150‑1 (MQ=255)
                                                                                                              aaaaTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGAc                                        <  1:1561344/150‑1 (MQ=255)
                                                                                                                   tCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTgcgc                                   <  2:432752/150‑1 (MQ=255)
                                                                                                                        ttGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCa                              <  1:194930/150‑1 (MQ=255)
                                                                                                                         tGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAt                             <  2:868551/150‑1 (MQ=255)
                                                                                                                            cGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtaca                          <  1:619977/150‑4 (MQ=255)
                                                                                                                             ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacag                         <  2:474164/150‑5 (MQ=255)
                                                                                                                             ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacag                         <  1:466966/150‑5 (MQ=255)
                                                                                                                             ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAg                         >  1:746844/1‑150 (MQ=255)
                                                                                                                              gATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCAtacagt                        <  2:1230592/150‑6 (MQ=255)
                                                                                                                               aTAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTc                       <  2:1653946/150‑1 (MQ=255)
                                                                                                                                 aaaaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTACGACATTACGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtt                     >  2:73187/1‑149 (MQ=255)
                                                                                                                                  aaaTATTTATCAGCCCAGCGCTCTGGGCTGATATTTCTAAATGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtgt                    >  1:524881/1‑150 (MQ=255)
                                                                                                                                    atatTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgg                  >  1:1177925/1‑150 (MQ=255)
                                                                                                                                     tatTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTggg                 >  1:1042514/1‑150 (MQ=255)
                                                                                                                                      atTTATCAGCCCAGCTCTCTGTGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgggg                <  2:371384/150‑1 (MQ=255)
                                                                                                                                       tttATCAGCCCAGCTCTGTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGgt               <  1:158751/150‑1 (MQ=255)
                                                                                                                                         tATCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCTATTTCACTGTATGACTGCGCAGTCATCGAGTGTGTGGGgtgt             >  1:173193/1‑150 (MQ=255)
                                                                                                                                          aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATTACTGCGCAGTAATCGAGTCTGTGGGGTGtt            >  1:1177248/1‑150 (MQ=255)
                                                                                                                                          aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            >  1:199499/1‑150 (MQ=255)
                                                                                                                                          aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            >  1:571316/1‑150 (MQ=255)
                                                                                                                                          aTCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCGGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            >  2:711888/1‑150 (MQ=255)
                                                                                                                                           tCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTCCGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTg           <  2:605548/150‑1 (MQ=255)
                                                                                                                                           tCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCGACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTg           >  1:832679/1‑150 (MQ=255)
                                                                                                                                            cAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGc          >  2:1041306/1‑150 (MQ=255)
                                                                                                                                            cAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCATTCAGAGGTAAGACATTAGGGCAATTACACTGTATGAATGCGCAGTCATCGAGTCTGTGGGGTCTTGc          >  2:597736/1‑150 (MQ=255)
                                                                                                                                             aGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCt         >  1:891620/1‑150 (MQ=255)
                                                                                                                                              gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGGCTGTGGGGTGTTGCTc        >  2:1400780/1‑150 (MQ=255)
                                                                                                                                              gCCCAGCTCTCTGGGCTGAGATTTCTATAAGATCGCACTAAATAGCTACGCGGAATAGGGCTCTACTTTCGGGTAGCCCTTTCCGAGGTAAGACCTTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGGAGCTc        >  2:259230/1‑150 (MQ=255)
                                                                                                                                                   gCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCtttt   <  1:1388878/150‑1 (MQ=255)
                                                                                                                                                   gCTCTATGGGCTGATATTTCTAAAATATCGCACTAAATAGCTACTCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCATTTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCtttt   <  2:304952/150‑1 (MQ=255)
                                                                                                                                                    ctctctGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg  <  2:629350/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757270‑2757567

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