Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,967,082 T→C K119R (AAA→AGA)  cheY ← chemotaxis regulator transmitting signal to flagellar motor component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,967,0820TC98.6% 284.5 / NA 72K119R (AAA→AGA) cheYchemotaxis regulator transmitting signal to flagellar motor component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/42);  minor base A (1/0);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGT  >  NC_000913/1966937‑1967227
                                                                                                                                                 |                                                                                                                                                 
tCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTtcc                                                                                                                                              >  2:2966/1‑151 (MQ=255)
tCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTtcc                                                                                                                                              >  1:264974/1‑151 (MQ=255)
  ccGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGCCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTtcct                                                                                                                                             <  1:1232420/150‑1 (MQ=255)
     cAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTtcctcc                                                                                                                                           >  1:316132/1‑149 (MQ=255)
       aaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCTTCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGAGTTTGTTGAGTCTTTCCTCCAGCg                                                                                                                                       <  2:703595/151‑1 (MQ=255)
       aaCTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCg                                                                                                                                       <  1:602630/151‑1 (MQ=255)
          tGTCGCGCAGCATACGCGTCAGGCTGCCGATGCTCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTATTTCCTCCAGCg                                                                                                                                       >  2:616866/1‑148 (MQ=255)
          tGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTcg                                                                                                                                    >  2:75025/1‑151 (MQ=255)
           gTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTc                                                                                                                                     <  2:954616/149‑1 (MQ=255)
             cgcgCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCgcgg                                                                                                                                 <  1:1164206/151‑1 (MQ=255)
                gcagcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGcggcgg                                                                                                                              <  1:478668/151‑1 (MQ=255)
                 cagcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGt                                                                                                                             >  1:27779/1‑151 (MQ=255)
                  agcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGt                                                                                                                             >  1:560346/1‑150 (MQ=255)
                   gcaTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTaa                                                                                                                           <  1:160571/151‑1 (MQ=255)
                     aTACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAAt                                                                                                                         >  1:631512/1‑151 (MQ=255)
                       aCGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATgg                                                                                                                       <  1:453172/151‑1 (MQ=255)
                       aCGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATgg                                                                                                                       <  2:540604/151‑1 (MQ=255)
                          cgTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCtt                                                                                                                    >  1:643754/1‑151 (MQ=255)
                            tCATGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTc                                                                                                                   >  1:1027757/1‑150 (MQ=255)
                             cAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTcac                                                                                                                 >  1:346484/1‑151 (MQ=255)
                                   gCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACcaca                                                                                                             >  1:887776/1‑149 (MQ=255)
                                   gCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATa                                                                                                           <  2:911597/151‑1 (MQ=255)
                                   gCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATa                                                                                                           <  2:704249/151‑1 (MQ=255)
                                         gcgcgcAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACt                                                                                                     <  1:709420/151‑1 (MQ=255)
                                               aaTGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGGCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGc                                                                                                  <  2:439459/148‑1 (MQ=255)
                                                     atCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTg                                                                                         <  2:939727/151‑1 (MQ=255)
                                                          cAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGcgccg                                                                                    <  1:78779/151‑1 (MQ=255)
                                                           aGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGcgccgc                                                                                   >  2:105313/1‑151 (MQ=255)
                                                           aGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGcgccgc                                                                                   >  1:395094/1‑151 (MQ=255)
                                                                 aTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAat                                                                             >  2:949640/1‑151 (MQ=255)
                                                                 aTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAat                                                                             <  2:239821/151‑1 (MQ=255)
                                                                 aTGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAat                                                                             <  2:233849/151‑1 (MQ=255)
                                                                   gCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatga                                                                           <  1:462825/151‑1 (MQ=255)
                                                                   gCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatg                                                                            >  1:288048/1‑150 (MQ=255)
                                                                      cGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAatgat                                                                          <  2:563820/149‑1 (MQ=255)
                                                                         cagcagGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtct                                                                     <  2:960382/151‑1 (MQ=255)
                                                                         cagcagGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTtct                                                                     >  2:256771/1‑151 (MQ=255)
                                                                           gcagGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTctct                                                                   <  2:867416/151‑1 (MQ=255)
                                                                             agGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCtcttt                                                                 <  1:761043/151‑1 (MQ=255)
                                                                                 ttgattgaTGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTttcttc                                                             >  1:237584/1‑151 (MQ=255)
                                                                                   gattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCg                                                            <  1:312423/150‑1 (MQ=255)
                                                                                   gattgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGc                                                           <  2:903108/151‑1 (MQ=255)
                                                                                    attgatGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCt                                                          <  2:672421/151‑1 (MQ=255)
                                                                                          ggTTGCATCATAGTCGCAACCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCCCCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCATGATGTTCTCTTTCTTCGCTTCTGCAg                                                   <  2:48364/151‑1 (MQ=255)
                                                                                           gTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                                   <  2:73404/151‑1 (MQ=255)
                                                                                            ttGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAg                                                   >  2:226907/1‑150 (MQ=255)
                                                                                              gCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCa                                                <  1:530579/151‑1 (MQ=255)
                                                                                              gCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCa                                                <  1:855252/151‑1 (MQ=255)
                                                                                                atcatAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCAcc                                              >  1:316268/1‑151 (MQ=255)
                                                                                                 tcatAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCAcc                                              >  1:394262/1‑150 (MQ=255)
                                                                                                  catAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCAt                                            <  2:316125/151‑1 (MQ=255)
                                                                                                  catAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCAt                                            <  2:504868/151‑1 (MQ=255)
                                                                                                   atAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCAtt                                           >  1:1049846/1‑151 (MQ=255)
                                                                                                      gTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAac                                        <  1:1245520/151‑1 (MQ=255)
                                                                                                        cGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAacac                                      >  1:225516/1‑151 (MQ=255)
                                                                                                        cGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAacac                                      <  1:998509/151‑1 (MQ=255)
                                                                                                        cGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAaca                                       <  2:849734/150‑1 (MQ=255)
                                                                                                        cGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGAAGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAact                                       <  2:397291/150‑2 (MQ=255)
                                                                                                           aTCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTgg                                   >  1:319334/1‑151 (MQ=255)
                                                                                                            tCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGc                                  >  1:652596/1‑151 (MQ=255)
                                                                                                               tCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAAt                               >  1:671431/1‑151 (MQ=255)
                                                                                                                     gCCCAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGAc                         <  1:899681/151‑1 (MQ=255)
                                                                                                                      cccAGTTTCTCAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGAc                         >  1:1216575/1‑150 (MQ=255)
                                                                                                                             tctcAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGTCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAAGGATGTTCTCTGTCTTCGCTTCTGCAGTCACCATTATCACTGGCAATGCCGACATcgcac                  <  2:105477/150‑3 (MQ=255)
                                                                                                                                cAAAGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCa                <  2:653175/149‑1 (MQ=255)
                                                                                                                                 aaaGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATcc             <  1:824885/151‑1 (MQ=255)
                                                                                                                                 aaaGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATcc             >  2:112461/1‑151 (MQ=255)
                                                                                                                                   aGATTTTGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGc           >  2:737123/1‑151 (MQ=255)
                                                                                                                                      ttttGTTGAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCCTGCGCCGCAGCAATGATGTTCTCTTTTTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACg        >  2:89751/1‑151 (MQ=255)
                                                                                                                                        ttgttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGa       <  2:536854/150‑1 (MQ=255)
                                                                                                                                          gttgAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAAt     <  2:349956/150‑1 (MQ=255)
                                                                                                                                             gAGTCTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGt  <  2:589261/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TCCCGCAAACTGTCGCGCAGCATACGCGTCAGGCTGCCGATGCGCGCAATGATATCGCCAGCTGAATGCTCGTCAGCAGGTTTGATTGATGGTTGCATCATAGTCGCATCCTCACATGCCCAGTTTCTCAAAGATTTTGTTGAGTTTTTCCTCCAGCGTCGCGGCGGTAAATGGCTTCACCACATAGCCACTGGCCCCCGCTTGCGCCGCAGCAATGATGTTCTCTTTCTTCGCTTCTGCAGTCACCATTAACACTGGCAATGCCGACATCGCGCCATCCGCACGAATTGT  >  NC_000913/1966937‑1967227

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: