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.90.8% 520.8 / 47.5 142coding (308/444 nt)fliJflagellar protein
Reads supporting (aligned to +/- strand):  ref base A (4/9);  new base . (71/58);  total (75/67)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348
                                                                                                   |                                                                                               
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 <  2:2699205/100‑1 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:391185/1‑100 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACATAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:3273931/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:3632114/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              >  1:755497/1‑100 (MQ=255)
    gTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaa                                                                                             >  1:1649523/1‑100 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  2:754671/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAAGCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  2:3605582/100‑1 (MQ=255)
       tATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaac                                                                                          <  1:1727612/100‑2 (MQ=255)
        aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaaca                                                                                         <  2:4002647/100‑3 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  2:1407155/1‑97 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  1:2506310/1‑97 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  1:1398608/1‑97 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTACTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAaccac                                                                                       >  1:3506840/1‑98 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                     <  2:3060542/100‑1 (MQ=255)
             gACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTT                                                                                   >  2:1534668‑M1/1‑100 (MQ=255)
              aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTT                                                                                  >  1:1657555‑M1/1‑100 (MQ=255)
              aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTT                                                                                  >  2:3355628‑M1/1‑100 (MQ=255)
                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAgg                                                                             >  2:1025138/1‑100 (MQ=255)
                    ggaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAgg                                                                             >  1:684259/3‑99 (MQ=255)
                     aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGcc                                                                           >  2:1844108/1‑100 (MQ=255)
                     aaaaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGAAAAAAAAAAAAAAACAACGTTTGCAGGcc                                                                           <  1:4336271/100‑1 (MQ=255)
                       aaaGGCATTACTCAGCATCGCCAGCAACTTAATCAGGGGACGCAGAAAGTGGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:2110788/100‑1 (MQ=255)
                       aaaGCCATTACTCAGCATCGCCAGCAACCTAATCAGTAGACGCAGAAAGTTGACATTTCCCTGAACAAAAAAAAAAAAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:3803184/100‑1 (MQ=255)
                        aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTgg                                                                        >  2:1772873/1‑100 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:660241/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:1448719/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAG‑‑AAAAAACCACGTTTGCAGGCCttg                                                                        >  2:3931903/1‑96 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                   >  1:266852/1‑100 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  2:2960815/1‑100 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  2:2432177/1‑100 (MQ=255)
                                 cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                               >  2:424503/1‑100 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:3165617/100‑1 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTAGAGAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:1011245/100‑1 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             <  2:2031035/100‑1 (MQ=255)
                                      cATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGa                                                          >  2:1241421/1‑100 (MQ=255)
                                      cATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGa                                                          <  2:1297060/100‑1 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGGAGGAACg                                                       <  2:2321312/100‑1 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        <  1:2765261/99‑1 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGg                                                      <  1:1780328‑M2/100‑2 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGg                                                      >  2:2215346‑M2/1‑99 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGgc                                                     >  2:248157‑M2/1‑98 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGg                                                      <  1:2972191‑M2/99‑2 (MQ=255)
                                               cAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATc                                                 >  2:800846/1‑100 (MQ=255)
                                               cAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATc                                                 >  1:2360948/1‑100 (MQ=255)
                                                aaCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATcc                                                <  1:1582072/100‑1 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  1:2425907/1‑100 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  1:2754415/1‑100 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  1:1257638/1‑100 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  2:275555/1‑100 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCGGAACAGTGGGCGAGAAAAAAACAACGTTTGCcgttcctgcagtgtctgccaggcctgcaaacgtt                                              <  1:1527922‑M1/100‑35 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCcgttcctgcagtgtctgccaggcctgcaaacgttgt                                            <  1:1345606‑M1/100‑37 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgn                                            <  1:941176/100‑2 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  2:2226115/1‑100 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           <  1:1388063/100‑1 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCcgttcctgcagtgtctgccaggcctgcaaacgttgttt                                          <  1:1022472‑M1/100‑39 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          >  2:573733/1‑100 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          >  1:2250551/1‑100 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          >  1:2335123/1‑100 (MQ=255)
                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         >  2:3888194/1‑100 (MQ=255)
                                                       tCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         <  1:3353868/100‑1 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcgg                                         <  1:237364/99‑1 (MQ=255)
                                                         aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcggtt                                       >  2:3724514/1‑98 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                      <  1:2528655/100‑1 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGAGGCAc                                      >  2:4127156/1‑100 (MQ=255)
                                                                gCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCcgttcctgcagtgtctgccaggcctgcaaacgttgtttttttctcgcc                                <  2:957499‑M1/100‑49 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  1:1505292/100‑1 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  2:3172975/100‑1 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAAAAAAAAAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:3202087/99‑1 (MQ=255)
                                                                    aaaGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                            <  1:3355608/100‑1 (MQ=255)
                                                                        tttcgatctcccGGAACAGTTGGCGAGAAAAAATCAACGTTTGCAGGCCGGGCAGACACTGCATGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        <  2:4094432/91‑1 (MQ=255)
                                                                        ttccgatctcccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        <  2:4094433/91‑1 (MQ=255)
                                                                        ttGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        >  1:1409661/1‑100 (MQ=255)
                                                                        ttGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAc                        >  2:4053058/1‑100 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                       <  2:4243983/100‑1 (MQ=255)
                                                                          ggcATTGGCCTGAACAGATGGCGAGAAAAAAACAACGTTGGCAGGGCGGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAGAAc                        <  1:930510/96‑1 (MQ=255)
                                                                          gACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                      >  1:2139338/1‑100 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAATGTTTGCAGGCCTGGAAGACACTGTAGGAACGGAAATCCACGGCGGCACTGCTTGCAGAAAACCGc                     >  1:1560458/1‑100 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                     <  2:3632139/100‑1 (MQ=255)
                                                                             aTTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                   >  1:3918735/1‑100 (MQ=255)
                                                                              ttGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                  >  2:1383088/1‑100 (MQ=255)
                                                                                gCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa                >  2:1595845/1‑100 (MQ=255)
                                                                                 cccTGACCAGTAGGCGAGAAAAAAACAACGCTTGCAGGCCTGGCAGACCCTGCAAGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcagatcggaa               >  1:4094402/1‑91 (MQ=255)
                                                                                 cccTGAACAGTTGGCGCGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCGATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               >  1:3503863/1‑100 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcagatcggaa               >  1:4094403/1‑91 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               >  1:2843877/1‑100 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  2:2406868/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              <  1:1213893/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              >  2:1993204/1‑100 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGACTCGAt               >  1:572169/1‑99 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             >  1:1794499/1‑100 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa             <  2:2600816/100‑1 (MQ=255)
                                                                                    tGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg            >  1:2212024/1‑100 (MQ=255)
                                                                                     naaCAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           >  1:4254086/2‑100 (MQ=255)
                                                                                     gAACAGTTGGAGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGa           <  2:54315/100‑1 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      >  1:3997720/1‑100 (MQ=255)
                                                                                            tGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    <  1:4308498/100‑1 (MQ=255)
                                                                                             ggCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   >  2:431624/1‑100 (MQ=255)
                                                                                             ggCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   >  2:1728791/1‑100 (MQ=255)
                                                                                             ggCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGAGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatg   <  1:2072799/100‑1 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: