Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,018,260 (A)8→7 100% 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.98.3% 529.4 / 1.7 115coding (308/444 nt)fliJflagellar protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (51/62);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 129 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348
                                                                                                   |                                                                                               
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:1007279/1‑100 (MQ=255)
  cagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaa                                                                                               >  2:136719/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              <  1:1157994/100‑1 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:218889/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:2829015/1‑100 (MQ=255)
    gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTACTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             >  1:141750/1‑100 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  1:1024137/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  1:2165891/100‑1 (MQ=255)
       tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          <  2:5399/100‑2 (MQ=255)
        aTCCAGACGATAAAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGAAAGAGAaaaaaaaa                                                                                           <  2:3046595/98‑1 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        <  2:490131/100‑4 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  1:1442858/1‑97 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGCTGGCGAGAAAAGAaaca                                                                                        <  1:311469/100‑1 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  2:630908/1‑96 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  2:1939335/1‑96 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGAAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     >  1:2289338/1‑100 (MQ=255)
               cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTg                                                                                 <  1:2847373/100‑1 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               <  2:459499/100‑1 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               >  1:1215734/1‑100 (MQ=255)
                  tGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAg                                                                              >  1:2342433/1‑100 (MQ=255)
                    gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                            >  2:2969282/1‑100 (MQ=255)
                    gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                            >  1:1375546/1‑100 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:3075924/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGAAACTTAATCAGTGGCCGCAGAAAGTTGACATTGCCCTGAACAATTAAAGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:2189734/100‑1 (MQ=255)
                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGCACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:1998655/99‑1 (MQ=255)
                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCa                                                                      <  2:2530030/100‑1 (MQ=255)
                          gCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCa                                                                      >  1:1662325/1‑100 (MQ=255)
                            cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                    <  2:822051/100‑1 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                   >  2:431333/1‑100 (MQ=255)
                              ttccTCAGCATCGTCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGan                                                                   <  2:59695/96‑2 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGTAAAAAACAACGTTTGCATGCCTGGCAGacgc                                                                 <  1:2187066/100‑3 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 <  1:198459/100‑1 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 >  1:387068/1‑100 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 <  1:2666343/100‑1 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  2:781777/1‑100 (MQ=255)
                                 cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                               >  1:2513984/1‑100 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  2:1391105/100‑1 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:2096339/100‑1 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             >  1:2630856/1‑100 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             >  2:1653416/1‑100 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             <  2:2605486/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:2434406/100‑1 (MQ=255)
                                      cATCGCCAACAACAAAATCAGTTTACGCAGAAAGTTGAGATTGCCCTGAACAGATGAAAAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 <  2:2594890/93‑1 (MQ=255)
                                       aTCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                         >  1:620817/1‑100 (MQ=255)
                                        tCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        >  1:2303689/1‑100 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                      >  2:1288487/1‑100 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       <  2:806222/99‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:49407/100‑1 (MQ=255)
                                            cAGCAACTTAATCAGTGGACGCAGAAAGTTGACGTTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGGAGGAACGGCa                                                    >  2:2136807/1‑100 (MQ=255)
                                            cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                    <  2:1007274/100‑1 (MQ=255)
                                            cAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCa                                                    <  1:977473/100‑1 (MQ=255)
                                              gCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAAt                                                  >  2:2450027/1‑100 (MQ=255)
                                                aaCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCTGGAACGGCAATc                                                 <  1:2128943/99‑1 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  2:1671297/1‑100 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              <  1:591673/100‑1 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:3230164/1‑100 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:3027075/1‑100 (MQ=255)
                                                   ttAATCCGTGGACGCAGAAAGTTGACATTGCCCTTAATAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACggt                                                     <  2:1107921/92‑2 (MQ=255)
                                                   ttAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcg                                             <  2:1442850/100‑1 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  1:2717848/1‑100 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTACCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCATGGCAGAAACTGCAGGAACGGCAATCCAcg                                             >  2:2931711/1‑99 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  1:3132510/100‑1 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           >  2:1155180/1‑100 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  1:1793689/100‑1 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           >  2:2024729/1‑100 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          <  1:1920860/100‑1 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        >  2:1639250/1‑100 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        <  2:3083488/100‑1 (MQ=255)
                                                         aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       <  2:141749/100‑1 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                      >  1:2581146/1‑100 (MQ=255)
                                                              aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  1:1653427/100‑1 (MQ=255)
                                                              aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  1:1657664/100‑1 (MQ=255)
                                                               cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                                 >  1:3009209/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:2055975/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  2:191752/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  2:1021896/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:2853634/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:2952021/1‑100 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              >  1:2057845/1‑100 (MQ=255)
                                                                   gAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                             <  1:1424395/100‑1 (MQ=255)
                                                                    aaaGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                            <  1:1136878/100‑1 (MQ=255)
                                                                      aGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          >  1:650819/1‑100 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    <  2:2798934/100‑1 (MQ=255)
                                                                              ttGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  <  1:1229072/100‑1 (MQ=255)
                                                                               tGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                 <  1:2545263/100‑1 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  2:2342417/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:567955/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:800469/100‑1 (MQ=255)
                                                                                  ccTGAACAGTGGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  2:2816789/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  1:2494990/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  2:2303673/100‑1 (MQ=255)
                                                                                    ttaaCAGTTTGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  1:307642/97‑1 (MQ=255)
                                                                                    tGAACCGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             >  2:2236803/1‑99 (MQ=255)
                                                                                     gAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           >  2:3047189/1‑100 (MQ=255)
                                                                                     gAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           >  1:1920318/1‑100 (MQ=255)
                                                                                      aaCAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaa          >  1:3004741/1‑100 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        >  1:1192922/1‑100 (MQ=255)
                                                                                           ttGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga     >  2:1317168/1‑100 (MQ=255)
                                                                                              gCGAGAAAAAAACAACTTTTGCAGTCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    <  1:931856/98‑1 (MQ=255)
                                                                                              gCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  >  1:3097985/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 < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: