Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,113,768 T→C 4.0% R79G (AGA→GGA)  USA300HOU_1983 ← hypothetical bacteriophage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,113,7680TC 4.0% 306.1 / 2.7 105R79G (AGA→GGA) USA300HOU_1983hypothetical bacteriophage protein
Reads supporting (aligned to +/- strand):  ref base T (54/46);  new base C (2/3);  total (56/49)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 111 total aligned reads displayed.

CGGTTATTCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTC  >  CP000730/2113628‑2113906
                                                                                                                                            |                                                                                                                                          
cgGTTATTCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCgctct                                                                                                                                            >  2:608366/1‑141 (MQ=255)
   tttttCATTTGTGTCATCCTCCATTTGCCCTAAAAATTAGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTGCTTCCTTGTTCgctctgct                                                                                                                                         <  1:545823/138‑1 (MQ=255)
     aTTCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCgctctgctct                                                                                                                                       >  2:26628/1‑141 (MQ=255)
      ttCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCactc                                                                                                                                      >  2:82077/1‑141 (MQ=255)
       tCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCa                                                                                                                                     >  2:518392/1‑141 (MQ=255)
          tttGTGTCATCCTCCATTTGCCCTAAAAACTCGTAGAACTCTTGTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATc                                                                                                                                  <  1:229380/141‑1 (MQ=255)
            tgtgTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTg                                                                                                                                <  1:472682/141‑1 (MQ=255)
            tgtgTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTg                                                                                                                                <  2:623/141‑1 (MQ=255)
             ttgtCACCCTCCATTTGCCCTAAAAATTCGGAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGt                                                                                                                               <  2:109144/140‑1 (MQ=255)
              tgtCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTa                                                                                                                              <  1:289164/141‑1 (MQ=255)
              tgtCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTACTTTGTCCATTCGGCACAATATAACAGTTAAGTTGATTTCAGCTTTTCGTTTTTTAGCTACTTGCTTGTTCGCTCTTCTATCAATCTGtg                                                                                                                              >  1:66990/1‑140 (MQ=255)
                 cATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGtt                                                                                                                           <  2:33683/141‑1 (MQ=255)
                 cATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGtt                                                                                                                           <  1:39426/141‑1 (MQ=255)
                 cATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGtt                                                                                                                           <  2:400410/141‑1 (MQ=255)
                   tcctccATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCg                                                                                                                         >  2:432119/1‑141 (MQ=255)
                           ttGCGCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGAACAATATAACTCTTAAGTTGATTTCAGCTCTTCTATATATAGCTAGTTCCTTGTTCGCTCTGCTCTCAATCTGTAAGTCGCTTAGTCt                                                                                                                 >  1:21176/1‑141 (MQ=255)
                            tGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCCCACTATAACACTTACGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTa                                                                                                                >  1:70971/1‑141 (MQ=255)
                               ccTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTaaaa                                                                                                             >  1:759401/1‑141 (MQ=255)
                                ctcaaaaTTCGTAGAACTCATTTGTTCCGTCTAGTTAATCCCTTCGGCACAATATAACACTTAAGTTGATTTCGGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAAc                                                                                                            <  1:155004/138‑1 (MQ=255)
                                 tAAAAATTCTTAGAACCCATTTTTTCCGTCTAGTTTGTCCTTTCGGCACAATATAACACTTAAGGTGATTTCAGCTCTTCTATATCTCGCTAATTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACg                                                                                                           <  2:343579/141‑1 (MQ=255)
                                   aaaaTTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGt                                                                                                         <  1:513612/141‑1 (MQ=255)
                                      aTTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCCCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTaaa                                                                                                      >  1:326656/1‑141 (MQ=255)
                                        tCGTATACCTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAAt                                                                                                    <  1:178106/141‑1 (MQ=255)
                                        tCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAACACGGTAAAAt                                                                                                    <  2:831617/141‑1 (MQ=255)
                                         cGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAAtt                                                                                                   >  2:419401/1‑141 (MQ=255)
                                         cGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAAtt                                                                                                   <  1:666144/141‑1 (MQ=255)
                                         cGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGCTAAAAtt                                                                                                   >  1:755601/1‑141 (MQ=255)
                                         cGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTACTGGGTGCGCTCTACTATAAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAGtt                                                                                                   >  1:736953/1‑141 (MQ=255)
                                            aGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGt                                                                                                >  1:747617/1‑141 (MQ=255)
                                            aGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGt                                                                                                >  1:151762/1‑141 (MQ=255)
                                            aGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGt                                                                                                >  1:151758/1‑141 (MQ=255)
                                            aGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTAAAGTCAATCGGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGt                                                                                                >  1:781117/1‑141 (MQ=255)
                                               aCTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGtctc                                                                                             >  2:602753/1‑141 (MQ=255)
                                                     ttGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCaa                                                                                       >  1:499492/1‑141 (MQ=255)
                                                     ttGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCCGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCaa                                                                                       >  1:455454/1‑141 (MQ=255)
                                                     ttGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCCGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCaa                                                                                       >  1:41757/1‑141 (MQ=255)
                                                     ttGTTCCGTCTAGTTGGTCCATTCGGCACAATAGAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCaa                                                                                       <  2:670979/141‑1 (MQ=255)
                                                      tGTTCCGTCTAGTTTGGCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCaaa                                                                                      <  1:39940/141‑1 (MQ=255)
                                                       gTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATATGTCGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAAc                                                                                     >  2:556472/1‑141 (MQ=255)
                                                          ccGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAAt                                                                                  >  1:493840/1‑141 (MQ=255)
                                                           cGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAAtt                                                                                 <  1:60790/141‑1 (MQ=255)
                                                            gTCTAGTTTGTCCATTCGGGACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTc                                                                                >  1:625163/1‑141 (MQ=255)
                                                             tCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCa                                                                               >  1:106783/1‑141 (MQ=255)
                                                               tAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCAtt                                                                             <  1:26628/141‑1 (MQ=255)
                                                                     gTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAg                                                                       <  2:240573/141‑1 (MQ=255)
                                                                      tCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGc                                                                      >  1:683105/1‑141 (MQ=255)
                                                                      tCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGc                                                                      >  2:651657/1‑141 (MQ=255)
                                                                           tCGGCACAATATAACACTTAAGTTTATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTgcg                                                                 >  1:535717/1‑141 (MQ=255)
                                                                            cGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTACTACTAAAACGGTAAAATTCGTTTCTTCCATACAATTCATTTTTAAGCGTgcgc                                                                >  1:23412/1‑141 (MQ=255)
                                                                              gCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCTccc                                                              >  1:640355/1‑140 (MQ=255)
                                                                                acaATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAATACGGTAAAAGTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCcaca                                                            <  1:789104/141‑1 (MQ=255)
                                                                                acaATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTCTCCTCCTAGCAATTCATTTTTAAGCGTGCGccccc                                                            >  2:760424/1‑138 (MQ=255)
                                                                                   atatAACACTTAAGTTGCTTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGt                                                         >  1:627114/1‑141 (MQ=255)
                                                                                        acacTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCCGCACTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTcc                                                    <  1:275417/141‑1 (MQ=255)
                                                                                         cacTTAAGTTGTTTTCAGCTCTTCTATATATAGCTACTTCCTTGGTCGCTCTGCTCTCAATCTTTCGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCa                                                   <  2:261277/141‑1 (MQ=255)
                                                                                         cacTTAAGTTGATTTCAGCTCTTCTAAATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCa                                                   <  2:516112/141‑1 (MQ=255)
                                                                                         cacTGAGGTGGATTTCAGCTCTTCTATATATAGCTACTTCCCTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTACACACAAGTCATTTTTAAGCGTGCGCCACATGTTCTCCa                                                   <  2:672606/141‑1 (MQ=255)
                                                                                          cctcacgTTGATTTCAGCTCTTCTATAGATAGCTCCTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAg                                                  <  2:66990/135‑1 (MQ=255)
                                                                                              aaGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCGCATGTTCTCCAGctct                                              <  1:496301/141‑1 (MQ=255)
                                                                                               aGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCCCATGTTCTCCAGCTCtt                                             >  2:249968/1‑141 (MQ=255)
                                                                                                 ttGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAACATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACCTTTCCTCCAGACCTTCg                                           >  2:441433/1‑141 (MQ=255)
                                                                                                 ttGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCg                                           <  2:90107/141‑1 (MQ=255)
                                                                                                 ttGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCg                                           >  2:650491/1‑141 (MQ=255)
                                                                                                 ttGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTAATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTct                                           >  2:395382/1‑140 (MQ=255)
                                                                                                   gATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGtt                                         <  1:466346/141‑1 (MQ=255)
                                                                                                   gATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGtt                                         >  1:246690/1‑141 (MQ=255)
                                                                                                    aTTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTa                                        <  2:501870/141‑1 (MQ=255)
                                                                                                       tCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCATACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGc                                     <  1:556472/141‑1 (MQ=255)
                                                                                                       tCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGc                                     >  2:828612/1‑141 (MQ=255)
                                                                                                       tCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCAACATGTTCTCCAGCTCTTCGTTACGc                                     >  1:540310/1‑141 (MQ=255)
                                                                                                       tCAGCTCTTCTATAGATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAAGTCGTATCTTCCAGACAATTCATTTTTAAGCGTGCGCCCCATGTTCTCCAGCTCTTCGTTCCGc                                     >  2:574392/1‑141 (MQ=255)
                                                                                                        cAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGcc                                    >  2:415947/1‑140 (MQ=255)
                                                                                                             cttctATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCtt                               <  2:85157/141‑1 (MQ=255)
                                                                                                              ttctATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTa                              <  2:750350/141‑1 (MQ=255)
                                                                                                              ttctATATATAGCTACGTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTa                              <  2:502248/141‑1 (MQ=255)
                                                                                                               tctATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCCCCAGCTCTTCGTTACGCCCTCTCaa                             >  2:708811/1‑141 (MQ=255)
                                                                                                                  ttatataGCTGCGTCCTTGTTAGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACtt                          <  2:207781/140‑1 (MQ=255)
                                                                                                                  ttatataGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACtt                          <  2:523067/140‑1 (MQ=255)
                                                                                                                    atataGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAg                        <  1:172465/141‑1 (MQ=255)
                                                                                                                       taGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCta                     <  1:334830/141‑1 (MQ=255)
                                                                                                                       aaGCTACTTCCTTGTGCGATCTGCACTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCta                     <  1:518392/140‑1 (MQ=255)
                                                                                                                        aGCTACTTCCTTGTTCGCTCCGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCtat                    <  2:455454/141‑1 (MQ=255)
                                                                                                                        aGCTACTTCCTTGTTCGCTCCGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCtat                    <  2:41757/141‑1 (MQ=255)
                                                                                                                        aGCTACTTCCTGGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCtat                    <  2:714233/141‑1 (MQ=255)
                                                                                                                         gCTACGTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAGCGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCtata                   <  2:143800/141‑1 (MQ=255)
                                                                                                                            aCTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACACTTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATcc                >  1:45705/1‑141 (MQ=255)
                                                                                                                             cTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATccc               >  2:170500/1‑141 (MQ=255)
                                                                                                                               tCCTTGTTCGCTCTGCTCTCGATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTCTCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTAAGCTCTCTTAACTTAGCTATATCCCCa             >  1:511060/1‑141 (MQ=255)
                                                                                                                               tCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCa             <  2:765549/141‑1 (MQ=255)
                                                                                                                                 cTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAAt           >  1:777843/1‑141 (MQ=255)
                                                                                                                                 cTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACCCTCTCTTAACTTAGCTATATCCCCaaa           >  1:514841/1‑140 (MQ=255)
                                                                                                                                  ttGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCCAACAATTCATGTTTCAGAGTGCGAAAAATTTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATa          >  1:649492/1‑141 (MQ=255)
                                                                                                                                  ttGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTCTATCCCCAATa          >  2:144408/1‑141 (MQ=255)
                                                                                                                                  ttGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATa          >  1:521084/1‑141 (MQ=255)
                                                                                                                                  ttGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATa          >  1:454656/1‑141 (MQ=255)
                                                                                                                                  ttGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATa          >  1:98816/1‑141 (MQ=255)
                                                                                                                                    gTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAg        <  2:535717/141‑1 (MQ=255)
                                                                                                                                     ttCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGc       >  1:836703/1‑141 (MQ=255)
                                                                                                                                        gctctgctctCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCg    <  1:154711/141‑1 (MQ=255)
                                                                                                                                          tctgctctCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTc  >  1:488599/1‑141 (MQ=255)
                                                                                                                                            |                                                                                                                                          
CGGTTATTCATTTGTGTCATCCTCCATTTGCCCTAAAAATTCGTAGAACTCATTTGTTCCGTCTAGTTTGTCCATTCGGCACAATATAACACTTAAGTTGATTTCAGCTCTTCTATATATAGCTACTTCCTTGTTCGCTCTGCTCTCAATCTGTAGTTCGCTAAGTCTAAAACGGTAAAATTCGTATCTTCCAAACAATTCATTTTTAAGCGTGCGCCACATGTTCTCCAGCTCTTCGTTACGCTCTCTTAACTTAGCTATATCCCCAATAAGCTCGTC  >  CP000730/2113628‑2113906

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