Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,107,541 A→T 4.3% V256D (GTT→GAT)  USA300HOU_1970 ← possible bacteriophage portal protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,107,5410AT 4.3% 277.0 / 3.3 94V256D (GTT→GAT) USA300HOU_1970possible bacteriophage portal protein
Reads supporting (aligned to +/- strand):  ref base A (43/47);  new base T (2/2);  total (45/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

TATGAGTTTCGCGTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTAAAAGTAT  >  CP000730/2107402‑2107679
                                                                                                                                           |                                                                                                                                          
tATGAGTTTCGCGTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAAc                                                                                                                                           <  1:436817/141‑1 (MQ=255)
tATGAGTTTCGCGTTTAATTCGTTCTGAATCTTTTTTAATAAAGATGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCTAGTCAGCTGTTTCTCCGTAAAGCAAACCTGGAGGTATACAAAACATCAACGCAAc                                                                                                                                           >  2:367586/1‑141 (MQ=255)
   gAGTTTCGCGTCTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACAtt                                                                                                                                        <  1:415683/141‑1 (MQ=255)
        tccccgttAATTCGTTCCTGATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTa                                                                                                                                   <  1:45428/135‑1 (MQ=255)
        tCGCGTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGCATACCCAACATCAACGCAACATTTTTTa                                                                                                                                   >  2:53747/1‑141 (MQ=255)
         tgcgTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTAt                                                                                                                                  <  1:407268/140‑1 (MQ=255)
           cgTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATTCAATCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCACTCATCAACGCAACATTTTTTTTTg                                                                                                                                >  2:367325/1‑141 (MQ=255)
           cgTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACCTTTTTTATTg                                                                                                                                >  2:618859/1‑141 (MQ=255)
           cgTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTg                                                                                                                                <  1:5260/141‑1 (MQ=255)
           cgTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTg                                                                                                                                >  2:200634/1‑141 (MQ=255)
           cgTTTAATTCGTTCAGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAACAATCAACGCAACATTTTTTATTg                                                                                                                                >  1:437080/1‑141 (MQ=255)
            gTTTAATTCGTTCTGATTCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGc                                                                                                                               <  1:273213/141‑1 (MQ=255)
            gTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGc                                                                                                                               <  1:440864/141‑1 (MQ=255)
            gTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGc                                                                                                                               <  2:120396/141‑1 (MQ=255)
               tAATTCGTTGTGAATCTTTTTTAATAAAGGTGTTTTGGGACATTACTCGAATACTAGATAGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATCCCAATCATAAACGCAACATTTTTTATTGCGtc                                                                                                                            >  1:488786/1‑141 (MQ=255)
                 aTTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCAtctc                                                                                                                          >  1:212310/1‑141 (MQ=255)
                  ttCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCAtctct                                                                                                                         >  1:9720/1‑141 (MQ=255)
                  ttCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAACGTGTGCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCCTCAACTCAACATTTTTTATTGCAtctat                                                                                                                         >  2:98508/1‑139 (MQ=255)
                     gTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCATCATTTTTTATTGCATCTCTCAt                                                                                                                      >  2:295588/1‑141 (MQ=255)
                     gTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCAt                                                                                                                      >  1:442495/1‑141 (MQ=255)
                     gTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACAGCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCATCATTTTTTATTGCATCTCTCAt                                                                                                                      >  2:439713/1‑141 (MQ=255)
                      gtCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCAtt                                                                                                                     <  2:127492/140‑1 (MQ=255)
                            aTCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGctca                                                                                                               <  1:118516/141‑1 (MQ=255)
                             tCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTCTTATTGCATCTCTCATTAGctcac                                                                                                              <  1:487609/141‑1 (MQ=255)
                                tttttCATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGctcactca                                                                                                           <  1:200634/141‑1 (MQ=255)
                                   ttAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGTGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAAtt                                                                                                        <  1:403120/141‑1 (MQ=255)
                                    taataaAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTc                                                                                                       >  2:112574/1‑141 (MQ=255)
                                      ataaAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAg                                                                                                     <  1:143243/141‑1 (MQ=255)
                                      ataaAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAg                                                                                                     <  2:589510/141‑1 (MQ=255)
                                        aaaGTTGTTCAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCCGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGaa                                                                                                   <  2:361884/141‑1 (MQ=255)
                                               ttAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCCAAATAAGGCa                                                                                            <  1:489957/141‑1 (MQ=255)
                                               ttAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTGTTTATTGCACCTCCCATAAGCTCACACAATTCAGAAAAAGGCa                                                                                            >  1:10135/1‑141 (MQ=255)
                                                tAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAAAACGTCGCCTATCTGTAAATCAAACCTGGAGGGATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCAt                                                                                           <  1:151421/141‑1 (MQ=255)
                                                 aaaCAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATg                                                                                          <  1:333597/141‑1 (MQ=255)
                                                  aaCAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGt                                                                                         >  2:319747/1‑141 (MQ=255)
                                                  aaCAGAACTTCTCAAATACAAGCGTGTTCTTTTCACAATCAGCTGTTACTCCGTCAATCAAACATGGAAGCATAGCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCCGAACAAGGCATGt                                                                                         >  2:111235/1‑141 (MQ=255)
                                                    cAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTa                                                                                       >  1:417784/1‑141 (MQ=255)
                                                    cAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCACCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATCTTa                                                                                       >  1:606086/1‑141 (MQ=255)
                                                       aaCTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTa                                                                                    <  2:240117/141‑1 (MQ=255)
                                                          ttCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTCTCTCCGTAAATCAAACCTGGAGGTATACCAATCAACAACGCAACATTTTTTATTGCATCTCTCATTAGCTAAAAAAAAAAAAAAAAAGGCATGTTACTATTc                                                                                 <  2:503305/141‑1 (MQ=255)
                                                          ctcccAAATACAAGCGTGTTCTTCAACACATCAGCTGGCTCCCATTATATCAAACCTGCGTGTATACCAATCATCAACGCGATATTTTTTATTGCATCTCTCATTCGCTTAATCAACTCAGAAAAAGGCATGTTACTATTc                                                                                 <  2:564052/137‑1 (MQ=255)
                                                           tctcAAATCCAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCt                                                                                <  1:14154/141‑1 (MQ=255)
                                                           tctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCt                                                                                >  2:382951/1‑141 (MQ=255)
                                                           tctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCt                                                                                >  1:63164/1‑141 (MQ=255)
                                                           tctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCCCACTCAATTCAGACAAAGGCATGTTACTATTCt                                                                                >  2:336244/1‑141 (MQ=255)
                                                           tctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATAGCATGTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCt                                                                                >  1:124067/1‑141 (MQ=255)
                                                           tctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTACTCCGTAAAGCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCt                                                                                <  2:69661/141‑1 (MQ=255)
                                                            ctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAGGCCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCtt                                                                               <  1:549140/141‑1 (MQ=255)
                                                            ctcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCtt                                                                               >  1:498736/1‑141 (MQ=255)
                                                             tcAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTa                                                                              >  2:351687/1‑141 (MQ=255)
                                                             cccaacaCAACCGTGTTCTTTTCCAAATCAGCTGTTCCTCCGTAAATCAATCCTGGAGGGATACCAATCATCAACGCAACATTTGTTCTTGCATCTCGCCTTAGCGCACTCAATTCAGCAAAAGGCATGTTACTATTCTTa                                                                              <  2:562255/135‑1 (MQ=255)
                                                                  caCAAGCGTGTCCTTTTCCAAACCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCCCATTAGCTCCCACAATTCATAAAAAGGCATGTTACTATTCTTaccacc                                                                         <  1:660990/140‑1 (MQ=255)
                                                                      aGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTa                                                                     <  1:499639/141‑1 (MQ=255)
                                                                       gCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTaa                                                                    >  1:49547/1‑140 (MQ=255)
                                                                       gCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAg                                                                    >  1:90519/1‑141 (MQ=255)
                                                                       gCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGATAATATTCTTACCACCATTaa                                                                    >  2:211509/1‑140 (MQ=255)
                                                                        cGTGTTCTTTTCCAAATCAGCTGGTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGa                                                                   >  2:172678/1‑141 (MQ=255)
                                                                         gtgtTTTTTTCCAAATCAGCTGTTTCTCCGTAAACCAACCCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGAt                                                                  <  1:62360/141‑1 (MQ=255)
                                                                          tgctCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAACCCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATa                                                                 <  1:374982/138‑1 (MQ=255)
                                                                          ggtTCTTTGCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATa                                                                 <  1:81438/140‑1 (MQ=255)
                                                                             tCTTTTCCAAATCCGCTGTTGCTCCGTACAGCAAACCTGGAGGTTTCCCACTCATCACCGCAACATTTTTTATTGCATCTCTCATTTGCTCACTCAATGCAGCAAAAGGAATGGTACTATTCTTACCACCATTTGATAAtt                                                              <  2:182760/141‑1 (MQ=255)
                                                                             tCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAAtt                                                              >  1:653517/1‑141 (MQ=255)
                                                                             tCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTACATAAtt                                                              >  1:189237/1‑141 (MQ=255)
                                                                                tttCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCt                                                           >  1:510306/1‑141 (MQ=255)
                                                                                tttCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCt                                                           >  1:305235/1‑141 (MQ=255)
                                                                                  tCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCa                                                         <  2:343387/141‑1 (MQ=255)
                                                                                       aTCAGCTGTTTCTCCGTAAATCAAACCTGTAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATc                                                    <  1:193433/141‑1 (MQ=255)
                                                                                          aGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATCATTCCTCATAATCaaa                                                 >  2:21721/1‑141 (MQ=255)
                                                                                                  ctcCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCGATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTa                                         >  2:497275/1‑141 (MQ=255)
                                                                                                    ccGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATACCTCATAATCAAAACCTTCTATc                                       <  1:286105/141‑1 (MQ=255)
                                                                                                     cGTAAATCAAACCTGGAGGTATACCAATCATCAACGCATCATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCa                                      <  1:439713/141‑1 (MQ=255)
                                                                                                     cGTAAATCAAACCTGGAGGTATACCAATCATCAACGCATCATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCa                                      <  1:295588/141‑1 (MQ=255)
                                                                                                          aTCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGgc                                 <  1:653255/141‑1 (MQ=255)
                                                                                                              aaCCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGa                             >  2:236833/1‑141 (MQ=255)
                                                                                                              aaCCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCAGGTCGCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTCTCAAAGgcgcta                             >  2:439232/1‑139 (MQ=255)
                                                                                                                    gAGGTATACCAATCATCAACGCAACCTTGTTTGTTTCATGGCTCATTAGCTCACCCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTa                       <  1:88664/141‑1 (MQ=255)
                                                                                                                    gAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTa                       <  1:319747/141‑1 (MQ=255)
                                                                                                                        tataCCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAACCTTCTATCAAAGGCGCGATTGCTAGTTGa                  <  2:444378/141‑1 (MQ=255)
                                                                                                                         acaCCAATCATCAACGCAACATTTTTTATTTCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGa                  <  2:106705/139‑1 (MQ=255)
                                                                                                                             cAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGAtttt              >  1:62456/1‑141 (MQ=255)
                                                                                                                             cAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGAtttt              <  1:163089/141‑1 (MQ=255)
                                                                                                                              aaTCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGAttttt             <  1:282347/141‑1 (MQ=255)
                                                                                                                                 catcaACGCAACATTTTTTATTGCATCTAAATTTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGCCGCGATTGCTGGTTGATTTttatt          >  2:611918/1‑141 (MQ=255)
                                                                                                                                 catcaACAAAACATTTTTTATTGCGCCTTTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTttatt          <  2:498736/141‑1 (MQ=255)
                                                                                                                                   tcaACGCAACCTTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGTGCGATCGCTACTTGATTTTTATTaa        <  1:367586/141‑1 (MQ=255)
                                                                                                                                   tcaACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTaa        >  1:495898/1‑141 (MQ=255)
                                                                                                                                    caACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTaaa       >  1:308589/1‑141 (MQ=255)
                                                                                                                                    caACGCAACATTTTTTATTGCATCTCTCATTAGCTAAATCCATTCAGAAGAATGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTaaa       <  1:549301/141‑1 (MQ=255)
                                                                                                                                    caACGCAACATTTTTTATTGCAGCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTaaa       <  1:47809/141‑1 (MQ=255)
                                                                                                                                     aaCGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAACAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATGTTTATTaaaa      >  2:38038/1‑141 (MQ=255)
                                                                                                                                      aCGCAACATTTTTTATTGTATTTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTATTACCACCATTAGATATTTCGTCCTCATCCAAACCTTCTACATAAGTCTAGCTTACTTGTTGATTTTTCTTataat     >  2:647713/1‑140 (MQ=255)
                                                                                                                                       cGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCCGCAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCACAACCTTCTATCAAACGCGCGATTGCTAGTTGTTTTTTATTCACAGt    >  1:569251/1‑141 (MQ=255)
                                                                                                                                        gCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTAAAAGTa   <  1:644340/141‑1 (MQ=255)
                                                                                                                                         cAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTAAAAGTAt  >  1:425604/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                          
TATGAGTTTCGCGTTTAATTCGTTCTGAATCTTTTTTAATAAAGGTGTTAAACAGAACTTCTCAAATACAAGCGTGTTCTTTTCCAAATCAGCTGTTTCTCCGTAAATCAAACCTGGAGGTATACCAATCATCAACGCAACATTTTTTATTGCATCTCTCATTAGCTCACTCAATTCAGAAAAAGGCATGTTACTATTCTTACCACCATTAGATAATTCCTCATAATCAAAACCTTCTATCAAAGGCGCGATTGCTAGTTGATTTTTATTAAAAGTAT  >  CP000730/2107402‑2107679

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