Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 263,351 T→A 5.5% S107C (AGC→TGC)  ykfI ← CP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913263,3510TA 5.5% 494.8 / 11.3 154S107C (AGC→TGC) ykfICP4‑6 prophage; toxin of the YkfI‑YafW toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base T (82/62);  new base A (2/7);  total (85/69)
Fisher's exact test for biased strand distribution p-value = 7.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

Only 100 of 155 total aligned reads displayed.

TCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263202‑263499
                                                                                                                                                     |                                                                                                                                                    
tCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCt                                                                                                                                                      >  2:1022923/1‑150 (MQ=255)
   ggTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTg                                                                                                                                                   >  2:1833626/1‑150 (MQ=255)
   ggTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTg                                                                                                                                                   >  2:2128472/1‑150 (MQ=255)
        cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              >  1:1744263/1‑150 (MQ=255)
        cATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGca                                                                                                                                              >  1:1744262/1‑150 (MQ=255)
          tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaac                                                                                                                                            <  1:2367104/150‑1 (MQ=255)
          tGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaac                                                                                                                                            <  1:2120386/150‑1 (MQ=255)
           gATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaaca                                                                                                                                           >  1:2237860/1‑150 (MQ=255)
            aTTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaa                                                                                                                                          <  2:1229035/150‑1 (MQ=255)
            aTTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGcaacaa                                                                                                                                          >  1:1192621/1‑150 (MQ=255)
             ttgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaac                                                                                                                                         <  2:197263/150‑2 (MQ=255)
              tgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaacg                                                                                                                                        <  1:525724/150‑3 (MQ=255)
              tgtttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGc                                                                                                                                        <  2:615899/150‑1 (MQ=255)
                tttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCa                                                                                                                                      >  2:1439904/1‑150 (MQ=255)
                 ttgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTGCGGCAAAGCCGtaacaacgtag                                                                                                                                     <  1:521863/150‑6 (MQ=255)
                  tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                                                                    <  1:937988/150‑1 (MQ=255)
                  tgtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGt                                                                                                                                    >  2:1046603/1‑150 (MQ=255)
                   gtCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCACAACAAGCCAGtt                                                                                                                                   >  2:2027327/1‑150 (MQ=255)
                        tAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                                                              <  1:1416688/150‑1 (MQ=255)
                        tAAAAAGATAACGCGCCTGCCGGGAAGTAGCAGGCTCATTACGCATTAGGTAAACAAAGGAGCAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctg                                                                                                                              >  2:453106/1‑150 (MQ=255)
                          aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  1:1248280/1‑150 (MQ=255)
                          aaaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCCTTGCCGCAACAAGCCAGTtgcctgcc                                                                                                                            >  2:2142786/1‑150 (MQ=255)
                           aaaGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcg                                                                                                                           <  2:1986640/150‑1 (MQ=255)
                             aGATAACGCGCCTGCCGTGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCcgcg                                                                                                                         <  2:1661359/150‑1 (MQ=255)
                                tAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTc                                                                                                                      <  1:1141877/150‑1 (MQ=255)
                                        cTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATg                                                                                                              <  1:879606/150‑1 (MQ=255)
                                         tGCCGGGTAGTAGCAGGCACATTACGCAATAGGTAAACAAGGGGGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGt                                                                                                             >  2:710041/1‑150 (MQ=255)
                                            cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCACCAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTg                                                                                                          >  2:1347494/1‑150 (MQ=255)
                                            cGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTg                                                                                                          <  1:1606758/150‑1 (MQ=255)
                                              ggTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCa                                                                                                        >  1:440380/1‑150 (MQ=255)
                                                  gtagCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGccc                                                                                                    >  2:1512044/1‑150 (MQ=255)
                                                  gtagCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGccc                                                                                                    >  1:1353823/1‑150 (MQ=255)
                                                   tagCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCg                                                                                                   <  2:1164827/150‑1 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACTGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:1209501/150‑1 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:1567256/150‑1 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  2:383208/150‑1 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  2:1548391/1‑150 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  <  1:1461071/150‑1 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  1:1366531/1‑150 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  2:948997/1‑150 (MQ=255)
                                                    agCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCgg                                                                                                  >  1:1281410/1‑150 (MQ=255)
                                                     gCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGga                                                                                                 >  2:863043/1‑150 (MQ=255)
                                                          cgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                                                            >  2:578381/1‑150 (MQ=255)
                                                          cgcATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATaa                                                                                            >  2:1874634/1‑150 (MQ=255)
                                                            cATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          <  2:491918/150‑1 (MQ=255)
                                                              ttACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGgag                                                                                        <  2:1985040/150‑1 (MQ=255)
                                                                aCGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGAt                                                                                      <  1:1197253/150‑1 (MQ=255)
                                                                  gCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTg                                                                                    <  1:607898/150‑1 (MQ=255)
                                                                   cAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGt                                                                                   >  2:818086/1‑150 (MQ=255)
                                                                   cAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGt                                                                                   >  2:155083/1‑150 (MQ=255)
                                                                        ggTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTg                                                                              >  1:1016706/1‑150 (MQ=255)
                                                                         gTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGc                                                                             >  1:2362402/1‑150 (MQ=255)
                                                                           aaaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCa                                                                           >  2:1314957/1‑150 (MQ=255)
                                                                           aaaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCa                                                                           <  1:2056376/150‑1 (MQ=255)
                                                                           aaaCAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCa                                                                           <  2:1614576/150‑1 (MQ=255)
                                                                               aaGGGAGGAAGTTCAGAAATGTCAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTa                                                                       <  2:1617767/150‑1 (MQ=255)
                                                                               aaGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTa                                                                       <  2:661361/150‑1 (MQ=255)
                                                                                aGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTaa                                                                      >  1:1584955/1‑150 (MQ=255)
                                                                                     ggaAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATccc                                                                 >  2:724043/1‑150 (MQ=255)
                                                                                         gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             <  1:658343/150‑1 (MQ=255)
                                                                                         gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             <  2:2237860/150‑1 (MQ=255)
                                                                                         gTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTcc                                                             >  2:1564106/1‑150 (MQ=255)
                                                                                          ttCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCt                                                            >  1:1173640/1‑150 (MQ=255)
                                                                                           tCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTg                                                           >  1:2157049/1‑150 (MQ=255)
                                                                                            cAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGt                                                          >  1:1022195/1‑150 (MQ=255)
                                                                                            cAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGt                                                          >  2:1355165/1‑150 (MQ=255)
                                                                                             aGACATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGGTCGCAGATTGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTc                                                         <  2:741467/150‑1 (MQ=255)
                                                                                             aGAAATGTAAATCGTGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTc                                                         <  2:749808/150‑1 (MQ=255)
                                                                                              gAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCg                                                        <  1:1518154/150‑1 (MQ=255)
                                                                                              gAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCg                                                        >  1:2327100/1‑150 (MQ=255)
                                                                                                  tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    <  2:1744263/150‑1 (MQ=255)
                                                                                                  tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    <  1:710041/150‑1 (MQ=255)
                                                                                                  tGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGAGCGGAGATAAGTAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATAc                                                    <  2:1744262/150‑1 (MQ=255)
                                                                                                     aaaTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGaa                                                 >  1:90719/1‑150 (MQ=255)
                                                                                                       aTCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAAcc                                               <  1:277305/150‑1 (MQ=255)
                                                                                                                gTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCATCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAt                                      >  2:2183555/1‑150 (MQ=255)
                                                                                                                   gTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTAtttt                                   >  2:1707805/1‑150 (MQ=255)
                                                                                                                     aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          <  2:1975558/93‑1 (MQ=255)
                                                                                                                     aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAgg                                                                                          >  1:1975558/1‑93 (MQ=255)
                                                                                                                     aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAGCTGGTTCGTAtcgaca                                 <  1:2236995/150‑6 (MQ=255)
                                                                                                                     aCGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTc                                 <  2:785788/150‑1 (MQ=255)
                                                                                                                        cAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAc                              >  1:431034/1‑150 (MQ=255)
                                                                                                                        cAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAc                              <  2:1173640/150‑1 (MQ=255)
                                                                                                                        cAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTAc                              >  2:1880702/1‑150 (MQ=255)
                                                                                                                              tCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGaaa                        >  2:870720/1‑150 (MQ=255)
                                                                                                                                aTCGTACTACGTTGGTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAAt                      >  2:1409098/1‑150 (MQ=255)
                                                                                                                                    tactacTTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATATGGAGATTGTTGCTGCCAGCTTAATCCCTTGCTGACGATACGAACCAGATCGTATTTTTCTACCAGAAAATTCAc                  <  2:2026306/150‑1 (MQ=255)
                                                                                                                                     actacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACg                 >  1:150941/1‑150 (MQ=255)
                                                                                                                                      ctggcttgttgCGGCAAAGCCGTAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                >  2:417379/6‑150 (MQ=255)
                                                                                                                                      ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTGGTAGaaaaatacgag                <  1:605230/150‑8 (MQ=255)
                                                                                                                                      ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                >  2:2334551/1‑150 (MQ=255)
                                                                                                                                      ctacGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACgg                <  1:2060308/150‑1 (MQ=255)
                                                                                                                                          gttgttACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATc            >  1:2334819/1‑150 (MQ=255)
                                                                                                                                                 cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCATGGCATCGGCTAgg     <  1:380083/150‑1 (MQ=255)
                                                                                                                                                 cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:994934/150‑1 (MQ=255)
                                                                                                                                                 cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:761184/150‑1 (MQ=255)
                                                                                                                                                 cGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAgg     <  1:1314957/150‑1 (MQ=255)
                                                                                                                                                  ggCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCGCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAggg    <  2:1310435/150‑1 (MQ=255)
                                                                                                                                                   gCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGt   >  1:2437436/1‑150 (MQ=255)
                                                                                                                                                    cTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAGAATTCACGGCATCGGCTAGGGTg  >  2:1100999/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTGGTTGCATGATTGTTTGTCCGTAAAAAGATAACGCGCCTGCCGGGTAGTAGCAGGCGCATTACGCAATAGGTAAACAAGGGAGGAAGTTCAGAAATGTAAATCGGGAAGGTTGTACGCAATGTTCATCGTACTACGTTGTTACGGCTTTGCCGCAACAAGCCAGTTGCCTGCCGCGCTCGCAGAATGTCTGCAGCCCGGAGATAAGGAGATTGTTCCTGCCAGCTAAATCCCTTCCTGTCGATACGAACCAGCTCGTATTTTTCTACCAGAAAATTCACGGCATCGGCTAGGGTG  >  NC_000913/263202‑263499

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