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.92.8% 494.0 / 26.8 125coding (308/444 nt)fliJflagellar protein
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base . (69/47);  total (73/52)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 141 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGATGG  >  NC_000913/2018154‑2018348
                                                                                                   |                                                                                               
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  2:3129290/1‑100 (MQ=255)
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAAAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                                 >  1:516976/1‑100 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                                >  2:3777499/1‑100 (MQ=255)
  cagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaa                                                                                               >  2:1241784/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                              <  2:3308730/100‑1 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaa                                                                                               >  2:838934/1‑99 (MQ=255)
     tttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaa                                                                                            <  2:790895/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  2:805390/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaa                                                                                           <  1:1006728/100‑1 (MQ=255)
      ttATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGAAGAAAAGaaaaaaa                                                                                           <  1:3510449/100‑1 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                        >  2:1784285/1‑97 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  1:2916101/1‑96 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACCTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                       >  1:2835220/1‑96 (MQ=255)
            aGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGt                                                                                    >  1:226836/1‑100 (MQ=255)
              aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGttt                                                                                  <  2:151476/100‑1 (MQ=255)
               cGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTg                                                                                 >  2:612162/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                               >  1:2176157/1‑100 (MQ=255)
                  tGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAg                                                                              >  2:524815/1‑100 (MQ=255)
                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAgg                                                                             <  1:3686479/100‑1 (MQ=255)
                    gAAAATGTCATTACTTAGCATCGCCAGCAACTTAATCAGTGGACGCATAAAGTTGACAGTGCCCTGAACAGTTGGCGAGCAAAAAGCAACGTTTGCa                                                                               >  1:1617855/1‑97 (MQ=255)
                    gAAAAAGCCAATATTCAGCATCGCCAGCAACTTAATCAGTGGACACAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                            >  1:3504169/1‑100 (MQ=255)
                       aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTg                                                                         <  2:910354/100‑1 (MQ=255)
                       aaaGCCATAAATCAGCATCGCCAGCAACATAATCAGTGGAAAAAGAAAGTTGACATTGCCCTGAAAAAAAAAAAAAAAAAAAACAACGTTTGCAGGCCTg                                                                         <  1:3694946/100‑1 (MQ=255)
                        aaGCCATTACCCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTgg                                                                        >  1:387134/1‑100 (MQ=255)
                         agcaataaatCAGAATCACCAACAACTTAAACAGTAAAAAAAGAAAGTTGACATTACCCTTAAAAAAAAAAAAAAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  1:1368691/91‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       >  2:1849973/1‑100 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       >  1:3034046/1‑100 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                       <  2:396884/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACACTTGGAGAGAACAACACCACGTATGAAGGCCtcca                                                                       >  2:1640459/1‑97 (MQ=255)
                            cATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGa                                                                    >  2:2239547/1‑100 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                   <  2:1317667/100‑1 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                                 <  1:1845787/100‑1 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCGCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  1:778556/1‑100 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                                >  1:398726/1‑100 (MQ=255)
                                 cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                               >  2:2046729/1‑100 (MQ=255)
                                  tCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  2:1935245/100‑1 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                             >  1:803578/1‑100 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAAAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            >  1:2885326/1‑100 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            >  2:3445528/1‑100 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            >  1:1961243/1‑100 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:1304818/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTCACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACCTGTGCAGCCTCGGGAGacac                                                                 >  1:2885185/1‑95 (MQ=255)
                                    aGCATCGCCAGCAACCTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                            <  2:3463127/100‑1 (MQ=255)
                                    aGCATCGCCAGCAAATTAATCAGGGGAGGCAGAAAGTTGACAGTGCCCGGAACAGATGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGc                                                              <  1:2962979/98‑1 (MQ=255)
                                       aTCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAACAACAACGTTTGCAGGCCTGGCAGACACTGCAGGaa                                                         >  1:1714390/1‑100 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                        >  1:315148/1‑99 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       <  2:2362753/100‑1 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       <  2:1580047/100‑1 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgatcgg                                                       >  1:821770/1‑95 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                       >  1:1886630/1‑99 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     <  2:2027407/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                     >  2:1780514/1‑100 (MQ=255)
                                                 aCTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                               >  2:3094021/1‑100 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              <  1:52305/100‑1 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:2194134/1‑100 (MQ=255)
                                                   ttAATCGGAGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACACCGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                              >  1:2661137/1‑99 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  1:233336/1‑100 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  2:2403528/1‑100 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            <  1:10435/100‑1 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  1:1901501/1‑100 (MQ=255)
                                                    tAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcgg                                            >  1:1543244/1‑100 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggc                                           >  2:3031655/1‑100 (MQ=255)
                                                      aTCAGTGGACGCAGATAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          <  2:3753809/100‑1 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                          <  1:1117943/100‑1 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        >  1:1797057/1‑100 (MQ=255)
                                                        cAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACggcggc                                        <  1:2760699/100‑1 (MQ=255)
                                                         aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       >  2:335568/1‑100 (MQ=255)
                                                         aGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCa                                       >  1:2102736/1‑100 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                      >  1:677234/1‑100 (MQ=255)
                                                            ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                    >  1:2822105/1‑100 (MQ=255)
                                                            ggACGAAGAAAGTTGACATTGCCCTGAACAGTTGGCGAAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                    <  1:510464/100‑1 (MQ=255)
                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                   <  2:654861/100‑1 (MQ=255)
                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                   >  2:2888665/1‑100 (MQ=255)
                                                              aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                  <  2:3451784/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  2:2439970/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  2:363427/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               >  1:1470467/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCGGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                               <  1:2246711/100‑1 (MQ=255)
                                                                  aGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCa                              <  2:2864155/100‑1 (MQ=255)
                                                                   gAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                             >  2:3201174/1‑100 (MQ=255)
                                                                    aaaGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGa                            <  2:2182004/100‑1 (MQ=255)
                                                                       gTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  1:955990/100‑1 (MQ=255)
                                                                       gTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                          <  1:3288937/99‑1 (MQ=255)
                                                                       attcaCATTGCCCTGAACAGCTGGCGATAAAAAAACATCGTTTGCAGGCGTGGCAGACACTGCAGGAACGCCAAGCCACGGCCGCACTGCTTGCCGa                            >  1:1763246/5‑97 (MQ=255)
                                                                        ttGACAGTGCCCTGAACAGTTGGCGATAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaaa                         <  1:3282356/99‑1 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                     <  1:321865/100‑1 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                    >  2:1649862/1‑100 (MQ=255)
                                                                               tGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                 <  1:2893405/100‑1 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  2:1883210/100‑1 (MQ=255)
                                                                                 cccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAt               <  2:12923/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc              >  2:651444/1‑100 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCATGCAGAAAACCGCCTCGATc              <  1:2326856/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg                 <  1:3537315/96‑1 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        <  2:2848148/100‑1 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa        >  1:583377/1‑100 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg      >  2:3684747/1‑100 (MQ=255)
                                                                                            tGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    <  2:567155/100‑1 (MQ=255)
                                                                                             ggCGAGAAAAAAACAACGTTTGCAGGCCTGGCATACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat    <  1:1999668/99‑1 (MQ=255)
                                                                                             ggCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCATAAAACCGCCTCGATCAGAAAAAgatg   >  1:2585922/1‑100 (MQ=255)
                                                                                              gCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgatgg  <  2:2093124/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: