Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,386,339 T→C 4.7% K41E (AAA→GAA)  USA300HOU_2245 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,386,3390TC 4.7% 248.8 / 2.0 87K41E (AAA→GAA) USA300HOU_2245hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/51);  new base C (2/2);  total (33/54)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-01

AGCCTTGCTCATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTA  >  CP000730/2386201‑2386476
                                                                                                                                          |                                                                                                                                         
aGCCTTGCTCATACTTCTCTAGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTGCTTCATCTAGTAGCTCAATATAAGTTTCATTGATATTACCTAATTTAt                                                                                                                                         >  1:639508/1‑141 (MQ=255)
      tcTCATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATg                                                                                                                                   <  1:826082/140‑1 (MQ=255)
      gCTCACACGTCTCGGGGACGATTTGAGGCGCACAAGGGGCTCCGCCTTATATCCTTTTTGCCATCCTTTTCAACTTTTCATCTTTTTCTAAATCTCGTAGCTCAATATAATTTTTATTGATATAACCTAATTTATTTAATg                                                                                                                                   <  2:390117/141‑1 (MQ=255)
       cTCCTACTTCCCCTGAACAATTTGAGTAGCAAAAGCCGCTCCGCCTTCTATCGTTTTTGCCATCCTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAATCTAATTTATTGAATGt                                                                                                                                  <  1:464779/141‑1 (MQ=255)
        tCATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGtt                                                                                                                                 >  2:126985/1‑141 (MQ=255)
         cATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTc                                                                                                                                <  2:627083/141‑1 (MQ=255)
         cATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTc                                                                                                                                <  2:36056/141‑1 (MQ=255)
         cATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTc                                                                                                                                <  2:370582/141‑1 (MQ=255)
          aTACGTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTTTCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATAGAATTTTCATTGATATAACCTACTTTATTGAATGTTcc                                                                                                                               >  1:641538/1‑141 (MQ=255)
            aCTTCTCTTGAGCAATTTGAATAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATC‑TTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTGAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCata                                                                                                                            >  1:640513/1‑141 (MQ=255)
             cTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTCATTGAATGTTCCata                                                                                                                            >  2:381995/1‑141 (MQ=255)
             cTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTCATTGAATGTTCCata                                                                                                                            >  2:128429/1‑141 (MQ=255)
               tctcTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATAtt                                                                                                                          >  2:492218/1‑141 (MQ=255)
               tctcTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTAATGAATGTTCCATAtt                                                                                                                          <  1:677460/141‑1 (MQ=255)
               tctcTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTATCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATCATTTTCATTGATATAACCTAATTAATTGAATGTTCCAGAtt                                                                                                                          >  2:608986/1‑141 (MQ=255)
                ctctTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCGAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATAttt                                                                                                                         <  1:133315/141‑1 (MQ=255)
                   ttGAACAATTTGAGTCGCAAAAGCGAGTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTGGTAGCGCAATATATTTTTCATTGATATAACCTAATTTCTTGAATGTTCCATATTTatg                                                                                                                      <  1:552765/141‑1 (MQ=255)
                     gAACAATTTGAGTAGCAAAAGCGCCTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTatgat                                                                                                                    <  2:523201/141‑1 (MQ=255)
                     gAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTatgat                                                                                                                    <  1:87147/141‑1 (MQ=255)
                     gAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTatgat                                                                                                                    >  1:23602/1‑141 (MQ=255)
                      aaCAATTTGAGTAGCAAAAGCGCCTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTATCTCAATATAATTTTAATTGATATAACCTAATTTATTGAATGTTCCATATTTatgatg                                                                                                                   <  2:232231/141‑1 (MQ=255)
                      aaCAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATTTAATTTTCATTGATATAACCTAATTTATTGCATGTTCCATATTTatgatg                                                                                                                   <  2:671057/141‑1 (MQ=255)
                      aaCAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTatgatg                                                                                                                   <  2:618451/141‑1 (MQ=255)
                       aCAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTAtgatga                                                                                                                  >  1:634614/1‑141 (MQ=255)
                       aCAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTAtgatga                                                                                                                  >  1:174937/1‑141 (MQ=255)
                       aCAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTAtgatga                                                                                                                  >  1:196152/1‑141 (MQ=255)
                           tgtgggTAGCAACAGCGACTCCGCCTTCTATCGTTTTTATCATCTTTTTCCACTTTCCATCATTTTCTACATCTAGTATCTAACTGTAAGTGTCATTGTTATAACCTAATTTATTAAATTTTCCAGATGTATGATGATAcc                                                                                                              <  2:21963/136‑1 (MQ=255)
                              gAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCAc                                                                                                           >  2:515781/1‑141 (MQ=255)
                              gAGTAGCAAAAGCGACTCCGCCTTCTATCGGTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCAc                                                                                                           >  2:516332/1‑141 (MQ=255)
                               cgTAGCAAAAACTACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCAcc                                                                                                          <  2:206619/140‑1 (MQ=255)
                               aGTAGCAAAAGCGTCTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCAcc                                                                                                          <  2:686504/141‑1 (MQ=255)
                               aGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTGATTGAATGTTCCATATTTATGATGATACCCAcc                                                                                                          <  2:59029/141‑1 (MQ=255)
                                 tAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCCAATTTATTGAATGTTCCATATTTATGATGATACCCACCTg                                                                                                        <  2:572874/141‑1 (MQ=255)
                                        aGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTaa                                                                                                 <  1:247515/141‑1 (MQ=255)
                                            aCTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAAtttt                                                                                             >  2:911089/1‑141 (MQ=255)
                                            aCTCCGCCTTCTATCGTTTTTGCCATCTTTGTTAAGTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGGTCCATATTTATGATGATACCCACCTGAATGGAAtttt                                                                                             >  1:216997/1‑141 (MQ=255)
                                             cTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTta                                                                                            >  1:220707/1‑141 (MQ=255)
                                                cGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGGATTTACTTTTATaa                                                                                         >  1:634642/1‑141 (MQ=255)
                                                cGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATaa                                                                                         <  1:806063/141‑1 (MQ=255)
                                                  ccTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACa                                                                                       <  2:107316/141‑1 (MQ=255)
                                                      cTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACAtctc                                                                                   <  1:526336/141‑1 (MQ=255)
                                                       tATCGTTTTTGCCATCTTTTACAACTTTTCCTTATTTTCTACATCGAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATTACATCTcc                                                                                  <  2:872027/141‑1 (MQ=255)
                                                        aTCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCt                                                                                 <  2:322192/141‑1 (MQ=255)
                                                        aTCGTTCTTGCTATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTACTTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCt                                                                                 <  2:896942/141‑1 (MQ=255)
                                                         tCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTg                                                                                >  2:613972/1‑141 (MQ=255)
                                                             ttttGTCATCTTGTTCAACTTTTCATTATTTACTACAGCTAGTAGATCAATATAATTTTCTTTGATATAACCAAATTTAATGAATAATCCATATTTCTGATTATATCCATCTGAATGTAATTTTATAACATCTCCTCGaaa                                                                            <  2:179638/141‑1 (MQ=255)
                                                             ttttGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATACTTTTCATTGATATAACCTAATTTATTGCATGTTCCCTATTTATTATGATACCCACCTGAATGTAATTTTATAACATCTCCTGCaaa                                                                            <  2:179751/141‑1 (MQ=255)
                                                                tGCCATCTTTTTCCACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACt                                                                         <  2:767640/141‑1 (MQ=255)
                                                                   cATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTGATTTTATTACATCTCCTGGAAAACTaaa                                                                      <  1:589529/141‑1 (MQ=255)
                                                                   cATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTaaa                                                                      <  1:134416/141‑1 (MQ=255)
                                                                    aTCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAAc                                                                     >  1:824077/1‑141 (MQ=255)
                                                                         tttCAACTTTTCTATATTTTCGACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATc                                                                <  1:245676/141‑1 (MQ=255)
                                                                         tttCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAGCTAAACCGCTc                                                                <  2:199590/141‑1 (MQ=255)
                                                                         tttCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATc                                                                <  2:279501/141‑1 (MQ=255)
                                                                         tttCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATc                                                                <  1:796733/141‑1 (MQ=255)
                                                                         tttCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATc                                                                <  2:521313/141‑1 (MQ=255)
                                                                          ttCACCTTTTCATTATTTTCTGCATCTATTAGCTCAATATAATTTTCACTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGTATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCt                                                               <  1:123613/141‑1 (MQ=255)
                                                                           tCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTa                                                              <  2:143466/141‑1 (MQ=255)
                                                                             aaCTCTTCCATATGTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTACTTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAAc                                                            <  1:15018/141‑1 (MQ=255)
                                                                               cTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACATAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACAGCTCCTGGAAAACTAAACCGATCTAACTg                                                          >  1:77560/1‑141 (MQ=255)
                                                                                   tCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGCTATAACCTAATTCATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTAGAACATCTCCTGGAAAACTAAACCGATCTAACTGATCt                                                      <  1:128429/141‑1 (MQ=255)
                                                                                   tCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTCATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCt                                                      <  1:381995/141‑1 (MQ=255)
                                                                                    tcttattTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCta                                                     <  2:6628/139‑1 (MQ=255)
                                                                                    cATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCta                                                     <  2:434668/141‑1 (MQ=255)
                                                                                       tattTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCtatat                                                  <  1:284528/141‑1 (MQ=255)
                                                                                        attTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCtatata                                                 >  1:157327/1‑141 (MQ=255)
                                                                                        attTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTAAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCtatata                                                 >  1:605559/1‑141 (MQ=255)
                                                                                         ttttCTACATCTAGTAGCTCAATATAATCTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTCTAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATaa                                                <  1:515781/141‑1 (MQ=255)
                                                                                              tACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATaatgaat                                           >  1:877614/1‑141 (MQ=255)
                                                                                                 aTCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGCTCATTATAATAAatgat                                        >  2:738508/1‑141 (MQ=255)
                                                                                                 aTCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATATTGGaggat                                        >  2:635626/1‑141 (MQ=255)
                                                                                                 aTCGAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGCATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCACCTGGAAAACTAAACCGATCTAACTGATCTATATAACGAatgat                                        <  2:641538/141‑1 (MQ=255)
                                                                                                    tagtagCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGAtgtg                                     <  2:717408/141‑1 (MQ=255)
                                                                                                      ttatcTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATTTAATTTTATAACATCTCCTTGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGAt                                   <  2:600698/137‑1 (MQ=255)
                                                                                                      gtagCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAatgattttat                                   >  2:13780/1‑136 (MQ=255)
                                                                                                             aaTATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAAttt                            <  1:126985/141‑1 (MQ=255)
                                                                                                             aaTATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAAttt                            >  1:188100/1‑141 (MQ=255)
                                                                                                                 taATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCact                        >  1:887289/1‑139 (MQ=255)
                                                                                                                   aTTTTCATTGATATAACCTAATTTATTGAATGTTACATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAatat                      >  2:374709/1‑141 (MQ=255)
                                                                                                                        taTTGATATAACCTAATTTATTGAATGTTCCATATTGATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTa                 <  2:482454/140‑1 (MQ=255)
                                                                                                                        cATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGACACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTa                 <  2:845449/141‑1 (MQ=255)
                                                                                                                             atatAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATc            >  1:586574/1‑141 (MQ=255)
                                                                                                                              tataACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATcc           <  2:634614/141‑1 (MQ=255)
                                                                                                                              tataACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATcc           <  2:174937/141‑1 (MQ=255)
                                                                                                                                 aaCCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGAGATCAAATTTCAATATCTCTAGTATccccc        >  2:47273/1‑141 (MQ=255)
                                                                                                                                  aCCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCt       >  1:864692/1‑141 (MQ=255)
                                                                                                                                       aTTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTa  <  2:371384/141‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AGCCTTGCTCATACTTCTCTTGAACAATTTGAGTAGCAAAAGCGACTCCGCCTTCTATCGTTTTTGCCATCTTTTTCAACTTTTCATTATTTTCTACATCTAGTAGCTCAATATAATTTTCATTGATATAACCTAATTTATTGAATGTTCCATATTTATGATGATACCCACCTGAATGTAATTTTATAACATCTCCTGGAAAACTAAACCGATCTAACTGATCTATATAATGAATGATGTGATCAAATTTCAATATCATTAGTATCCCCCTATTTA  >  CP000730/2386201‑2386476

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