Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,120,639 A→G 9.0% L46L (TTA→CTA)  USA300HOU_1994 ← hypothetical bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,120,6390AG 9.0% 242.7 / 13.1 89L46L (TTA→CTA) USA300HOU_1994hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base A (45/36);  new base G (5/3);  total (50/39)
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.47e-01

TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  CP000730/2120500‑2120778
                                                                                                                                           |                                                                                                                                            
tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTACAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:784846/141‑1 (MQ=255)
tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:989805/141‑1 (MQ=255)
tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             >  2:700825/1‑141 (MQ=255)
tctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAc                                                                                                                                             <  1:61745/141‑1 (MQ=255)
 ctcATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAAcc                                                                                                                                            >  1:748148/1‑141 (MQ=255)
   cATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCg                                                                                                                                          >  2:160075/1‑141 (MQ=255)
    aTGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGt                                                                                                                                         >  2:909957/1‑141 (MQ=255)
          ataCCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTg                                                                                                                                   >  1:522359/1‑141 (MQ=255)
           tacctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGa                                                                                                                                  <  1:238647/141‑1 (MQ=255)
           tacctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACCAGTAACCCGTCAATTGa                                                                                                                                  >  1:603800/1‑141 (MQ=255)
            acctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCNCATCAGAAACTANTAACCCGTCAATTGaa                                                                                                                                 <  1:809794/141‑1 (MQ=255)
              ctaccATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAAtt                                                                                                                               <  1:449886/141‑1 (MQ=255)
                   aTTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCtt                                                                                                                          <  2:1030522/141‑1 (MQ=255)
                    tttCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTc                                                                                                                         <  2:890044/141‑1 (MQ=255)
                      tCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCAt                                                                                                                       <  2:385348/141‑1 (MQ=255)
                       cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      <  1:886707/141‑1 (MQ=255)
                       cATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATg                                                                                                                      >  2:723607/1‑141 (MQ=255)
                        aTGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATgg                                                                                                                     >  2:587370/1‑141 (MQ=255)
                           aCTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCt                                                                                                                  <  1:592804/141‑1 (MQ=255)
                           aCTAAGTTAATTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCt                                                                                                                  <  2:1076496/141‑1 (MQ=38)
                               aGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaa                                                                                                              <  1:503076/141‑1 (MQ=255)
                                gttaaTTAGTCTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAaac                                                                                                             >  2:1080817/5‑141 (MQ=38)
                                  tCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCGTAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaacaa                                                                                                           <  1:852966/141‑1 (MQ=255)
                                   cATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGaacaac                                                                                                          >  1:714974/1‑141 (MQ=255)
                                       aGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCGGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAtt                                                                                                      <  1:281820/141‑1 (MQ=255)
                                        gTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCAttt                                                                                                     >  1:872486/1‑141 (MQ=255)
                                             tCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACg                                                                                                >  2:1065026/1‑141 (MQ=255)
                                              ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGAcac                                                                                               >  2:705909/1‑139 (MQ=255)
                                              ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  1:663286/1‑141 (MQ=255)
                                              ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                               >  2:760113/1‑141 (MQ=255)
                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:538008/1‑141 (MQ=255)
                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  2:923004/1‑141 (MQ=255)
                                                  ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                           >  1:644608/1‑141 (MQ=255)
                                                   tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:526314/1‑140 (MQ=255)
                                                   tCATCTGTGTTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACtatat                                                                                          >  1:734959/1‑140 (MQ=255)
                                                   tCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATAc                                                                                          <  2:1003342/141‑1 (MQ=255)
                                                       cTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTc                                                                                      <  1:429292/141‑1 (MQ=255)
                                                         gtgtTATTTTCAATCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCgc                                                                                    <  1:455632/141‑2 (MQ=255)
                                                         gtgtTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGt                                                                                    >  2:635553/1‑141 (MQ=255)
                                                            ttATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTc                                                                                 <  2:629719/141‑1 (MQ=255)
                                                              aTTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCtt                                                                               <  1:299598/141‑1 (MQ=255)
                                                                   cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          <  1:309741/141‑1 (MQ=255)
                                                                   cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                          >  1:372840/1‑141 (MQ=255)
                                                                    aaTCCATTCATAAATACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCGTCAGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAActct                                                                         <  2:82870/141‑1 (MQ=255)
                                                                     aTCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCtt                                                                        <  1:846708/141‑1 (MQ=255)
                                                                      tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:992887/1‑141 (MQ=255)
                                                                      tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                       >  1:475000/1‑141 (MQ=255)
                                                                          ttCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTg                                                                   <  1:641314/141‑1 (MQ=255)
                                                                              tAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGaaaa                                                              >  2:1046637/1‑141 (MQ=255)
                                                                                aaCACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATc                                                            <  2:730749/141‑1 (MQ=255)
                                                                                  cacTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGGTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCaa                                                          >  1:432525/1‑141 (MQ=255)
                                                                                   acTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCaaa                                                         >  2:188678/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  1:940092/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  2:101574/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  1:1050983/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  2:1054135/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  1:219361/1‑141 (MQ=255)
                                                                                      ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                      >  2:944847/1‑141 (MQ=255)
                                                                                        gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCAt                                                    <  1:929619/141‑1 (MQ=255)
                                                                                              aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                               >  1:897473/1‑141 (MQ=255)
                                                                                                atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACatta                                            <  1:460266/141‑1 (MQ=255)
                                                                                                atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACatta                                            <  2:669992/141‑1 (MQ=255)
                                                                                                 tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaa                                           <  1:797050/141‑1 (MQ=255)
                                                                                                   tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaatt                                         <  1:631767/141‑1 (MQ=255)
                                                                                                   tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaatt                                         <  2:46529/141‑1 (MQ=255)
                                                                                                      aaaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAattatt                                      <  1:522513/141‑1 (MQ=255)
                                                                                                       aaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTAttt                                     >  2:1058932/1‑141 (MQ=255)
                                                                                                        aGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTAtttt                                    <  1:576778/141‑1 (MQ=255)
                                                                                                         gCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTAttttt                                   >  1:712845/1‑141 (MQ=255)
                                                                                                            gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGa                                >  2:409349/1‑141 (MQ=255)
                                                                                                               tataGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACAt                             <  1:439432/141‑1 (MQ=255)
                                                                                                                       gTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:617052/141‑1 (MQ=255)
                                                                                                                        ttCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAccc                    >  2:60905/1‑141 (MQ=255)
                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:625438/1‑141 (MQ=255)
                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:883004/1‑141 (MQ=255)
                                                                                                                         tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                   >  2:927714/1‑141 (MQ=255)
                                                                                                                          ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723508/1‑141 (MQ=255)
                                                                                                                          ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                  >  2:723512/1‑141 (MQ=255)
                                                                                                                           tcgtcaGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGtt                 >  1:22892/4‑141 (MQ=255)
                                                                                                                             atcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGtttt               >  1:29322/1‑141 (MQ=255)
                                                                                                                               caGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGt             <  1:34641/141‑1 (MQ=255)
                                                                                                                                aGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCTTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTa            >  1:340122/1‑141 (MQ=255)
                                                                                                                                 gAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTaa           <  2:379579/141‑1 (MQ=255)
                                                                                                                                      tagtagCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACAt      >  1:715115/1‑141 (MQ=255)
                                                                                                                                      tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     >  2:1074934/1‑126 (MQ=255)
                                                                                                                                      tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                     <  1:1074934/126‑1 (MQ=255)
                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  1:1014721/1‑141 (MQ=255)
                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     <  1:923004/141‑1 (MQ=255)
                                                                                                                                       aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt     >  2:993461/1‑141 (MQ=255)
                                                                                                                                        gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg    >  2:176986/1‑141 (MQ=255)
                                                                                                                                          aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt  >  1:26639/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                            
TCTCATGCTTATACCTACCATTTCATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTT  >  CP000730/2120500‑2120778

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