Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,018,260 (A)8→7 coding (315/444 nt) fliJ → flagellar protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,018,2530A.95.0% 408.7 / 13.6 99coding (308/444 nt)fliJflagellar protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base . (49/45);  total (51/48)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 112 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348
                                                                                                   |                                                                                               
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:1123951/1‑100 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  1:1139056/1‑100 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCAGTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:2461033/1‑100 (MQ=255)
    gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             <  1:1277514/100‑1 (MQ=255)
     tttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGTCATTGCCCTGAACAGTTGGCGAGaaaaaa                                                                                            <  2:2556306/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  1:2388747/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  1:873351/100‑1 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        <  2:354248/100‑4 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  2:3105520/1‑97 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        <  2:282828/100‑4 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  2:2106138/1‑96 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     >  2:770226/1‑100 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     >  2:2737215/1‑100 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     >  1:720690/1‑100 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGACAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGGCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     >  1:3182980/1‑100 (MQ=255)
                gCTGGAAGAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGc                                                                                >  1:2540726/1‑100 (MQ=255)
                gCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCCGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGc                                                                                >  2:1835585/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               <  1:2838283/100‑1 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               >  1:6888/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               <  1:2106147/100‑1 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               >  1:1432026/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               >  2:438198/1‑100 (MQ=255)
                  tGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAg                                                                              >  2:2444359/1‑100 (MQ=255)
                     aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGcc                                                                           <  2:1123949/100‑1 (MQ=255)
                      aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCt                                                                          >  1:16343/1‑100 (MQ=255)
                      aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCt                                                                          >  1:3285615/1‑100 (MQ=255)
                      aaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGCGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                            >  2:750072/1‑98 (MQ=255)
                       aaaGCCATTACTCAGCATCGCCAGCAACTTAATCTGTGGACGCAGAAAGTTGACTTTGCCCTGAACAGTTGACAAAAAAAAAACAACGTTTGCAGGCCTg                                                                         <  2:1138422/100‑1 (MQ=255)
                       aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTg                                                                         <  2:459247/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:2418392/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:516304/100‑1 (MQ=255)
                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGTACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCa                                                                      <  2:2007657/100‑1 (MQ=255)
                           ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAg                                                                     <  2:2093291/100‑1 (MQ=255)
                            cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                    >  2:797912/1‑100 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTTAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCATac                                                                   <  1:714992/100‑1 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                   <  2:2500303/100‑1 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 >  2:1535347/1‑100 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  2:707607/1‑100 (MQ=255)
                                 cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                               >  1:1543400/1‑100 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  2:3170472/100‑1 (MQ=255)
                                   tccgaTCTCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             <  2:1485024/96‑1 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             <  1:1913279/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:2759551/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            >  1:93316/1‑100 (MQ=255)
                                       aTCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                         >  2:625606/1‑100 (MQ=255)
                                       aTCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                         <  1:1705111/100‑1 (MQ=255)
                                        tCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        >  1:2022107/1‑100 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAGACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       <  2:593946/100‑1 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      <  2:121578/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  1:2558101/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2575250/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2572068/100‑1 (MQ=255)
                                              gCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAt                                                  >  1:263776/1‑100 (MQ=255)
                                              gCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAt                                                  >  2:1847009/1‑100 (MQ=255)
                                                aaCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATcc                                                >  2:1509650/1‑100 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          >  1:2938930/1‑100 (MQ=255)
                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCCGGCCCGGCAGACACTGCCGGAACGGCCATCCAcgc                                            >  2:2940889/1‑96 (MQ=255)
                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         >  2:1704309/1‑100 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        >  1:1091843/1‑100 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        >  2:708542/1‑100 (MQ=255)
                                                         aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       <  1:3118090/100‑1 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                      <  2:1881246/100‑1 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                      <  1:2461044/100‑1 (MQ=255)
                                                           tGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     >  2:583759/1‑100 (MQ=255)
                                                           tGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACt                                     <  1:1163129/100‑1 (MQ=255)
                                                            ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGACGACAAAAAATACAAAGTTTGCAGGCCTGGCAGAAACGGCAGGAACAGCAATTCACGGGGGCa                                       >  2:2472458/1‑98 (MQ=255)
                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                   >  2:1486988/1‑100 (MQ=255)
                                                               cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                 <  2:3075257/100‑1 (MQ=255)
                                                               cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                 >  1:1729/1‑100 (MQ=255)
                                                                gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTg                                <  1:1635244/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:55709/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  2:3184239/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  2:3082287/1‑100 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:625607/100‑1 (MQ=255)
                                                                   gAAAGTTGACAGTGCCCGGAACAGTGGGCGAGATAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:503423/99‑1 (MQ=255)
                                                                    aaaGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                            <  1:2818317/100‑1 (MQ=255)
                                                                      aGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  1:2355359/1‑100 (MQ=255)
                                                                       gTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  1:907567/100‑1 (MQ=255)
                                                                        ttGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        <  2:2485459/100‑1 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       <  2:1707779/100‑1 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       >  2:6900/1‑100 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       <  2:2484830/100‑1 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAATGCTTGCAGAAAAcc                       >  1:3299749/1‑100 (MQ=255)
                                                                          gACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                      >  1:1651247/1‑100 (MQ=255)
                                                                          gACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                      <  1:8861/100‑1 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                     >  1:1071268/1‑100 (MQ=255)
                                                                             aTTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                   >  2:1238054/1‑100 (MQ=255)
                                                                              ttGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  <  2:40404/100‑1 (MQ=255)
                                                                                 cccGGAACAGTTCGCGAAAAAAAAACAACGTTTGCAGTCCGTTCAGACACTGCAGAAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCCAt               <  2:3182962/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:49220/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:931547/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  1:3169328/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             >  1:699063/1‑100 (MQ=255)
                                                                                    tGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  1:876411/100‑1 (MQ=255)
                                                                                     gAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           <  1:6900/100‑1 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        >  1:392194/1‑100 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaaa       >  1:2285222/1‑100 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      <  1:1333095/100‑1 (MQ=255)
                                                                                           ttGGCGAGAAAAAAGCAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      <  2:1334890/99‑1 (MQ=255)
                                                                                            tGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    >  1:1043771/1‑100 (MQ=255)
                                                                                                   |                                                                                               
AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348

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