Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,974,817 C→G 5.1% G170R (GGG→CGG)  cheA ← fused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,974,8170CG 5.1% 266.1 / 3.7 79G170R (GGG→CGG) cheAfused chemotactic sensory histidine kinase in two‑component regulatory system with CheB and CheY: sensory histidine kinase/signal sensing protein
Reads supporting (aligned to +/- strand):  ref base C (32/42);  new base G (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

ATCACAAAACAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGG  >  NC_000913/1974678‑1974965
                                                                                                                                           |                                                                                                                                                    
aTCACAACAGAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGGAATTTTGGCGCGAGCGAATCCGACCCGTTCACCACGTCTGTTTACGTTGTCAGATGTGCGTGTTCTTCGTCCAGCAGGGGGGCTTCCCGGGCCtgccg                                                                                                                                             >  8:319848/1‑145 (MQ=255)
aTCACAAAACAGAGTACCGCTGTGAGGTCATCGTCGGCGATGTCGCCCGGGAATGTTGCCGGGAGCGATTCCGCCCCGTTCCCGCCGTCAGTTAACGTTGTCAGATGTCCCTGTTCGTCGTCCAGGAGGTCGACGTCCCGGGccggcag                                                                                                                                             >  8:325398/1‑144 (MQ=255)
    caAAACAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGgcgc                                                                                                                                         >  8:122486/1‑149 (MQ=255)
        aCAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGaaa                                                                                                                                     <  1:244050/149‑1 (MQ=255)
        aCAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGaaa                                                                                                                                     <  1:343108/149‑1 (MQ=255)
        aCAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGaaa                                                                                                                                     <  2:135177/149‑1 (MQ=255)
        aCAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGa                                                                                                                                       >  7:248205/1‑147 (MQ=255)
         cAGAGTACCGCTGTGATGTCATCTTCGGCGATGGCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAg                                                                                                                                    <  5:326943/149‑1 (MQ=255)
           gagTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGa                                                                                                                                  <  3:367706/149‑1 (MQ=255)
           gagTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCAACACGTCAGTTAACGTTGTCAGATGTACCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGa                                                                                                                                  <  4:98160/149‑1 (MQ=255)
                ccGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAAt                                                                                                                              <  5:80837/148‑1 (MQ=255)
                ccGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAAt                                                                                                                              >  5:147785/1‑148 (MQ=255)
                ccGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCATCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAAtt                                                                                                                             <  3:93399/149‑1 (MQ=255)
                 cGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATtc                                                                                                                            <  7:398565/149‑1 (MQ=255)
                 cGCTGTGATGTCATCTTCGGCGATGGCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATtc                                                                                                                            <  6:51529/149‑1 (MQ=255)
                      tgATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATtcgtc                                                                                                                         <  3:44034/147‑1 (MQ=255)
                           tCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCgac                                                                                                                 <  6:185079/149‑1 (MQ=255)
                            cATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCga                                                                                                                  <  7:44892/148‑1 (MQ=255)
                            cATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCga                                                                                                                  >  8:44892/1‑148 (MQ=255)
                             aTCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCgact                                                                                                                >  5:320314/1‑149 (MQ=255)
                                  cGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCgactga                                                                                                              <  8:248205/146‑1 (MQ=255)
                                  cGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCgactga                                                                                                              <  1:26771/146‑1 (MQ=255)
                                  cGGCGATGTCGCCAGGTACTATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCATCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGGCGCGGCgactgact                                                                                                            <  4:105635/148‑1 (MQ=255)
                                          tCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTg                                                                                                   <  7:115945/149‑1 (MQ=255)
                                            gCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCATCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGc                                                                                                  <  3:100571/148‑1 (MQ=255)
                                             cccGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGgcg                                                                                                                                          >  2:90136/1‑107 (MQ=255)
                                             cccGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGgcg                                                                                                                                          <  1:90136/107‑1 (MQ=255)
                                             cccGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCt                                                                                                 >  4:382069/1‑148 (MQ=255)
                                              ccGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCATCCAGCAGGTCGACTTCCCCGGCCTTCAGTCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCa                                                                                               <  4:42633/149‑1 (MQ=255)
                                                 gTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCt                                                                                            >  4:257885/1‑149 (MQ=255)
                                                       ttGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTATGCTCATCTTGCGGt                                                                                      >  4:38618/1‑149 (MQ=255)
                                                              gagCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACttt                                                                               >  3:172120/1‑149 (MQ=255)
                                                               agCGAATCCGCCCCTTTCACCACGTCAGTTAACGGTGTCAGATGTCCCAGTTCTTCTGCCAGCTGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACtttt                                                                              <  5:32047/149‑1 (MQ=255)
                                                                gCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTg                                                                             <  4:369943/149‑1 (MQ=255)
                                                                 cGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTgg                                                                            >  4:243084/1‑149 (MQ=255)
                                                                 cGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTgg                                                                            >  4:248197/1‑149 (MQ=255)
                                                                     tCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAAc                                                                        <  6:349188/149‑1 (MQ=255)
                                                                        gCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACca                                                                      >  7:99428/1‑148 (MQ=255)
                                                                          cccTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACGCTGCTCATCTTGCGGTTCACTTTTGGCAACCacac                                                                   >  1:227537/1‑149 (MQ=255)
                                                                           ccTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACt                                                                  <  4:183372/149‑1 (MQ=255)
                                                                           ccTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACt                                                                  >  3:159975/1‑149 (MQ=255)
                                                                             tttCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACtt                                                                 >  4:42169/1‑148 (MQ=255)
                                                                             tttCACCACGTCAGTTAACGTTCACAGATGTCCCAGTTCTACTTCCAGCAGCTCGACTTCCCCGGCCTTCACGCGCGAAAGGATAATTCGTCGCGGCGCCTGACTGCGACACTGCTCATCTTGCGGTTCACTTTTGGCAACCACACtt                                                                 >  2:360535/1‑148 (MQ=255)
                                                                                 accacGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATc                                                             <  7:122486/148‑1 (MQ=255)
                                                                                  ccacGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCgg                                                           <  3:382069/149‑1 (MQ=255)
                                                                                  ccacGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCgg                                                           <  2:244302/149‑1 (MQ=255)
                                                                                     cGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGt                                                         >  5:175890/1‑148 (MQ=255)
                                                                                        cAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACt                                                     >  7:7770/1‑149 (MQ=255)
                                                                                           ttAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCg                                                  >  4:280233/1‑149 (MQ=255)
                                                                                              aCGTTGTCAGATGTCCCTGTTCTTCTTCCAGCAGGTCGACTTCCCGGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGa                                                <  7:319848/148‑1 (MQ=255)
                                                                                              aCGTTGTCAGATGTCCCTGTTCTTCTTCCAGCAGGTCGACTTCCCGGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGa                                                <  7:325398/148‑1 (MQ=255)
                                                                                              aCGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGa                                                <  5:29789/148‑1 (MQ=255)
                                                                                              aCGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGa                                                <  5:9775/148‑1 (MQ=255)
                                                                                              aCGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGAt                                               >  6:147182/1‑149 (MQ=255)
                                                                                                       gATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCCGATGGCGTTTCg                                      >  1:339697/1‑149 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  6:39293/1‑149 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  4:94488/1‑149 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  1:99691/1‑149 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  3:2325/1‑149 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                <  3:42169/149‑1 (MQ=255)
                                                                                                             ccAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  3:323402/1‑149 (MQ=255)
                                                                                                                gTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                >  6:208862/1‑146 (MQ=255)
                                                                                                                gTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTc                                <  5:208862/146‑1 (MQ=255)
                                                                                                                    ttcttcAAGCAGGTCGACTTACCAGGCCTTCAGGCGCGAAAGGATATTTCGTCGCGGCGACTGACTGCGCCTCTGCTCCTCTTGCGGTTCACTTTTGGGCACCACAGGTAATCGGGTCACTGAGGATGGCGTTTCGACGTTCGCTTCt                          >  6:296650/1‑148 (MQ=255)
                                                                                                                     tcttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        >  6:34219/1‑149 (MQ=255)
                                                                                                                      cttcCAGCAGGTCGACTTCCCCGGGCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        <  4:381564/148‑1 (MQ=255)
                                                                                                                      cttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        <  4:172120/148‑1 (MQ=255)
                                                                                                                      cttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        <  3:280233/148‑1 (MQ=255)
                                                                                                                      cttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        <  3:248197/148‑1 (MQ=255)
                                                                                                                      cttcCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTaa                        <  3:243084/148‑1 (MQ=255)
                                                                                                                         ccAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAAt                                                              <  6:117631/107‑1 (MQ=255)
                                                                                                                         ccAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAAt                                                              >  5:117631/1‑107 (MQ=255)
                                                                                                                         ccAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAAt                       <  8:7770/146‑1 (MQ=255)
                                                                                                                            gcagGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGt                 >  8:358329/1‑149 (MQ=255)
                                                                                                                             cagGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCTACTGACTGCGACTCTGATCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGtt                <  8:383328/149‑1 (MQ=255)
                                                                                                                                gTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGAc             <  1:28153/149‑1 (MQ=255)
                                                                                                                                    aCTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCaa         <  6:320314/149‑1 (MQ=255)
                                                                                                                                     cTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAg        <  5:147182/149‑1 (MQ=255)
                                                                                                                                     cTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAg        >  3:129747/1‑149 (MQ=255)
                                                                                                                                        ccccGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGcc     >  2:5260/1‑149 (MQ=255)
                                                                                                                                        ccccGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGcc     >  8:323383/1‑149 (MQ=255)
                                                                                                                                         cccGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCt    <  5:39293/149‑1 (MQ=255)
                                                                                                                                          ccGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTg   >  1:25700/1‑149 (MQ=255)
                                                                                                                                           cGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTgg  >  2:160711/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
ATCACAAAACAGAGTACCGCTGTGATGTCATCTTCGGCGATGTCGCCCGGTAATATTGCCGAGAGCGAATCCGCCCCTTTCACCACGTCAGTTAACGTTGTCAGATGTCCCAGTTCTTCTTCCAGCAGGTCGACTTCCCCGGCCTTCAGGCGCGAAAGGATAATTCGTCGCGGCGACTGACTGCGACTCTGCTCATCTTGCGGTTCACTTTTGGCAACCACACTTAATCGGGTCACTGCGGATGGCGTTTCGCCTTTCGCTTCTAATGCCAGTTGACGCAAGGCCTGG  >  NC_000913/1974678‑1974965

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: