Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,120,697 T→C 6.1% K26K (AAA→AAG USA300HOU_1994 ← hypothetical bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,120,6970TC 6.1% 271.1 / 5.7 95K26K (AAA→AAGUSA300HOU_1994hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base T (52/37);  new base C (3/3);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 6.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

GTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACATCAGC  >  CP000730/2120557‑2120837
                                                                                                                                            |                                                                                                                                             
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:475000/1‑141 (MQ=255)
             tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                                                                                  >  1:992887/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                                                                                                                 >  2:944847/1‑141 (MQ=255)
                             ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                                                                                 >  2:101574/1‑141 (MQ=255)
                             ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                                                                                 >  1:219361/1‑141 (MQ=255)
                             ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                                                                                 >  1:940092/1‑141 (MQ=255)
                             ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                                                                                 >  1:1050983/1‑141 (MQ=255)
                             ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGcc                                                                                                                 >  2:1054135/1‑141 (MQ=255)
                               gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCAt                                                                                                               <  1:929619/141‑1 (MQ=255)
                                     aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                                                                                          >  1:897473/1‑141 (MQ=255)
                                       atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACatta                                                                                                       <  2:669992/141‑1 (MQ=255)
                                       atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACatta                                                                                                       <  1:460266/141‑1 (MQ=255)
                                        tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaa                                                                                                      <  1:797050/141‑1 (MQ=255)
                                          tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaatt                                                                                                    <  2:46529/141‑1 (MQ=255)
                                          tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACattaatt                                                                                                    <  1:631767/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:927714/1‑141 (MQ=255)
                                                                tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                                              >  2:883004/1‑141 (MQ=255)
                                                                tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                                              >  2:625438/1‑141 (MQ=255)
                                                                 ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                                             >  2:723512/1‑141 (MQ=255)
                                                                 ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                                             >  2:723508/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                                                                >  2:993461/1‑141 (MQ=255)
                                                                              aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                                <  1:923004/141‑1 (MQ=255)
                                                                              aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                                                >  1:1014721/1‑141 (MQ=255)
                                                                               gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg                                                               >  2:176986/1‑141 (MQ=255)
                                                                                 aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt                                                             >  1:26639/1‑141 (MQ=255)
                                                                                      gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                                                        >  1:572862/1‑141 (MQ=255)
                                                                                      gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                                                        >  1:789354/1‑141 (MQ=255)
                                                                                        cAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcc                                                      >  2:400888/1‑141 (MQ=255)
                                                                                         aaTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcct                                                     >  2:780143/1‑141 (MQ=255)
                                                                                          aTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATAAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctc                                                    <  1:590673/141‑1 (MQ=255)
                                                                                           ttGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctcc                                                   <  1:635701/141‑1 (MQ=255)
                                                                                            tGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTcctcct                                                  <  2:205492/141‑1 (MQ=255)
                                                                                             gAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCtt                                                 >  2:114196/1‑141 (MQ=255)
                                                                                              aaTTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTg                                                >  1:1084372/1‑141 (MQ=255)
                                                                                                tttCCTCCCTAGTCTAGAACAACCATTTCGACGCAATACGCTCGCTTCTTAACTCTTAATTGAAAAGCAAAGCCATGTACCTTAGTTCTTTTTCGCCATTCGTCACCCGTTTCGTAACGCATTGTTCTCGTCCTCCTTGTc                                              <  2:777534/141‑1 (MQ=255)
                                                                                                 ttCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGGCATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCa                                             <  1:842338/141‑1 (MQ=255)
                                                                                                  tCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCCTCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTATTACATTGTTCTCGTCCTCCTTGTCAt                                            >  1:217886/1‑141 (MQ=255)
                                                                                                  tCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCAt                                            >  1:962909/1‑141 (MQ=255)
                                                                                                   ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTGATACATTGTTCTCGTCCTCCTTGTCAtt                                           >  1:755978/1‑141 (MQ=255)
                                                                                                   ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCAtt                                           >  1:996421/1‑141 (MQ=255)
                                                                                                       cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAACATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTAGTGTAGTACAGTGTTCCTGACGTCTTTGTCATTATCt                                       >  2:271327/1‑141 (MQ=255)
                                                                                                       cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCt                                       <  1:766597/141‑1 (MQ=255)
                                                                                                        aTGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGCAAATCAAAGCCATCTACATTAATTATTTTTCGACATACTTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCta                                      <  2:754503/141‑1 (MQ=255)
                                                                                                         tGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCtat                                     >  1:861688/1‑141 (MQ=255)
                                                                                                               aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATNCATNGTTCTCGTCCTCCTTNTCATTATCTATAGCTAg                               <  1:927714/141‑1 (MQ=255)
                                                                                                               aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg                               <  1:249101/141‑1 (MQ=255)
                                                                                                               aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg                               <  1:700825/141‑1 (MQ=255)
                                                                                                                   cgactaTTTCGACACTATATGCTCGCTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTTTAGTCCTCCTTGTCGTCATCTATACCGAGAAtt                           <  1:705909/136‑1 (MQ=255)
                                                                                                                   caacCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAtt                           >  1:21411/1‑141 (MQ=255)
                                                                                                                   caacCATNTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAtt                           <  2:816252/141‑1 (MQ=255)
                                                                                                                     actaTTTCGACACTATATGCTCGCTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTTTAGTCCTCCTTGTCGTCATCTATACCGAGAAtttt                         >  2:153837/4‑141 (MQ=255)
                                                                                                                      gcATTTCGACGCTGTACGCTCGTTTCTTAACTCTTAATGGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAAttttt                        <  2:638307/140‑1 (MQ=255)
                                                                                                                         tttCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGTCTTACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAAc                     >  2:38862/1‑141 (MQ=255)
                                                                                                                         tttCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAAc                     >  1:66218/1‑141 (MQ=255)
                                                                                                                          ttCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACg                    <  1:189671/141‑1 (MQ=255)
                                                                                                                           tCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGc                   >  2:540546/1‑141 (MQ=255)
                                                                                                                            cGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa                  >  2:366528/1‑141 (MQ=255)
                                                                                                                            cGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCa                  >  2:600363/1‑141 (MQ=255)
                                                                                                                             gACGCTATACGCTCGTTTCTTAACTCTTAATTGATAATCAAAGCCATCTACATTAATTATTTTTCGTCATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAt                 <  2:268482/141‑1 (MQ=255)
                                                                                                                             gACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAt                 <  1:671547/141‑1 (MQ=255)
                                                                                                                                gCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTTTTCATACTTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAtttt              >  2:316756/1‑141 (MQ=255)
                                                                                                                                gCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAtttt              <  1:436533/141‑1 (MQ=255)
                                                                                                                                gCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCAtttt              >  1:949581/1‑141 (MQ=255)
                                                                                                                                  tataCGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGNATTTTGa            >  2:816374/1‑141 (MQ=255)
                                                                                                                                   ataCGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACTTTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGAt           >  1:977693/1‑141 (MQ=255)
                                                                                                                                       gCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACa       >  1:638523/1‑141 (MQ=255)
                                                                                                                                          cGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACATCa    <  1:723508/141‑1 (MQ=255)
                                                                                                                                          cGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACATCa    <  1:723512/141‑1 (MQ=255)
                                                                                                                                            tttCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACATCAGc  >  2:642313/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                             
GTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAGAATTTTTAACGCATTTTGATAACATCAGC  >  CP000730/2120557‑2120837

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