Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,762,221 G→T 5.9% I440I (ATC→ATA) ‡ abpB ← CP4‑57 prophage; putative helicase YfjK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,762,2210GT 5.9% 284.7 / 5.6 85I440I (ATC→ATA) ‡abpBCP4‑57 prophage; putative helicase YfjK
Reads supporting (aligned to +/- strand):  ref base G (42/38);  new base T (3/2);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01

CAAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAGG  >  NC_000913/2762079‑2762369
                                                                                                                                              |                                                                                                                                                    
caAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGc                                                                                                                                               >  2:1170532/1‑150 (MQ=255)
caAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGc                                                                                                                                               <  1:384266/150‑1 (MQ=255)
  aaGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCcg                                                                                                                                             <  2:1005132/150‑1 (MQ=255)
    gCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctgg                                                                                                                                           <  2:953561/150‑4 (MQ=255)
      gATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAATCTTTTGcctggag                                                                                                                                         <  1:1366402/150‑6 (MQ=255)
      gATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCcgacga                                                                                                                                         <  2:1112570/150‑1 (MQ=255)
       aTCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGaa                                                                                                                                        >  1:1024128/1‑150 (MQ=255)
        tCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAg                                                                                                                                       <  2:448690/150‑1 (MQ=255)
        tCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAg                                                                                                                                       <  2:1024128/150‑1 (MQ=255)
          cTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTg                                                                                                                                     <  2:75182/150‑1 (MQ=255)
            gtgAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGaa                                                                                                                                   >  2:510988/1‑150 (MQ=255)
                aaaCTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGcc                                                                                                                               >  1:102111/1‑150 (MQ=255)
                    tCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATg                                                                                                                           <  1:881489/150‑1 (MQ=255)
                    tCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATg                                                                                                                           <  1:876299/150‑1 (MQ=255)
                     cGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATgc                                                                                                                          <  1:80020/150‑1 (MQ=255)
                       gACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATgcgc                                                                                                                        >  1:13535/1‑150 (MQ=255)
                        aCAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGcc                                                                                                                       <  1:1306290/150‑1 (MQ=255)
                           aGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCgg                                                                                                                    >  2:102318/1‑150 (MQ=255)
                              ggTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGgcgc                                                                                                                 >  2:835383/1‑150 (MQ=255)
                              ggTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGgcgc                                                                                                                 <  2:868787/150‑1 (MQ=255)
                              ggTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGgcgc                                                                                                                 >  2:1133238/1‑150 (MQ=255)
                               gTCTGCTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGcgcg                                                                                                                >  2:80344/1‑150 (MQ=255)
                                 cTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGAc                                                                                                              <  1:64393/150‑1 (MQ=255)
                                         aCGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATg                                                                                                      <  2:90112/150‑1 (MQ=255)
                                              ggCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTccc                                                                                                 >  2:994518/1‑150 (MQ=255)
                                                 tAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                              >  1:998405/1‑150 (MQ=255)
                                                 tAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                              <  2:411638/150‑1 (MQ=255)
                                                      ggTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGa                                                                                         >  1:881621/1‑150 (MQ=255)
                                                      ggTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGa                                                                                         >  2:1095637/1‑150 (MQ=255)
                                                      ggTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGa                                                                                         >  2:1349397/1‑150 (MQ=255)
                                                         tGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAAtt                                                                                      >  1:938725/1‑150 (MQ=255)
                                                         tGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAAtt                                                                                      <  2:109124/150‑1 (MQ=255)
                                                            ttGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGt                                                                                   >  1:688295/1‑150 (MQ=255)
                                                             tGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTc                                                                                  <  2:260574/150‑1 (MQ=255)
                                                              ggggTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCa                                                                                 <  2:998405/150‑1 (MQ=255)
                                                                     cAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAAtgctgc                                                                          >  1:557241/1‑150 (MQ=255)
                                                                     cAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAAtgctgc                                                                          >  2:946786/1‑150 (MQ=255)
                                                                      aaTGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTCTCAATCCTGCg                                                                         >  2:1013132/1‑150 (MQ=255)
                                                                       aTGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGa                                                                        <  1:1349397/150‑1 (MQ=255)
                                                                         gCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGt                                                                      >  1:1193852/1‑150 (MQ=255)
                                                                         gCCCTGAATACCGAGCGGAATATCCACCTCTCGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGt                                                                      >  1:1154201/1‑150 (MQ=255)
                                                                          cccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                              <  2:1243306/125‑1 (MQ=255)
                                                                          cccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATg                                                                                              >  1:1243306/1‑125 (MQ=255)
                                                                           ccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTc                                                                    >  1:999425/1‑150 (MQ=255)
                                                                           ccTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTc                                                                    <  2:501964/150‑1 (MQ=255)
                                                                            cTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTcc                                                                   >  1:566955/1‑150 (MQ=255)
                                                                                   ccGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCgg                                                            >  1:181346/1‑150 (MQ=255)
                                                                                   ccGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCgg                                                            >  2:1095425/1‑150 (MQ=255)
                                                                                     gAGCGGAATATCCACCTGTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGtt                                                          <  2:881621/150‑1 (MQ=255)
                                                                                      cgcgcAATATCCACCGCTTGATTGAGATGGTCTTCCGGGATTGCCCCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTg                                                         <  2:947184/149‑1 (MQ=255)
                                                                                       gCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGt                                                        <  2:1193852/150‑1 (MQ=255)
                                                                                       gCGGAATATCCACCTCTTGATTCCGCGTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGCGTGAATTTGTCAATGCTGCGAGCTCCATCGCGGTTGt                                                        <  2:1270805/150‑1 (MQ=255)
                                                                                        cGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTc                                                       <  1:1133238/150‑1 (MQ=255)
                                                                                           aaTATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTa                                                    <  1:80344/150‑1 (MQ=255)
                                                                                           aaTATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTa                                                    <  1:946786/150‑1 (MQ=255)
                                                                                           aaTATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTa                                                    <  1:510988/150‑1 (MQ=255)
                                                                                               tCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCTAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCTATGCTGGGAGTTCCATCGCGGTTGTCGTAGATc                                                >  2:920153/1‑150 (MQ=255)
                                                                                               tCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCGAAAGATTTTGCCTACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATc                                                >  2:220228/1‑150 (MQ=255)
                                                                                               tCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCg                                                                         >  1:1031684/1‑125 (MQ=255)
                                                                                               tCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCg                                                                         <  2:1031684/125‑1 (MQ=255)
                                                                                               tCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATc                                                <  2:13535/150‑1 (MQ=255)
                                                                                               tCCACATCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATc                                                >  2:694738/1‑150 (MQ=255)
                                                                                                ccACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATca                                               >  2:452397/1‑150 (MQ=255)
                                                                                                 cACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATcac                                              >  2:396101/1‑150 (MQ=255)
                                                                                                   ccTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACaa                                            >  2:448510/1‑150 (MQ=255)
                                                                                                        tgattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTc                                       >  1:1041990/1‑150 (MQ=255)
                                                                                                        tgattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTc                                       >  1:1037052/1‑150 (MQ=255)
                                                                                                        tgattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTc                                       <  1:694738/150‑1 (MQ=255)
                                                                                                         gattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCt                                      <  1:369202/150‑1 (MQ=255)
                                                                                                         gattgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGACCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCt                                      >  1:719168/1‑150 (MQ=255)
                                                                                                          attgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCtt                                     <  1:835383/150‑1 (MQ=255)
                                                                                                            tgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAg                                   <  2:833190/150‑1 (MQ=255)
                                                                                                            tgaGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAg                                   >  1:374204/1‑150 (MQ=255)
                                                                                                                   gTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGCCGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAgt                            >  2:801409/1‑150 (MQ=255)
                                                                                                                         cGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGAc                      >  2:351706/1‑150 (MQ=255)
                                                                                                                             aTTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCt                  >  2:1111503/1‑150 (MQ=255)
                                                                                                                                    gtcggCAAAATCTTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGa           >  1:1268061/4‑150 (MQ=255)
                                                                                                                                     tcggCAAAATCTTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGAt          >  1:19579/3‑150 (MQ=255)
                                                                                                                                     cAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGAt          <  2:938725/150‑1 (MQ=255)
                                                                                                                                        gCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGt       >  1:206492/1‑150 (MQ=255)
                                                                                                                                         cAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTg      <  1:452397/150‑1 (MQ=255)
                                                                                                                                          aaaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTgg     >  1:1263827/1‑150 (MQ=255)
                                                                                                                                           aaaGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGa    >  1:543834/1‑150 (MQ=255)
                                                                                                                                             aGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAgg  <  1:994518/150‑1 (MQ=255)
                                                                                                                                             aGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAgg  <  1:396101/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CAAAGCGATCCTGTGAAAACTCGGACAAGTGGTCTGGTTGAACGCTGGCTAACAGGTTGATTGGGGTATCAATGCCCTGAATACCGAGCGGAATATCCACCTCTTGATTGAGATTGTCTTCCGGGATTTCCTCCAGGCAAAAGATTTTGCCGACGAAGTGAACGCCCATGCGCCCGGCGCGACCTTTGATGTTCCCATGGGTGAATTTGTCAATGCTGCGAGTTCCATCGCGGTTGTCGTAGATCACAACGTTCTTAGCGATAGTGTTGACTCCTTCAATGATCGTGGAGG  >  NC_000913/2762079‑2762369

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: