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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,120,6690GA 8.4% 276.8 / 9.6 93L36L (CTA→TTA) USA300HOU_1994hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base G (50/35);  new base A (6/2);  total (56/37)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

AGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAG  >  CP000730/2120531‑2120808
                                                                                                                                          |                                                                                                                                            
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:923004/1‑141 (MQ=255)
                   ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                         >  1:644608/1‑141 (MQ=255)
                   ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                         >  2:538008/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:372840/1‑141 (MQ=255)
                                    cAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAActc                                                                                                        <  1:309741/141‑1 (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)
                                               tAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑aaaa                                                                                            >  2:1046637/1‑141 (MQ=255)
                                                 aaCACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATc                                                                                          <  2:730749/141‑1 (MQ=255)
                                                   cacTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGGTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaa                                                                                        >  1:432525/1‑141 (MQ=255)
                                                    acTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCaaa                                                                                       >  2:188678/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  2:944847/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  1:1050983/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  2:1054135/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  2:101574/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  1:219361/1‑141 (MQ=255)
                                                       ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                                    >  1:940092/1‑141 (MQ=255)
                                                         gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                                                  <  1:929619/141‑1 (MQ=255)
                                                               aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                                                             >  1:897473/1‑141 (MQ=255)
                                                                 atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                          <  2:669992/141‑1 (MQ=255)
                                                                 atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                          <  1:460266/141‑1 (MQ=255)
                                                                  tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaa                                                                         <  1:797050/141‑1 (MQ=255)
                                                                    tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                                                       <  1:631767/141‑1 (MQ=255)
                                                                    tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                                                       <  2:46529/141‑1 (MQ=255)
                                                                       aaaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAattatt                                                                    <  1:522513/141‑1 (MQ=255)
                                                                        aaGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttt                                                                   >  2:1058932/1‑141 (MQ=255)
                                                                         aGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAtttt                                                                  <  1:576778/141‑1 (MQ=255)
                                                                          gCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTAttttt                                                                 >  1:712845/1‑141 (MQ=255)
                                                                             gtgtATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGa                                                              >  2:409349/1‑141 (MQ=255)
                                                                                tataGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACAt                                                           <  1:439432/141‑1 (MQ=255)
                                                                                        gTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                   <  1:617052/141‑1 (MQ=255)
                                                                                         ttCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAccc                                                  >  2:60905/1‑141 (MQ=255)
                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                 >  2:883004/1‑141 (MQ=255)
                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                 >  2:927714/1‑141 (MQ=255)
                                                                                          tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                                 >  2:625438/1‑141 (MQ=255)
                                                                                           ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                >  2:723508/1‑141 (MQ=255)
                                                                                           ctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGt                                                >  2:723512/1‑141 (MQ=255)
                                                                                            tcgtcaGAAACTAGTAGCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGtt                                               >  1:22892/4‑141 (MQ=255)
                                                                                              atcaGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGtttt                                             >  1:29322/1‑141 (MQ=255)
                                                                                                caGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGt                                           <  1:34641/141‑1 (MQ=255)
                                                                                                 aGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCTTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTa                                          >  1:340122/1‑141 (MQ=255)
                                                                                                  gAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTaa                                         <  2:379579/141‑1 (MQ=255)
                                                                                                       tagtagCCCGTCAATTGAATTTCCTTCATGATCTAAAACGACTATTTCGACACTATATGCTCGCTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAAATATTTTTCGACATACGTCACCCGTTTTGTAATACAt                                    >  1:715115/1‑141 (MQ=255)
                                                                                                       tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                   >  2:1074934/1‑126 (MQ=255)
                                                                                                       tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                                   <  1:1074934/126‑1 (MQ=255)
                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                   >  2:993461/1‑141 (MQ=255)
                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                   <  1:923004/141‑1 (MQ=255)
                                                                                                        aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                                   >  1:1014721/1‑141 (MQ=255)
                                                                                                         gTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTg                                  >  2:176986/1‑141 (MQ=255)
                                                                                                           aaCCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGtt                                >  1:26639/1‑141 (MQ=255)
                                                                                                                gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                           >  1:572862/1‑141 (MQ=255)
                                                                                                                gTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGt                           >  1:789354/1‑141 (MQ=255)
                                                                                                                  cAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcc                         >  2:400888/1‑141 (MQ=255)
                                                                                                                   aaTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcct                        >  2:780143/1‑141 (MQ=255)
                                                                                                                    aTTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATAAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctc                       <  1:590673/141‑1 (MQ=255)
                                                                                                                     ttGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGtcctcc                      <  1:635701/141‑1 (MQ=255)
                                                                                                                      tGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTcctcct                     <  2:205492/141‑1 (MQ=255)
                                                                                                                       gAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCtt                    >  2:114196/1‑141 (MQ=255)
                                                                                                                        aaTTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTg                   >  1:1084372/1‑141 (MQ=255)
                                                                                                                          tttCCTCCCTAGTCTAGAACAACCATTTCGACGCAATACGCTCGCTTCTTAACTCTTAATTG‑AAAAGCAAAGCCATGTACCTTAGTTCTTTTTCGCCATTCGTCACCCGTTTCGTAACGCATTGTTCTCGTCCTCCTTGTc                 <  2:777534/141‑1 (MQ=255)
                                                                                                                           ttCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGGCATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCa                <  1:842338/141‑1 (MQ=255)
                                                                                                                            tCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCCTCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTATTACATTGTTCTCGTCCTCCTTGTCAt               >  1:217886/1‑141 (MQ=255)
                                                                                                                            tCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCAt               >  1:962909/1‑141 (MQ=255)
                                                                                                                             ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTGATACATTGTTCTCGTCCTCCTTGTCAtt              >  1:755978/1‑141 (MQ=255)
                                                                                                                             ccTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCAtt              >  1:996421/1‑141 (MQ=255)
                                                                                                                                 cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AACATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTAGTGTAGTACAGTGTTCCTGACGTCTTTGTCATTATCt          >  2:271327/1‑141 (MQ=255)
                                                                                                                                 cATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCt          <  1:766597/141‑1 (MQ=255)
                                                                                                                                  aTGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑CAAATCAAAGCCATCTACATTAATTATTTTTCGACATACTTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCta         <  2:754503/141‑1 (MQ=255)
                                                                                                                                   tGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCtat        >  1:861688/1‑141 (MQ=255)
                                                                                                                                         aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATNCATNGTTCTCGTCCTCCTTNTCATTATCTATAGCTAg  <  1:927714/141‑1 (MQ=255)
                                                                                                                                         aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg  <  1:249101/141‑1 (MQ=255)
                                                                                                                                         aGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAg  <  1:700825/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                            
AGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTATAGCTAG  >  CP000730/2120531‑2120808

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