Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,757,419 T→G 7.8% intergenic (+19/+225) 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,4190TG 7.8% 338.8 / 18.8 116intergenic (+19/+225)intA/yfjHCP4‑57 prophage; integrase/CP4‑57 prophage; protein YfjH
Reads supporting (aligned to +/- strand):  ref base T (54/53);  new base G (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 = 7.49e-01

Only 100 of 122 total aligned reads displayed.

GCGCTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTG  >  NC_000913/2757270‑2757567
                                                                                                                                                     |                                                                                                                                                    
gcgcTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGct                                                                                                                                                      >  2:1543896/1‑150 (MQ=255)
   cTTATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctctc                                                                                                                                                   >  2:1416657/1‑150 (MQ=255)
    ttATAACCGCAGTGATTATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGctctct                                                                                                                                                  >  1:629123/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)
         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)
                   ttATCTTGAGTAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   <  1:801045/150‑1 (MQ=255)
                   ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTTTCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   <  2:1427693/150‑1 (MQ=255)
                   ttATCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTa                                                                                                                                   >  1:1148975/1‑150 (MQ=255)
                      tCTTGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTaaaa                                                                                                                                <  1:187969/150‑1 (MQ=255)
                         tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             >  2:210080/1‑150 (MQ=255)
                         tGAGCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGAt                                                                                                                             <  2:1481308/150‑1 (MQ=255)
                            gCAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          >  1:802814/1‑150 (MQ=255)
                            ccAACGTCGTCCGATGATGCAATGGGGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGc                                                                                                                          <  2:339712/149‑1 (MQ=255)
                             cAACGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCa                                                                                                                         >  1:431029/1‑150 (MQ=255)
                              caCGTCGCCCGATCATGCCATCGCCGGCTGACCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        <  2:294064/149‑1 (MQ=255)
                              aaCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCAc                                                                                                                        >  2:780674/1‑150 (MQ=255)
                               aCGTCGTCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACt                                                                                                                       <  2:431029/150‑1 (MQ=255)
                                cgtcgtCCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTa                                                                                                                      <  2:193296/150‑1 (MQ=255)
                                  tcgtcCGATGATGCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTaaa                                                                                                                    <  1:1457621/150‑1 (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:429379/150‑1 (MQ=255)
                                          tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            >  2:1550940/1‑150 (MQ=255)
                                          tgatgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTAcg                                                                                                            <  1:565075/150‑1 (MQ=255)
                                             tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                         >  1:1004881/1‑150 (MQ=255)
                                             tgCAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCgg                                                                                                         >  1:998802/1‑150 (MQ=255)
                                               cAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGaa                                                                                                       >  2:1578818/1‑150 (MQ=255)
                                               cAATGGTGGGCTGATCTCGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCAAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGaa                                                                                                       >  2:527912/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)
                                                                cGTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGt                                                                                      >  2:1379716/1‑150 (MQ=255)
                                                                 gTAAAAGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGtt                                                                                     >  1:1277412/1‑150 (MQ=255)
                                                                   aaaaGCAGCAGATAGTGGTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTAGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTa                                                                                   >  2:593346/1‑150 (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)
                                                                                   ggTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCt                                                                                                               >  1:33773/1‑106 (MQ=255)
                                                                                    gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  1:1054753/150‑1 (MQ=255)
                                                                                    gTAGCATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAg                                                                  <  1:1000314/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)
                                                                                        cATCGTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTc                                                                           >  2:1402765/1‑137 (MQ=255)
                                                                                            gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          >  2:281213/1‑150 (MQ=255)
                                                                                            gTTTTAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAgg                                                          >  1:226025/1‑150 (MQ=255)
                                                                                             ttttAACTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAggg                                                         <  1:104207/150‑1 (MQ=255)
                                                                                                  aataaTTAGAGCAAAATTCGTCTTGTCCGATAAAATATTTATCAGCCCAGCTCGCTGGGCTTAAAATTAAAAAACATCCCACTAAAAAGCAACGCCTAATACTGATCTACTTTCGTTTACCCCTTTCAGAGGAAAGAAATAATGCCAAtt                                                    >  2:179337/5‑150 (MQ=255)
                                                                                                  aCTCATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAAtt                                                    <  1:539459/150‑1 (MQ=255)
                                                                                                    tCATTTGACCAAAATTGGTCTTGGCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTc                                                  <  2:321441/150‑1 (MQ=255)
                                                                                                     cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAATACATTAGGGCAATTTCa                                                 >  2:1463535/1‑150 (MQ=255)
                                                                                                     cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCa                                                 >  2:633023/1‑150 (MQ=255)
                                                                                                     cATTTGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCa                                                 >  1:137093/1‑150 (MQ=255)
                                                                                                       tttGAGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACt                                               <  1:154309/150‑1 (MQ=255)
                                                                                                        ttGGGCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTg                                              <  2:249909/150‑1 (MQ=255)
                                                                                                          gAGCAAAATTCGTCTTGTCGCATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTa                                            <  1:1317656/150‑1 (MQ=255)
                                                                                                            gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:1550940/150‑1 (MQ=255)
                                                                                                            gCAAAATTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          >  2:322501/1‑150 (MQ=255)
                                                                                                              aaaaTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGAc                                        >  2:1075324/1‑150 (MQ=255)
                                                                                                                 aTTCGTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATg                                          <  1:516389/145‑1 (MQ=255)
                                                                                                                     gTCTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCGCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAg                                 <  2:1142125/150‑1 (MQ=255)
                                                                                                                       cTTGTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTc                               >  1:356193/1‑150 (MQ=255)
                                                                                                                          gTCGGATAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATc                            >  1:463435/1‑150 (MQ=255)
                                                                                                                               aTAAAATATTTATCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTc                       >  2:597514/1‑150 (MQ=255)
                                                                                                                                  aaaTATTTATCAGCCAAGCTATCTGGGCTTATATTTATAAAAGATCGAACTAAATAGCTACGCGTAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGAAATTAGGCAAATTTCACTGTATGACTGCGCAGTCATCGAGTCtgt                    >  2:1476467/1‑150 (MQ=255)
                                                                                                                                     gaaatATCAGCCCAGAGAGCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTggg                 >  2:1447643/5‑150 (MQ=255)
                                                                                                                                         tATCAGCCCCGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGgtgt             >  2:1307620/1‑150 (MQ=255)
                                                                                                                                          aTCAGCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGGGTGGTGtt            <  2:1105931/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:802814/150‑1 (MQ=255)
                                                                                                                                              gCCCAGCTCTCTGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTc        >  2:1121481/1‑150 (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)
                                                                                                                                                    ctctctGGGCTGATATTTCTAAAAGATCGCACTAAATAGCTACGCGGAATAGTGCTCTACTTTCGTTTAGCCCTTTCAGAGGTAAGACATTAGGGCAATTTCACTGTATGACTGCGCAGTCATCGAGTCTGTGGGGTGTTGCTCCTTTTg  <  1:996807/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 < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: