Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 2,648,469 T→A 8.1% G79G (GGA→GGT USA300HOU_2501 ← MFS family major facilitator transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007302,648,4690TA 8.1% 190.9 / 2.8 67G79G (GGA→GGTUSA300HOU_2501MFS family major facilitator transporter
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base A (4/2);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACTTGAAAATCG  >  CP000730/2648330‑2648604
                                                                                                                                           |                                                                                                                                       
tATCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTc                                                                                                                                        <  1:506511/141‑1 (MQ=255)
tATCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTc                                                                                                                                        <  2:232889/141‑1 (MQ=255)
 aTCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGATTAATACCGCAAACAATTGCTTGAAATAACCTTCATATAGCACAAAATATAAAGTGATGTGCCACTGCACAACTAAAGCTAACTATGATAACTAttgc                                                                                                                                       >  2:752351/1‑139 (MQ=255)
 aTCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTcc                                                                                                                                       >  1:253062/1‑141 (MQ=255)
 aTCTGTACTTAGTTAGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGGCAGCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTcc                                                                                                                                       >  1:452188/1‑141 (MQ=255)
  tCTGTACGTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGTTTGAAATAACCTTGATATAGCAAATAATATAAAGCTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTccc                                                                                                                                      <  2:250949/141‑1 (MQ=255)
      tACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATc                                                                                                                                  <  2:270613/141‑1 (MQ=255)
         ttAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCg                                                                                                                               <  1:751341/141‑1 (MQ=255)
          tAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATAATTACTATTCCCAATccgtt                                                                                                                              <  1:343961/141‑4 (MQ=255)
             tttGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATCACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCAAAAATACAGCAAGATCTCATTAATCTTCCCAATACCATTaa                                                                                                                           >  2:455581/1‑141 (MQ=255)
               tGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTaaaa                                                                                                                         >  1:229355/1‑141 (MQ=255)
               tGGCAAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTaaaa                                                                                                                         <  1:645298/141‑1 (MQ=255)
                ggCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCACTTGCTTGAAATAACCTTGATATAGCAAATAATATAAATTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAg                                                                                                                        <  1:126261/141‑1 (MQ=255)
                ggCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAg                                                                                                                        <  1:576908/141‑1 (MQ=255)
                 gCTAATAAAGTTAAAGTCACGAGTGGAATAAATCCGACACTAATTGCTTGAAATAACCAAGATATAGCAAAAAATATATAGTTATGTGGCACGGGACACATAAAGCTCGATATGATTACTCTTAACAATCCACTTAAAAga                                                                                                                       >  2:45835/1‑140 (MQ=255)
                 gCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGt                                                                                                                       >  2:813789/1‑141 (MQ=255)
                   taataaAGTTAACGTCACGAGAGGAATAATTCCGCAACTATTTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTAc                                                                                                                     <  2:511769/141‑1 (MQ=255)
                   taataaAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGCTAAAAGTAc                                                                                                                     <  1:635232/141‑1 (MQ=255)
                   taataaAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATCATTACTATTCCCAATCCCGTTAAAAGTAc                                                                                                                     <  2:512025/141‑1 (MQ=255)
                     caaaagttaaagtCACGAGAGGAATAATTCCACAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCATTAAAAGTACtt                                                                                                                   <  2:527190/139‑1 (MQ=255)
                      taaagttaaagtCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACttt                                                                                                                  >  1:724870/1‑141 (MQ=255)
                            taaagtCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGa                                                                                                            >  2:445599/1‑141 (MQ=255)
                             aaagtCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGaa                                                                                                           >  1:816165/1‑141 (MQ=255)
                               agtCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATAGCCAATCCCGTTAAAAGTACTTTTTTCGATTc                                                                                                         >  1:106641/1‑141 (MQ=255)
                               agtCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATACCCAATCCCGTTAAAAGTACTTTTTTCGAATc                                                                                                         >  1:493112/1‑141 (MQ=255)
                                    cGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACtttttttttttttttat                                                                                                    >  2:779147/1‑130 (MQ=255)
                                    cGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACgag                                                                                                    <  2:660689/141‑1 (MQ=255)
                                      agagGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGaa                                                                                                  >  2:772713/1‑141 (MQ=255)
                                      agagGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGAGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAACGTACTTTTTTCGAATCACGAGaa                                                                                                  >  2:454796/1‑141 (MQ=255)
                                        acgAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATc                                                                                                <  2:363598/139‑1 (MQ=255)
                                         gggacaatGCCGCAACCAATTGCTTGAAATAACCCTTATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATGACTATTCCCAATCCCGTTAAAAGTACTGTTTTCGAATCACGAGAATCt                                                                                               <  2:283290/136‑1 (MQ=255)
                                         ggAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATCAAGGTCTGTGTCACTGACCAAATAAAGCTAGATATCATTACTATTACAAATCACATTAAAAGTAGTTGGGTCGAATCACCGGAAtat                                                                                               >  1:390012/1‑139 (MQ=255)
                                          gAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCta                                                                                              <  1:543956/141‑1 (MQ=255)
                                          gAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCta                                                                                              >  2:830481/1‑141 (MQ=255)
                                          gAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCta                                                                                              <  1:264510/141‑1 (MQ=255)
                                          gAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGATAGATATGATTACTATTCCCAATCCCGTTAAAAGTCCTTTTTTCGACTCCCGGGAAtata                                                                                              >  1:32204/1‑141 (MQ=255)
                                           aacaatTCCACAACCAATTCCTTGAAATAACCTTGATATATCAAATAATACAAAGTTATGTGTCACTGCACAAATACAACTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCtat                                                                                             <  2:4247/138‑1 (MQ=255)
                                                    gCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATCTGTCGCGTCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTAcc                                                                                    <  1:812367/141‑1 (MQ=255)
                                                      aaCCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCaa                                                                                  >  2:663871/1‑141 (MQ=255)
                                                           aTTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACa                                                                             >  2:565775/1‑141 (MQ=255)
                                                               cTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAAttt                                                                         <  2:248407/141‑1 (MQ=255)
                                                                 tGAAATAACCTAGATATAGCACATCATCTAAAGCTATAAGTCACTGCACAGATAAAGCTAGATATGATAACTCTACCCACTCCCGTTACAAGTACTTTTTTCGAATCACGCGAATCTATACTTCCACCAAATACAAtttcc                                                                       >  2:219291/1‑139 (MQ=255)
                                                                         ccTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCAATACTTCTACCAAGTACAATTTGCATAACTGc                                                               >  1:479612/1‑141 (MQ=255)
                                                                           ttGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCaa                                                             >  2:90884/1‑141 (MQ=255)
                                                                                 taGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATaa                                                       <  2:509413/141‑1 (MQ=255)
                                                                                    cAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTTTTTTTTTTTTTATCTAAAATTATACCAAGAACAATTTGCATGACTGCAAGGATaaaaa                                                    >  2:611596/1‑139 (MQ=255)
                                                                                           tataAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGGTAAATATAAAAcc                                             >  1:367882/1‑141 (MQ=255)
                                                                                              aaaGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTTTTTTTTTTTGAATCTCTACTCCTCCTAAGTAAAATTTGCATGACTCCAGCGTTCATTATGAATCCGcc                                          >  2:704852/1‑141 (MQ=255)
                                                                                                    ttgttgcccctcacaAATATAGATAGATATGCTTACTAGTCCCAGTCCCGTTAAAAGTACTTTTTTCGAGTCACTACCATCTATACTTCTACCAAGTACACTTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGt                                    <  2:468286/130‑1 (MQ=255)
                                                                                                     tgtgTCACTGCACCAATAAAGCTAGATATGATTACCAGTCCTAATCACGTTAAAAGTACTTATTTCGAATCACGCGAATCTATCCTTCGACCAGGCACAATTTGCATAACTAACACGATAAATATAAAACCGACAACTGTc                                   >  2:388483/1‑141 (MQ=255)
                                                                                                      gtgtaaatgaaaaaaTAAAGCTAAATATGATTACTATACCCAATCCCGTCAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCg                                  <  2:106641/129‑1 (MQ=255)
                                                                                                      gtgtCACTGCACAAATAAAGCTAGATATGATTACTATACCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCg                                  <  2:493112/141‑1 (MQ=255)
                                                                                                       tgtCACTGCACAAATAAAGCTAGGTATGATTACTATTCCCAACCCCGTTAAATGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGc                                 <  2:914/141‑1 (MQ=255)
                                                                                                       tgtCACTGCACAAATAAAGCTAGATATGATAAAACAAAAAAAAACCGTTAAAAGTACTTTTTTGGAAGCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCACCGATAAAAACAAAACCGCCacccgttgc                                 >  1:508800/1‑133 (MQ=255)
                                                                                                         tCATTGCACAAATAAAGCTCGAGTTGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTAGAATTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCAt                               <  2:210251/141‑1 (MQ=255)
                                                                                                           aCTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCACCGATAAATATAAAACCGCCAACTGTCCCAttt                             >  1:341629/1‑141 (MQ=255)
                                                                                                           aCTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCAttt                             >  1:341657/1‑141 (MQ=255)
                                                                                                                  aaaTAAAGCTAGATATGATTCCTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCTCGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAGACCGCCAACTGTCGCATTTATCCaaa                      >  1:527210/1‑141 (MQ=255)
                                                                                                                   aaTAGAGCTAGATATGATTACTATTCCCAATCCCGTTAACAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCaaaa                     <  2:483691/141‑1 (MQ=255)
                                                                                                                      aaaGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAg                  >  1:825684/1‑141 (MQ=255)
                                                                                                                        aGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCCACGATAAATGTTAAACCGCCAACTGTCGCATTTATCCAAAATAGTg                >  1:585991/1‑141 (MQ=255)
                                                                                                                           tAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTa             >  1:457448/1‑141 (MQ=255)
                                                                                                                              atatGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCATCGATAAATACAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACtt          <  1:593053/141‑1 (MQ=255)
                                                                                                                               tatGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACttt         >  2:469456/1‑140 (MQ=255)
                                                                                                                               tatGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACTTg         <  1:813789/141‑1 (MQ=255)
                                                                                                                                 tGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTGATCCAAAATAGTGGTACTTGaa       >  1:686194/1‑141 (MQ=255)
                                                                                                                                     tACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACTTGaaaacc   >  2:864969/1‑139 (MQ=255)
                                                                                                                                      aCTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACTTGAAAATCg  >  1:746173/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                       
TATCTGTACTTAGTTTGGCTAATAAAGTTAAAGTCACGAGAGGAATAATTCCGCAACCAATTGCTTGAAATAACCTTGATATAGCAAATAATATAAAGTTATGTGTCACTGCACAAATAAAGCTAGATATGATTACTATTCCCAATCCCGTTAAAAGTACTTTTTTCGAATCACGAGAATCTATACTTCTACCAAGTACAATTTGCATAACTGCAACGATAAATATAAAACCGCCAACTGTCGCATTTATCCAAAATAGTGGTACTTGAAAATCG  >  CP000730/2648330‑2648604

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: