Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,019,158 G→A 73.5% S258N (AGT→AAT)  fliK → flagellar hook‑length control protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,019,1580GA73.5% 175.9 / 82.3 117S258N (AGT→AAT) fliKflagellar hook‑length control protein
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base A (38/48);  total (54/63)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 117 total aligned reads displayed.

TCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCA  >  NC_000913/2019062‑2019254
                                                                                                |                                                                                                
tCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGc                                                                                               <  2:3226376/100‑1 (MQ=255)
 cgccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCa                                                                                              >  1:3940804/1‑100 (MQ=255)
  gccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCag                                                                                             >  1:682829/1‑100 (MQ=255)
  gccgcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCag                                                                                             <  2:1014079/100‑1 (MQ=255)
  gccgcACCTGTTTTGAGTGCACCGCTCGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATCTCGCTGTTCACCCGCCAGGGGCAACACAATGCa                                                                                              >  2:2766541/1‑99 (MQ=255)
     gcACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGt                                                                                          <  2:2316891/100‑1 (MQ=255)
       aCCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTg                                                                                        <  2:1459726/100‑1 (MQ=255)
       aCCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTg                                                                                        <  1:3331297/100‑1 (MQ=255)
       aCCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTg                                                                                        <  1:2374427/100‑1 (MQ=255)
        ccTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGc                                                                                       >  2:119623/1‑100 (MQ=255)
          tGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGTATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGt                                                                                     >  2:3220445/1‑100 (MQ=255)
          tGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGt                                                                                     >  1:1056171/1‑100 (MQ=255)
             tttGAGTGCACCGCTGGGTTCTCACTAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTg                                                                                  <  2:3366816/100‑1 (MQ=255)
             tttGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCagacgcaactctg                                                                                  <  2:3855124/100‑11 (MQ=255)
                 aGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCAc                                                                               <  1:3946282/99‑1 (MQ=255)
                  gTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCAccc                                                                             >  1:1204901/1‑100 (MQ=255)
                  gTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCAccc                                                                             >  2:1813329/1‑100 (MQ=255)
                  gTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCAccc                                                                             <  1:2258201/100‑1 (MQ=255)
                   tGCATCACTGTGTTCTCACGAATGGCAACTATCATTAAGCTAGCATATTTCGCTGTTCACCCGCCATGGGCAACTATATGCAGAGTTGCGTCTGCACCCg                                                                            <  2:1505889/100‑1 (MQ=255)
                    gCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGc                                                                           >  1:1754826/1‑100 (MQ=255)
                    gCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGc                                                                           >  2:2945497/1‑100 (MQ=255)
                     cACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCa                                                                          >  1:3406355/1‑100 (MQ=255)
                       ccGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAgg                                                                        <  1:196191/100‑1 (MQ=255)
                        cGCTGGATTCTCACGAATGCCAACATTAATTAAGCCCGCATCTTTCGCGGTTCACCCGCCAGGCGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAgg                                                                        <  1:2483521/99‑1 (MQ=255)
                         gCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGAt                                                                      >  1:3677767/1‑100 (MQ=255)
                          cTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGAtt                                                                     <  2:2242315/100‑1 (MQ=255)
                            gggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTa                                                                   >  2:2972846/1‑100 (MQ=255)
                            gggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTa                                                                   <  2:2796219/100‑1 (MQ=255)
                             ggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAg                                                                  <  2:1204896/100‑1 (MQ=255)
                             ggTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTCCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAg                                                                  <  2:3889659/100‑1 (MQ=255)
                             ggTTCTCACGAATGGCAACAATCATTAAACCAGCATATTTCGCTGTTCACCCGCCATGGACAACAAAAAGCAGAGTTGCGTCTGCACCCGCAGGATTTAg                                                                  >  2:1819954/1‑100 (MQ=255)
                              gTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAgg                                                                 >  2:2026332/1‑100 (MQ=255)
                              gTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAgg                                                                 <  1:425198/100‑1 (MQ=255)
                              gTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAgg                                                                 <  2:1258464/100‑1 (MQ=255)
                              gTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTAAgg                                                                 >  2:554820/1‑100 (MQ=255)
                               ttCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGt                                                                >  2:2121399/1‑100 (MQ=255)
                                tctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTg                                                               >  1:2433443/1‑100 (MQ=255)
                                 ctcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGCGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTg                                                               >  2:3541177/1‑99 (MQ=255)
                                  tctcGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGaa                                                             <  1:1634408/97‑1 (MQ=255)
                                  tcACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGaa                                                             >  1:1346823/1‑100 (MQ=255)
                                   cACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAg                                                            >  2:2026819/1‑100 (MQ=255)
                                   cACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAg                                                            <  1:3298007/100‑1 (MQ=255)
                                   cACGAATGGCAACAATAATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAAAAAAGTTAAGAGTTGCGTCTGCACCCGCAGGATTTAGGTg                                                               >  2:144475/1‑97 (MQ=255)
                                    aCGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGt                                                           <  2:1324329/100‑1 (MQ=255)
                                     cGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTg                                                          <  1:2972863/100‑1 (MQ=255)
                                     cGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGaaaga                                                          >  2:1634399/1‑97 (MQ=255)
                                         tGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCaaa                                                      <  2:1518297/100‑1 (MQ=255)
                                          ggCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAt                                                     >  2:1747194/1‑100 (MQ=255)
                                           gCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAtc                                                    >  2:390277/1‑100 (MQ=255)
                                            caacaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAtct                                                   <  2:150181/100‑1 (MQ=255)
                                             aacaaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAAtctc                                                  >  1:3601100/1‑100 (MQ=255)
                                               caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                >  2:4120585/1‑100 (MQ=255)
                                               caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                <  2:2021110/100‑1 (MQ=255)
                                               caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                <  1:3198081/100‑1 (MQ=255)
                                               caaTCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTccc                                                <  2:39795/100‑1 (MQ=255)
                                                      tAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGt                                         <  2:233262/100‑1 (MQ=255)
                                                         gCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGga                                      <  1:3541204/100‑1 (MQ=255)
                                                          ccAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAGGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGTAAATCTCCCTCAAAGTGgat                                     >  1:1625817/1‑100 (MQ=255)
                                                          ccAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgat                                     <  1:4160328/100‑1 (MQ=255)
                                                          ccAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgat                                     <  1:1208462/100‑1 (MQ=255)
                                                            aGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatga                                   >  2:137884/1‑100 (MQ=255)
                                                             gCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGgatgat                                  <  1:1902853/100‑1 (MQ=255)
                                                             gCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCTAATCTCCCTCAAAGTGgatgat                                  <  2:3180938/100‑1 (MQ=255)
                                                               atatTTCGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATaa                                <  2:3141065/100‑1 (MQ=255)
                                                                  tgtCGCTTTTCTCCCGCCTGGGCCAACAAAGTGCAGAGTTGAGGCTGCTCCGGCAGGATTTATGGGAAGTGCAAATCTCCCTCAAAGTGTAGGATAAc                               <  2:1359200/96‑1 (MQ=255)
                                                                   ttCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAg                            >  2:1448586/1‑100 (MQ=255)
                                                                     cGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgc                          <  2:2000683/100‑1 (MQ=255)
                                                                     cGCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgc                          >  2:1701260/1‑100 (MQ=255)
                                                                      gCTGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgcg                         >  1:829803/1‑100 (MQ=255)
                                                                      gCTGTGCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCGGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGgc                          <  2:3400607/99‑1 (MQ=255)
                                                                        tGTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCa                       <  2:2894310/100‑1 (MQ=255)
                                                                         gTTCACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCaa                      >  2:2007987/1‑100 (MQ=255)
                                                                            cACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATc                   <  1:4005981/100‑1 (MQ=255)
                                                                            cACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATc                   >  1:524935/1‑100 (MQ=255)
                                                                            cACCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATc                   >  2:1286641/1‑100 (MQ=255)
                                                                             aCCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATcc                  >  1:789148/1‑100 (MQ=255)
                                                                             aCCCGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATc                   >  2:2441402/1‑99 (MQ=255)
                                                                              cccGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCa                 >  2:2651222/1‑100 (MQ=255)
                                                                              cccGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCa                 >  1:3337051/1‑100 (MQ=255)
                                                                              cccGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCa                 <  1:2180937/100‑1 (MQ=255)
                                                                              cccGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTAATGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCa                 >  2:1417716/1‑100 (MQ=255)
                                                                               ccGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCATGCGCAAATCCAg                <  1:2736960/100‑1 (MQ=255)
                                                                                cGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGa               >  2:2004072/1‑100 (MQ=255)
                                                                                cGCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGa               <  1:2796710/100‑1 (MQ=255)
                                                                                 gCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAACTCTCCCTCCAAGTGGATGATAACCAGGCGCAAATCCAGAt              >  1:179087/1‑100 (MQ=255)
                                                                                 gCCAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGAt              >  2:3482213/1‑100 (MQ=255)
                                                                                  ccAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATg             <  1:2407942/100‑1 (MQ=255)
                                                                                  ccAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATg             <  1:2278366/100‑1 (MQ=255)
                                                                                   cAGGGGCAACAAAGTGCAGAGTTGCGTCTGAACCAGCAGGATTTAGATGAAGCGAAAATCTCCCACAAAGTGGATGATAACCAGGCGCAAATCCAg                <  2:615524/96‑1 (MQ=255)
                                                                                   cAGGGGCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATgg            <  1:3743996/100‑1 (MQ=255)
                                                                                     ggggCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGtt          <  1:716183/100‑1 (MQ=255)
                                                                                     ggggCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGtt          >  1:158243/1‑100 (MQ=255)
                                                                                      gggCAACAAAGTGCAGAGTTGCGTCTACACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGttt         >  1:2544215/1‑100 (MQ=255)
                                                                                       ggCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTc        <  1:2650034/100‑1 (MQ=255)
                                                                                        gCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa       <  2:668327/100‑1 (MQ=255)
                                                                                        gCAACAAAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCa       <  2:2197520/100‑1 (MQ=255)
                                                                                         caacaaAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCAc      <  2:3774793/100‑1 (MQ=255)
                                                                                           acaaAATGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCAcc     >  2:1979075/1‑99 (MQ=255)
                                                                                            caaAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGc   >  1:936181/1‑100 (MQ=255)
                                                                                             aaaGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCa  <  1:2367992/100‑1 (MQ=255)
                                                                                                |                                                                                                
TCGCCGCACCTGTTTTGAGTGCACCGCTGGGTTCTCACGAATGGCAACAATCATTAAGCCAGCATATTTCGCTGTTCACCCGCCAGGGGCAACAAAGTGCAGAGTTGCGTCTGCACCCGCAGGATTTAGGTGAAGTGCAAATCTCCCTCAAAGTGGATGATAACCAGGCGCAAATCCAGATGGTTTCACCGCA  >  NC_000913/2019062‑2019254

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