Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,418 C→A 7.7% intergenic (+18/+226) intA → / ← yfjH CP4‑57 prophage; integrase/CP4‑57 prophage; protein YfjH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,757,4180CA 7.7% 383.2 / 13.4 116intergenic (+18/+226)intA/yfjHCP4‑57 prophage; integrase/CP4‑57 prophage; protein YfjH
Reads supporting (aligned to +/- strand):  ref base C (54/53);  new base A (4/5);  total (58/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 123 total aligned reads displayed.

CGCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757269‑2757567
                                                                                                                                                     |                                                                                                                                                     
cgcgCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGc                                                                                                                                                       <  2:333025/150‑1 (MQ=255)
    cTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctctc                                                                                                                                                   >  2:1416657/1‑150 (MQ=255)
       ataACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTgg                                                                                                                                                >  1:1286766/1‑150 (MQ=255)
       ataACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGAAGCAATGGTGGGCTGATCTCGTAAAAGCAGAAGATAGTGGTAGCATCGTTTTAACTCATTTTAGAAAAATTCGTCTTGTCGGATAAAATATTTATCAGCACAGCTCTCTgg                                                                                                                                                >  2:535723/1‑150 (MQ=255)
        taACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTggg                                                                                                                                               >  2:353385/1‑150 (MQ=255)
        taACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATAATTATCAGCCCAGCTCTCTggg                                                                                                                                               >  2:16345/1‑150 (MQ=255)
          aCCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCt                                                                                                                                             <  1:859279/150‑1 (MQ=255)
            cGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGa                                                                                                                                           <  1:700464/150‑1 (MQ=255)
             gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGGGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGat                                                                                                                                          >  1:249909/1‑150 (MQ=255)
             gCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGat                                                                                                                                          <  1:356501/150‑1 (MQ=255)
              cAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGata                                                                                                                                         <  1:1614458/150‑1 (MQ=255)
              cAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGata                                                                                                                                         <  1:585924/150‑1 (MQ=255)
               aGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGAGAGCTGGGCTGataa                                                                                                                                        <  2:436701/150‑2 (MQ=255)
                   aTTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCTTCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCt                                                                                                                                    <  2:131162/150‑1 (MQ=255)
                    ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTTTCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   <  2:1427693/150‑1 (MQ=255)
                    ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   >  1:1148975/1‑150 (MQ=255)
                          tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             <  2:1481308/150‑1 (MQ=255)
                          tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             >  2:210080/1‑150 (MQ=255)
                             gCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          >  1:802814/1‑150 (MQ=255)
                             ccAACGTCGTCCGATGATGCAATGGGGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          <  2:339712/149‑1 (MQ=255)
                               caCGTCGCCCGATCATGCCATCGCCGGCTGACCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        <  2:294064/149‑1 (MQ=255)
                               aaCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        >  2:780674/1‑150 (MQ=255)
                                aCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACt                                                                                                                       <  2:431029/150‑1 (MQ=255)
                                   tcgtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTaaa                                                                                                                    <  1:1457621/150‑1 (MQ=255)
                                    cgtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGAAAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAAt                                                                                                                   >  1:445045/1‑150 (MQ=255)
                                     gtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATa                                                                                                                  <  2:1098319/150‑1 (MQ=255)
                                      tcCGATGATGGAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAg                                                                                                                 <  2:1509636/150‑1 (MQ=255)
                                      tcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCGAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGa                                                                                                                                           <  2:1108879/124‑1 (MQ=255)
                                      tcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCGAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGa                                                                                                                                           >  1:1108879/1‑124 (MQ=255)
                                           tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  2:431409/150‑1 (MQ=255)
                                           tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  1:565075/150‑1 (MQ=255)
                                           tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  1:429379/150‑1 (MQ=255)
                                           tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            >  2:1550940/1‑150 (MQ=255)
                                             atgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcgcg                                                                                                          <  1:353385/150‑1 (MQ=255)
                                              tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                         >  1:998802/1‑150 (MQ=255)
                                              tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                         >  1:1004881/1‑150 (MQ=255)
                                                cAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGaa                                                                                                       >  2:1578818/1‑150 (MQ=255)
                                                  aTGGTGTTCTTATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATa                                                                                                     <  1:492672/150‑1 (MQ=255)
                                                   tggtggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAg                                                                                                    <  1:1416657/150‑1 (MQ=255)
                                                      tggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGc                                                                                                 <  2:289249/150‑1 (MQ=255)
                                                      tggGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGc                                                                                                 >  1:215041/1‑150 (MQ=255)
                                                                 cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                      >  2:1379716/1‑150 (MQ=255)
                                                                  gTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGtt                                                                                     >  1:1277412/1‑150 (MQ=255)
                                                                    aaaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTAGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTa                                                                                   >  2:593346/1‑150 (MQ=255)
                                                                     aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  <  1:674922/150‑1 (MQ=255)
                                                                     aaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAg                                                                                  >  2:273448/1‑150 (MQ=255)
                                                                       agcagcAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGcc                                                                                >  1:1534124/1‑150 (MQ=255)
                                                                          agcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCGCtt                                                                             <  1:593346/150‑1 (MQ=255)
                                                                           gcagATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCttt                                                                            >  1:154248/1‑150 (MQ=255)
                                                                              gATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCag                                                                         <  1:252076/150‑1 (MQ=255)
                                                                                tAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCagag                                                                       <  2:498714/150‑1 (MQ=255)
                                                                                tAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCagag                                                                       <  2:147787/150‑1 (MQ=255)
                                                                                 aGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAgg                                                                      >  2:132951/1‑150 (MQ=255)
                                                                                  gTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGt                                                                     >  1:1487765/1‑150 (MQ=255)
                                                                                    ggTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCt                                                                                                               <  2:33773/106‑1 (MQ=255)
                                                                                     gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  1:1000314/150‑1 (MQ=255)
                                                                                     gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  1:1054753/150‑1 (MQ=255)
                                                                                      tAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGa                                                                 <  1:1233198/150‑1 (MQ=255)
                                                                                       aGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGAc                                                                >  1:250619/1‑150 (MQ=255)
                                                                                         cATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCTCTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACAt                                                              >  2:218972/1‑150 (MQ=255)
                                                                                         cATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTc                                                                           <  1:1402765/137‑1 (MQ=255)
                                                                                             gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          >  1:226025/1‑150 (MQ=255)
                                                                                              ttttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAggg                                                         <  1:104207/150‑1 (MQ=255)
                                                                                                  aaCTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGATCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAAt                                                     >  2:266009/1‑150 (MQ=255)
                                                                                                   aataaTTAGAGCAAAATTCGTCTTGTCCGATAAAATATTTATCAGCCCAGCTCGCTGGGCTTAAAATTAAAAAACATCCCACTAAAAAGCAACGCCTAATACTGATCTACTTTCGTTTACCCCTTTCAGAGGAAAGAAATAATGCCAAtt                                                    >  2:179337/5‑150 (MQ=255)
                                                                                                     tCATTTGACCAAAATTGGTCTTGGCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTc                                                  <  2:321441/150‑1 (MQ=255)
                                                                                                      cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAATACATTAGGGCAATTTCa                                                 >  2:1463535/1‑150 (MQ=255)
                                                                                                      cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCa                                                 >  1:137093/1‑150 (MQ=255)
                                                                                                      cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCa                                                 >  2:633023/1‑150 (MQ=255)
                                                                                                       aTTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCAc                                                <  2:1486918/150‑1 (MQ=255)
                                                                                                         ttGGGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                              <  2:249909/150‑1 (MQ=255)
                                                                                                          tGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGt                                             <  1:1543896/150‑1 (MQ=255)
                                                                                                             gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:1550940/150‑1 (MQ=255)
                                                                                                             gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          >  2:1118159/1‑150 (MQ=255)
                                                                                                             gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:580961/150‑1 (MQ=255)
                                                                                                              cAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTATCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGa                                         >  2:1436529/1‑150 (MQ=255)
                                                                                                               aaaaTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGAc                                        >  2:1075324/1‑150 (MQ=255)
                                                                                                                  aTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          >  2:516389/1‑145 (MQ=255)
                                                                                                                  aTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:516389/145‑1 (MQ=255)
                                                                                                                      gTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCGCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAg                                 <  2:1142125/150‑1 (MQ=255)
                                                                                                                        cTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTc                               >  1:356193/1‑150 (MQ=255)
                                                                                                                           gTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTTCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATc                            <  1:1360230/150‑1 (MQ=255)
                                                                                                                           gTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATc                            >  1:463435/1‑150 (MQ=255)
                                                                                                                              ggATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAg                         >  2:1366081/1‑150 (MQ=255)
                                                                                                                                aTAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTc                       >  2:597514/1‑150 (MQ=255)
                                                                                                                                  aaaaTATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCtg                     <  2:1012391/150‑1 (MQ=255)
                                                                                                                                      tatTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTggg                 >  2:304796/1‑150 (MQ=255)
                                                                                                                                      gaaatATCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTggg                 >  2:1447643/5‑150 (MQ=255)
                                                                                                                                       atTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTgggg                <  1:805018/150‑1 (MQ=255)
                                                                                                                                           aTCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGtt            >  1:581023/1‑150 (MQ=255)
                                                                                                                                            tCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGACATCGAGTCTGTGGGGTGTTg           >  1:990537/1‑150 (MQ=255)
                                                                                                                                             cAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGc          <  1:1578818/150‑1 (MQ=255)
                                                                                                                                             cAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTCCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGc          >  2:573895/1‑150 (MQ=255)
                                                                                                                                               gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTc        >  2:1121481/1‑150 (MQ=255)
                                                                                                                                               gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTc        <  2:802814/150‑1 (MQ=255)
                                                                                                                                                 ccAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCt      <  1:1494023/150‑1 (MQ=255)
                                                                                                                                                    gCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCtttt   <  1:79931/150‑1 (MQ=255)
                                                                                                                                                    gCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCtttt   <  2:654833/150‑1 (MQ=255)
                                                                                                                                                     ctctctGGGCTGATATTTCTAAAATATCTCACTAAATAGCTACGCGGAATAGTGATCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGAAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg  <  2:1081404/150‑1 (MQ=255)
                                                                                                                                                     ctctctGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg  <  1:770175/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757269‑2757567

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