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.96.8% 403.8 / 5.2 95coding (308/444 nt)fliJflagellar protein
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base . (39/53);  total (40/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGAT  >  NC_000913/2018154‑2018346
                                                                                                   |                                                                                             
agcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCgaga                                                                                               >  2:304201/1‑100 (MQ=255)
 gcagTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaa                                                                                              >  1:2897885/1‑100 (MQ=255)
   agTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaa                                                                                            >  2:270368/1‑100 (MQ=255)
        agccAGGCGGTGGAAAAAGCCATTACTCAGCATCGCAAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTTCCCTGAACAGTTGGCGAGaaaaaaa                                                                                         <  2:290406/96‑1 (MQ=255)
        aTCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaaca                                                                                       >  2:1761938/1‑98 (MQ=255)
         tCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                      >  2:1361481/1‑97 (MQ=255)
         nccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaa                                                                                      >  2:1499478/2‑97 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGaaaaaaacaac                                                                                     <  1:788282/100‑5 (MQ=255)
          ccAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAAAAGTTGGCGAGacaaaaacaac                                                                                     >  2:127126/1‑96 (MQ=255)
           cAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACg                                                                                   <  1:1615594/100‑1 (MQ=255)
            aGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGt                                                                                  >  1:1597787/1‑100 (MQ=255)
             gACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGtt                                                                                 <  2:54398/100‑1 (MQ=255)
              aCGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGttt                                                                                <  2:76404/100‑1 (MQ=255)
                tcTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCACTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGc                                                                              <  2:2150075/99‑1 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                             >  1:2809305/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                             >  1:2150104/1‑100 (MQ=255)
                 cTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCa                                                                             <  1:821533/100‑1 (MQ=255)
                  tGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAg                                                                            >  1:2302463/1‑100 (MQ=255)
                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAgg                                                                           >  1:907597/1‑100 (MQ=255)
                   ggAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAgg                                                                           <  1:270069/100‑1 (MQ=255)
                    gAAAAAGCCATTGCGCAGCATCGGCAGCAACTTGATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                          <  1:2653345/100‑1 (MQ=255)
                    gAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGc                                                                          <  2:1505388/100‑1 (MQ=255)
                       aaaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTg                                                                       <  1:1649050/100‑1 (MQ=255)
                        aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAATACAACGTTTGCAGGCCTgg                                                                      <  1:1160203/100‑1 (MQ=255)
                        aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTgg                                                                      <  2:1942121/100‑1 (MQ=255)
                        aaGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTgg                                                                      <  1:2657296/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                     <  1:1794527/100‑1 (MQ=255)
                         aGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGc                                                                     <  1:2596781/100‑1 (MQ=255)
                          gcaaTTACTCAGCATCGCTAGCAACTTAATCAGTGGACGCAGAAGGTTGGTGTTGCCCTGAACAGTTGGCGAGAAGAAAACAACGTTTGCAGGCCTGGc                                                                     <  2:1473592/96‑1 (MQ=255)
                           ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAg                                                                   <  2:880236/100‑1 (MQ=255)
                           ccATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAg                                                                   <  1:2155961/100‑1 (MQ=255)
                             nttACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                 >  2:2410866/2‑100 (MQ=255)
                             aTTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                 <  1:283251/100‑1 (MQ=255)
                              ttACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGac                                                                 <  1:696004/100‑1 (MQ=255)
                               tACTCAGCATTGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                               <  2:2010781/100‑1 (MQ=255)
                               tACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAAacac                                                               <  2:2035445/100‑1 (MQ=255)
                                aCTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACt                                                              >  2:1383769/1‑100 (MQ=255)
                                 cTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTg                                                             >  1:1273031/1‑100 (MQ=255)
                                   cAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCa                                                           <  1:162059/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAg                                                          <  1:2816608/100‑1 (MQ=255)
                                    aGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAAAAAAAAAGAAAAAAACAACGTTTGCAGGCCTGGCAGacac                                                               <  2:520961/95‑1 (MQ=255)
                                     gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                         >  1:2550749/1‑100 (MQ=255)
                                     gCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAgg                                                         >  2:2595987/1‑100 (MQ=255)
                                        tCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAAc                                                      <  1:304208/100‑1 (MQ=255)
                                         cGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACg                                                     >  2:1444611/1‑100 (MQ=255)
                                          gCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACgg                                                    >  1:1495134/1‑100 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                   <  1:1041842/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                   <  1:362820/100‑1 (MQ=255)
                                           ccAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGc                                                   <  1:1663935/100‑1 (MQ=255)
                                              gCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGGGAAAAGACACACGTTTGCAGGACTGGCAGAAAACGCAGGAACGGc                                                   >  2:2418232/1‑97 (MQ=255)
                                                 atctAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCa                                             <  2:1797538/97‑1 (MQ=255)
                                                  cTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAc                                            >  2:1991857/1‑100 (MQ=255)
                                                   ttAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcg                                           <  1:2359270/100‑1 (MQ=255)
                                                     aaTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCaaga                                          >  1:1797558/1‑96 (MQ=255)
                                                      aTCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCAcggcg                                        <  2:84666/100‑1 (MQ=255)
                                                          gTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCAc                                    <  2:2606873/100‑1 (MQ=255)
                                                            ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                  >  2:186777/1‑100 (MQ=255)
                                                            ggACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTg                                  >  1:3017823/1‑100 (MQ=255)
                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                 <  1:1383782/100‑1 (MQ=255)
                                                             gACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGc                                 >  2:1799382/1‑100 (MQ=255)
                                                              aCGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCt                                <  1:279473/100‑1 (MQ=255)
                                                               cGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCtt                               >  2:3053488/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                             <  2:2809266/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                             >  1:2748992/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                             >  2:2251299/1‑100 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                             <  1:1224823/100‑1 (MQ=255)
                                                                 cAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGc                             <  1:1188760/100‑1 (MQ=255)
                                                                   gAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                           >  1:2672140/1‑100 (MQ=255)
                                                                   gAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                           <  2:333941/100‑1 (MQ=255)
                                                                   gAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAg                           <  1:2881053/100‑1 (MQ=255)
                                                                      aGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGaaa                        <  1:2989535/100‑1 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCCATCCACGGCGGCACTGCTTGCAGAAAAcc                     >  2:485537/1‑100 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                     >  2:1335393/1‑100 (MQ=255)
                                                                         tGACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTCGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAAcc                     <  1:1055836/100‑1 (MQ=255)
                                                                          gACATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCg                    >  1:270903/1‑100 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                   >  1:191339/1‑100 (MQ=255)
                                                                           aCATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGc                   >  1:2328187/1‑100 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGGCACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                  >  1:2998073/1‑100 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                  <  1:76548/100‑1 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                  <  1:208440/100‑1 (MQ=255)
                                                                            cATTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGcc                  <  1:2136705/100‑1 (MQ=255)
                                                                             aTTGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCt                 >  1:1068403/1‑100 (MQ=255)
                                                                              ttGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTc                <  2:770127/100‑1 (MQ=255)
                                                                               tGCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCg               >  1:2175128/1‑100 (MQ=255)
                                                                                gCCCTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGa              >  1:677500/1‑100 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc            <  1:1486205/100‑1 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc            >  2:49800/1‑100 (MQ=255)
                                                                                  ccTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATc            <  1:2074831/100‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAACCCGCCTCGAt             <  1:39286/98‑1 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa           >  1:1442300/1‑100 (MQ=255)
                                                                                   cTGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCa           >  1:2611871/1‑100 (MQ=255)
                                                                                    tGAACAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAg          <  1:1619557/100‑1 (MQ=255)
                                                                                     gAACAGTTGGCGAGAAAAAAGCCACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGAGCAGa         >  2:1779894/1‑100 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa      <  2:2875020/100‑1 (MQ=255)
                                                                                        cAGTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGaaaa      >  1:1928544/1‑100 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg    <  1:2870374/100‑1 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg    >  1:682871/1‑100 (MQ=255)
                                                                                          gTTGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAg    <  1:1499494/100‑1 (MQ=255)
                                                                                            tGGCGAGAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAga   >  2:2668685/1‑99 (MQ=255)
                                                                                            tGGCGAAAAAAAAACAACGTTTGCAGGACTGGCAGACACTGCAGGAACAGAAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAgat  <  2:186532/100‑1 (MQ=255)
                                                                                                   |                                                                                             
AGCAGTTTATCCAGACGCTGGAAAAAGCCATTACTCAGCATCGCCAGCAACTTAATCAGTGGACGCAGAAAGTTGACATTGCCCTGAACAGTTGGCGAGAAAAAAAACAACGTTTGCAGGCCTGGCAGACACTGCAGGAACGGCAATCCACGGCGGCACTGCTTGCAGAAAACCGCCTCGATCAGAAAAAGAT  >  NC_000913/2018154‑2018346

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: