Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,134,118 A→C Q106P (CAG→CCG)  flgF → flagellar component of cell‑proximal portion of basal‑body rod

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,134,1180AC98.7% 529.9 / ‑4.4 153Q106P (CAG→CCG) flgFflagellar component of cell‑proximal portion of basal‑body rod
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (62/89);  total (63/90)
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 153 total aligned reads displayed.

ACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTG  >  NC_000913/1134021‑1134213
                                                                                                 |                                                                                               
acCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCgg                                                                                               <  2:2591276/100‑1 (MQ=255)
acCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCgg                                                                                               <  1:3878483/100‑1 (MQ=255)
acCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCgg                                                                                               <  1:646488/100‑1 (MQ=255)
acCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCgg                                                                                               <  2:1937492/100‑1 (MQ=255)
 ccTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGAGGGCTGGCTGGCCGTGCAGACCGCTGACGGCATCGGAGGGTATACGCGTACTGGCAGCATTCCGGt                                                                                              >  2:714563/1‑100 (MQ=255)
  cTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGtt                                                                                             >  2:1994288/1‑100 (MQ=255)
   tCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTg                                                                                            <  1:1836496/100‑1 (MQ=255)
     gcgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTTATGGCAGCATTCCGGTTGAt                                                                                          <  1:3179010/100‑1 (MQ=255)
      ngcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATc                                                                                         >  1:3249039/2‑100 (MQ=255)
      cgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATc                                                                                         >  1:2923479/1‑100 (MQ=255)
      cgcCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGATATACGCGTAATGGCAGCATTCCGGTTGATc                                                                                         <  1:3277175/100‑1 (MQ=255)
         ccGCTGGACGTCGCGTGGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCa                                                                                      <  1:1920388/100‑1 (MQ=255)
          cgcgGGACGCCGCGTTGCAGCAGGATGGCTGGCTGGCTGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCAc                                                                                     <  2:2072400/100‑1 (MQ=255)
          cGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCAc                                                                                     >  1:396901/1‑100 (MQ=255)
           gCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCAcc                                                                                    <  1:2063197/100‑1 (MQ=255)
            cTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGTAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCg                                                                                   <  1:1077348/100‑1 (MQ=255)
            cTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGTAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCg                                                                                   <  1:762656/100‑1 (MQ=255)
            cTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTCTACGCGTAATGGCAGCATTCCGGTTGATCCCACCg                                                                                   >  2:1372335/1‑100 (MQ=255)
              ggACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCggg                                                                                 <  2:263717/100‑1 (MQ=255)
              ggACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCggg                                                                                 >  2:2647284/1‑100 (MQ=255)
                  gTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAAc                                                                             <  1:1689892/100‑1 (MQ=255)
                  gTCGCGTTGCAGCAGGAAGGCTGGCTGGGCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGGTCCCACCGGGCAAc                                                                             >  2:644818/1‑100 (MQ=255)
                   tCGCGTTGCAGTAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACt                                                                            <  2:3198499/100‑1 (MQ=255)
                   tCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACt                                                                            >  2:896555/1‑100 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                         <  1:2656102/100‑1 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                         <  2:3977321/100‑1 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGAc                                                                         <  2:3975154/100‑1 (MQ=255)
                      cgTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGAGAAGCGTATACGCGTAATAGCAGCATTCCGGCTGATCCCACCGGGCAACTGAc                                                                         >  2:2311457/1‑100 (MQ=255)
                       gTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACa                                                                        >  2:3559166/1‑100 (MQ=255)
                       gTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACa                                                                        <  1:1434388/100‑1 (MQ=255)
                        ttGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACaa                                                                       <  1:3212456/100‑1 (MQ=255)
                         tGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAAt                                                                      <  1:3477935/100‑1 (MQ=255)
                         tGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAAt                                                                      <  2:3017274/100‑1 (MQ=255)
                            agcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCa                                                                   <  2:2898724/100‑1 (MQ=255)
                            agcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCa                                                                   >  1:571942/1‑100 (MQ=255)
                             gcagGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAg                                                                  <  1:3540094/100‑1 (MQ=255)
                             gcagGATGGCTGGCTCGCCCTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAg                                                                  <  2:1005297/100‑1 (MQ=255)
                                ggATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAgggg                                                               <  2:164127/100‑1 (MQ=255)
                                 gATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGc                                                              <  2:531751/100‑1 (MQ=255)
                                 gATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGc                                                              <  1:1233250/100‑1 (MQ=255)
                                    tgctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATc                                                           <  1:3318264/99‑1 (MQ=255)
                                     gctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATcc                                                          <  1:3689612/100‑1 (MQ=255)
                                     gctggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATca                                                          >  2:3263509/1‑99 (MQ=255)
                                       tggctggcCGTGCAGACCGCTGAGGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCAACCgg                                                        <  2:3396995/100‑1 (MQ=255)
                                       tggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCgg                                                        >  1:3432160/1‑100 (MQ=255)
                                        ggctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGt                                                       >  2:1116784/1‑100 (MQ=255)
                                         gctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTg                                                      <  2:625415/100‑1 (MQ=255)
                                          ctggcCGTGCAGACCGCTGACGGCAGCGAATGGTATACGCGTAATAGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGa                                                     >  1:3543566/1‑100 (MQ=255)
                                          ctggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGa                                                     <  2:2498647/100‑1 (MQ=255)
                                           tggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGAt                                                    <  2:1483186/100‑1 (MQ=255)
                                           nggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGAt                                                    >  1:1646990/2‑100 (MQ=255)
                                            ggcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATa                                                   >  1:1489805/1‑100 (MQ=255)
                                             gcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGCCAATTCAGGGGCATCCGGTGATAg                                                  >  2:236775/1‑100 (MQ=255)
                                             gcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAg                                                  >  2:538356/1‑100 (MQ=255)
                                             gcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAg                                                  >  2:3118726/1‑100 (MQ=255)
                                             gcCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAg                                                  >  1:1576600/1‑100 (MQ=255)
                                               cGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGc                                                <  1:462384/100‑1 (MQ=255)
                                               cGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGc                                                <  2:656361/100‑1 (MQ=255)
                                               cGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGc                                                <  1:2379054/100‑1 (MQ=255)
                                                 tGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGa                                              >  1:182220/1‑100 (MQ=255)
                                                  gCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAg                                             <  1:925709/100‑1 (MQ=255)
                                                   cAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAgg                                            >  1:3309349/1‑100 (MQ=255)
                                                   cAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTAGATCCCACCGCGCAACTGACAATTCCGGGGCATCCGGTCACAGGCCAgg                                            >  2:3940582/1‑100 (MQ=255)
                                                       ccGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTgg                                        >  1:3058989/1‑100 (MQ=255)
                                                        cGCTGACGGCAGCGAAGGGTATACGCATAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTggg                                       <  2:1872644/100‑1 (MQ=255)
                                                         gCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGc                                      >  2:728821/1‑100 (MQ=255)
                                                         gCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGc                                      >  1:2017718/1‑100 (MQ=255)
                                                          tgcaCGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGcn                                     <  2:1478460/97‑2 (MQ=255)
                                                          cTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGcc                                     <  1:3496849/100‑1 (MQ=255)
                                                           tGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCa                                    >  2:627325/1‑100 (MQ=255)
                                                             aCGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAt                                  <  2:138963/100‑1 (MQ=255)
                                                              cGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAAtt                                 >  1:2812859/1‑100 (MQ=255)
                                                                gCAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGc                               >  2:2345373/1‑100 (MQ=255)
                                                                 cAGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCt                              >  1:690529/1‑100 (MQ=255)
                                                                 cAGCGAAGGGTATACGCGGAATGGCAGCATTCCAGTTGATCCCACCGGGCAACTGACAATTCAGAGGCATCCGGTGATAGGCGCGGCTGGGCCAATTGCt                              >  2:2835996/1‑100 (MQ=255)
                                                                  aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtg                             <  2:1623299/100‑1 (MQ=255)
                                                                  aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtg                             <  2:2068533/100‑1 (MQ=255)
                                                                  aGCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtg                             <  1:2778524/100‑1 (MQ=255)
                                                                   gCGAAGGGTATACGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCtgt                            >  1:2629024/1‑100 (MQ=255)
                                                                     gAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGc                          >  2:3713477/1‑100 (MQ=255)
                                                                           tataCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAg                    <  2:607339/100‑1 (MQ=255)
                                                                            ataCGCGTAATGGCAGCATTCCGGTTGCTCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgg                   <  1:745804/100‑1 (MQ=255)
                                                                             taCGCGTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAggg                  >  2:2456905/1‑100 (MQ=255)
                                                                              aCGCGTAAGGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAgggg                 <  2:3347877/100‑1 (MQ=255)
                                                                                gcgTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCg               <  2:2982802/100‑1 (MQ=255)
                                                                                 cgTAATGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCgg              <  2:391293/100‑1 (MQ=255)
                                                                                     aTGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAAt          <  1:202210/100‑1 (MQ=255)
                                                                                      tGGCAGCATTCCGGTTGATCCCACCGGGCAATTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGACATTTGCTGTGCCGGAAGGGGCGGAAATc         >  2:3805116/1‑100 (MQ=255)
                                                                                      tGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATc         <  1:61686/100‑1 (MQ=255)
                                                                                      tGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATc         >  1:1744138/1‑100 (MQ=255)
                                                                                      tGGCAGCATTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATc         <  1:1921718/100‑1 (MQ=255)
                                                                                       ggCAGCATTCCGGTTGATCCGACCGGGCAACTGACAATTCAGGGGCATTCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAAtga        >  2:3985828/1‑98 (MQ=255)
                                                                                        gcagcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCAc       <  1:3020431/100‑1 (MQ=255)
                                                                                        gcagcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCAc       <  1:3020466/100‑1 (MQ=255)
                                                                                        gcagcaTTCCGGGTGATCCCGCCGGGCGACTGGCCATTCAGGGGCATCGGGTGGTAGGCGTGGCTGGGCAAATGGCTGGGCCGAAAGGGGCGGAAATCAc       <  2:603711/100‑1 (MQ=255)
                                                                                        gatctaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCAc       <  2:2570891/95‑1 (MQ=255)
                                                                                         cagcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTCGGCCAATTGCTGTGACGGAGGGGGCGGAAATCACt      >  2:1777065/1‑100 (MQ=255)
                                                                                           gcaTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTAt    >  1:3205599/1‑100 (MQ=255)
                                                                                             aTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCAcagatc  >  1:2610740/1‑95 (MQ=255)
                                                                                             aTTCCGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATtg  >  1:2116573/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ACCTCGCGCCCGCTGGACGTCGCGTTGCAGCAGGATGGCTGGCTGGCCGTGCAGACCGCTGACGGCAGCGAAGGGTATACGCGTAATGGCAGCATTCAGGTTGATCCCACCGGGCAACTGACAATTCAGGGGCATCCGGTGATAGGCGAGGCTGGGCCAATTGCTGTGCCGGAAGGGGCGGAAATCACTATTG  >  NC_000913/1134021‑1134213

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: