Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,120,664 G→A 9.8% D37D (GAC→GAT)  USA300HOU_1994 ← hypothetical bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,120,6640GA 9.8% 264.4 / 14.1 92D37D (GAC→GAT) USA300HOU_1994hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base G (50/33);  new base A (6/3);  total (56/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTAT  >  CP000730/2120524‑2120802
                                                                                                                                            |                                                                                                                                           
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                                                                                                                       >  2:760113/1‑141 (MQ=255)
                      ccTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGc                                                                                                                       >  1:663286/1‑141 (MQ=255)
                          ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                   >  1:644608/1‑141 (MQ=255)
                          ttCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCtata                                                                                                                   >  2:923004/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:992887/1‑141 (MQ=255)
                                              tCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTa                                                                                               >  1:475000/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                                                                              >  1:940092/1‑141 (MQ=255)
                                                              ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                              >  2:101574/1‑141 (MQ=255)
                                                              ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                              >  2:1054135/1‑141 (MQ=255)
                                                              ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                              >  1:1050983/1‑141 (MQ=255)
                                                              ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                              >  2:944847/1‑141 (MQ=255)
                                                              ttGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGcc                                                                              >  1:219361/1‑141 (MQ=255)
                                                                gTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCAt                                                                            <  1:929619/141‑1 (MQ=255)
                                                                      aaaTATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTGAAAAATCAAAGCCATCTACa                                                                       >  1:897473/1‑141 (MQ=255)
                                                                        atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                    <  1:460266/141‑1 (MQ=255)
                                                                        atatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACatta                                                                    <  2:669992/141‑1 (MQ=255)
                                                                         tatCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaa                                                                   <  1:797050/141‑1 (MQ=255)
                                                                           tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                                                 <  2:46529/141‑1 (MQ=255)
                                                                           tCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACattaatt                                                                 <  1:631767/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:927714/1‑141 (MQ=255)
                                                                                                 tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                           >  2:625438/1‑141 (MQ=255)
                                                                                                 tctcATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCg                                           >  2:883004/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                                             <  1:1074934/126‑1 (MQ=255)
                                                                                                              tAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCAcc                                             >  2:1074934/1‑126 (MQ=255)
                                                                                                               aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                             >  1:1014721/1‑141 (MQ=255)
                                                                                                               aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                             >  2:993461/1‑141 (MQ=255)
                                                                                                               aGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACAtt                             <  1:923004/141‑1 (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)
                                                                                                                                            |                                                                                                                                           
ATGACTAAGTTCATTAGTTTGTCCTGTTCATCTGTGTTATTTTCAATCCATTCATAAACACTTTGTTTCAAAATATCTAAAGCTGTGTATAGATCGTTCTCATCAGAAACTAGTAACCCGTCAATTGAATTTCCTTCATGGTCTAGAACAACCATTTCGACGCTATACGCTCGTTTCTTAACTCTTAATTG‑AAAATCAAAGCCATCTACATTAATTATTTTTCGACATACGTCACCCGTTTTGTAATACATTGTTCTCGTCCTCCTTGTCATTATCTAT  >  CP000730/2120524‑2120802

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: ‑