Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,587,013 A→G 8.6% Q428R (CAG→CGG) ‡ narQ → sensory histidine kinase in two‑component regulatory system with NarP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,587,0130AG 8.6% 295.0 / 15.3 105Q428R (CAG→CGG) ‡narQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base A (38/57);  new base G (6/3);  total (44/61)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.84e-01

Only 100 of 114 total aligned reads displayed.

TTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCG  >  NC_000913/2586864‑2587161
                                                                                                                                                     |                                                                                                                                                    
ttacttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgca                                                                                                                                                      <  2:1107490/150‑1 (MQ=255)
 tccttacGTCTCCAGTTGACGTTACTGAAGCGGTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATCATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTATCCGCCTGACGCTgcag                                                                                                                                                     <  2:1622832/148‑1 (MQ=255)
 tccttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcag                                                                                                                                                     <  1:1156119/148‑1 (MQ=255)
 tacttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcag                                                                                                                                                     >  1:418940/1‑150 (MQ=255)
   cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcagca                                                                                                                                                   >  2:1164297/1‑150 (MQ=255)
   cttacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTgcagca                                                                                                                                                   >  1:523326/1‑150 (MQ=255)
     tacGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGAGTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAgg                                                                                                                                                 >  1:138618/1‑150 (MQ=255)
        gcaTCCAGTTGACGTTCCTGACGCGTTCGATCCCGGAAGCTCACGCCACCGCACAAAGGATCATGGCCGATTTTTCCCCGGCGTTGAATGAGGCTTCGCGGCAGGTACGCGAGCTGTTGACTACTTTCCGCCTGACGCGGCAGCAGGCgg                                                                                                                                              <  1:1158875/148‑1 (MQ=255)
        gTATCCAGTTGACGTTACTGCAGCGTTCGATACCCGAAGATACCGCCACCGCACAAAGTATCCTGGCCGATTTTTCCCAGGCGTTGAATGATGCATATCGGCAGTTACGCGAGCTGTGGACTACTTTCCGCCTGACGCTGCAGCAGGCgg                                                                                                                                              <  2:673205/150‑1 (MQ=255)
        gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGCTGCAGCGTCAGGCgg                                                                                                                                              <  1:1045772/150‑1 (MQ=255)
        gTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCgg                                                                                                                                              <  2:1484484/150‑1 (MQ=255)
          aTCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGCTGCAGCGTCAGGCGGaa                                                                                                                                            <  1:863873/150‑2 (MQ=255)
           tCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtc                                                                                                                                           <  2:1184012/150‑1 (MQ=255)
             cAGTTGACGTTACTGAAGCGTTCGATACTGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCGGTTGACTACTTTCCGCCTGCTGCAGCGTCAGGCGGaaagt                                                                                                                                         <  2:1667940/150‑5 (MQ=255)
             cAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtctc                                                                                                                                         >  2:377555/1‑150 (MQ=255)
             cAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTGGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGAtctc                                                                                                                                         <  1:42438/150‑1 (MQ=255)
              aGTTGACGTTACTGAAGCGTGCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTcc                                                                                                                                        <  2:1510311/150‑1 (MQ=255)
                    cgaaaCTGACGCGTTCGCTACCGGACGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCg                                                                                                                                  <  2:147269/146‑1 (MQ=255)
                    cGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCg                                                                                                                                  <  2:269287/150‑1 (MQ=255)
                     gTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGc                                                                                                                                 >  1:1830519/1‑150 (MQ=255)
                         cTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTg                                                                                                                             >  2:534912/1‑150 (MQ=255)
                          tGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGa                                                                                                                            >  2:777128/1‑150 (MQ=255)
                           gAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAg                                                                                                                           >  2:1436303/1‑150 (MQ=255)
                            caGCGTTCGATACCGGAAGATGACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCGTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAgg                                                                                                                          <  1:205333/149‑1 (MQ=255)
                             aGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATGTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAggg                                                                                                                         <  2:897708/150‑1 (MQ=255)
                              gCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATTGCCGATTTTTCCCAGGCGTTGATTGATGCTTATCGGCAGTAACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGCCGGATCTCCCCTCCGCATTGAGGGa                                                                                                                        >  1:134427/1‑150 (MQ=255)
                              gCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGa                                                                                                                        >  2:1767317/1‑150 (MQ=255)
                               cGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGaa                                                                                                                       >  2:138268/1‑150 (MQ=255)
                                gTTCGATACAGGAAGATAACGCCACCGCACAAAGTATCATCGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGAGAGCTGTTGAATACTTTCCGCCTGACCCTGCAGCAGGAGGATCTCCCATCCGCATTTAGGGaaa                                                                                                                      >  1:1588735/1‑150 (MQ=255)
                                 ttCGATACCGGAAGATAAAGCCACCTCACAAAGTATCTTGGCCGATTTTTGCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAAt                                                                                                                     <  1:1306737/150‑1 (MQ=255)
                                  tCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCATATCAGCACTTACCCGAGCTGTAGACAACTTTCCGCCTGACGCTGCAGCTGGAGGCTCTCCCCTCCGCATTGAGGGAAATg                                                                                                                    >  1:305279/1‑150 (MQ=255)
                                   cGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGc                                                                                                                   >  2:1490259/1‑150 (MQ=255)
                                     aTACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTg                                                                                                                 <  1:1770916/150‑1 (MQ=255)
                                      tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGCTGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCtcg                                                                                                                >  2:164800/1‑148 (MQ=255)
                                      tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTgg                                                                                                                <  1:219979/150‑1 (MQ=255)
                                      tACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTgg                                                                                                                <  1:1003471/150‑1 (MQ=255)
                                       aCCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGAGTTTTCCCAGGCGTTGAATGATGCATATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGa                                                                                                               <  2:143507/150‑1 (MQ=255)
                                        ccGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGAt                                                                                                              >  2:1607713/1‑150 (MQ=255)
                                             aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                         >  1:1748618/1‑150 (MQ=255)
                                             aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                         <  1:890805/150‑1 (MQ=255)
                                             aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                         <  2:1003009/150‑1 (MQ=255)
                                             aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                         >  2:392198/1‑150 (MQ=255)
                                                   cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaaa                                                                                                   >  1:1498081/1‑150 (MQ=255)
                                                   cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaaa                                                                                                   <  2:1573269/150‑1 (MQ=255)
                                                    gCCACCGCCCAAAGGATCATGGCCGATTTTTCCCCGGCGGTGAATGATGCTTAGCGGCAGTTACGCGAGCTGTTGCCTTCTTTGCGCCTGACGCTGCAGCAGGCGGATCTCGCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAAt                                                                                                  <  2:1134403/150‑1 (MQ=255)
                                                     ccACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATc                                                                                                 <  1:1471246/150‑1 (MQ=255)
                                                          gCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAAcc                                                                                            <  2:218942/150‑1 (MQ=255)
                                                                       tGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGa                                                                               <  2:789609/150‑1 (MQ=255)
                                                                        ggCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCACACTGAc                                                                              >  1:289114/1‑150 (MQ=255)
                                                                         gCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAcc                                                                             >  1:230471/1‑150 (MQ=255)
                                                                         gCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAcc                                                                             >  1:871514/1‑150 (MQ=255)
                                                                          ccGATTTTTACCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGAccc                                                                            <  1:123346/150‑1 (MQ=255)
                                                                           ccaTTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTCACGCTGCAGCAGGCGGATCTCCACACCGCATTGAGGGAAATGCTGTATACGTTCAAAAGCAAACCAGCGCCAAACTCACCCTc                                                                          >  2:343989/3‑150 (MQ=255)
                                                                                    ccAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCg                                                                  >  2:413341/1‑150 (MQ=255)
                                                                                     cAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCATCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGt                                                                 <  2:214920/150‑1 (MQ=255)
                                                                                      agtcGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTc                                                                <  2:1804826/147‑1 (MQ=255)
                                                                                      aGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTc                                                                <  2:523326/150‑1 (MQ=255)
                                                                                      aGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTACGCCTGACGCTTCAGCAGGCGGATCTCCCCTCCGAATTGAGGGAAATGCGGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTc                                                                <  1:132985/150‑1 (MQ=255)
                                                                                         cGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAACCTGACCCTCGACTGCCGTCTGc                                                             >  2:422831/1‑150 (MQ=255)
                                                                                         cGTTGAATGATGCGTATCGGCAGTTACGCTAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGc                                                             <  2:314250/150‑1 (MQ=255)
                                                                                          gTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGcc                                                            <  2:1044070/150‑1 (MQ=255)
                                                                                                 gatgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTg                                                     <  1:1431451/150‑1 (MQ=255)
                                                                                                  atgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTgg                                                    <  1:534912/150‑1 (MQ=255)
                                                                                                   tgCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTAACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGc                                                   <  1:99594/150‑1 (MQ=255)
                                                                                                    gCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCa                                                  <  1:1436303/150‑1 (MQ=255)
                                                                                                      ttATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACt                                                <  1:851844/150‑1 (MQ=255)
                                                                                                       tATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTg                                               <  2:1830519/150‑1 (MQ=255)
                                                                                                        aTCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTgg                                              >  1:846539/1‑150 (MQ=255)
                                                                                                          cGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGAt                                            <  1:1164297/150‑1 (MQ=255)
                                                                                                             cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcg                                         <  1:400300/150‑1 (MQ=255)
                                                                                                             cAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATgcg                                         >  1:1874096/1‑150 (MQ=255)
                                                                                                               gTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCa                                       <  1:259841/150‑1 (MQ=255)
                                                                                                               gTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCa                                       <  1:799763/150‑1 (MQ=255)
                                                                                                                   cgcgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATg                                   <  2:912418/150‑1 (MQ=255)
                                                                                                                        gCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGt                              <  1:223386/150‑1 (MQ=255)
                                                                                                                           gTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACTTTACAAAATCAAACCAGCGCCAAACTGACCCTCCACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCa                           >  2:1006087/1‑150 (MQ=255)
                                                                                                                           gTTGACTACTTTCCGCCTGACGCGGCAGCAGGCGGATCTCCCCTCCGCATTGTGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCa                           <  1:779687/150‑1 (MQ=255)
                                                                                                                               actactTTCCGCCTGACGCTTCAGAAGGCGGATCTACACTCCGCATTGAGGTAAATGGTGGATACTTTACAAAATCAAACCAGCGCCAAACTGCCCCTCGACAGCCGTCTGCCACCCCTGGCACTGGATGCGCAATAGCAGGTGACTttg                       >  2:1157906/1‑150 (MQ=255)
                                                                                                                               actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCAtttt                       >  2:1726138/1‑149 (MQ=255)
                                                                                                                               actactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttg                       >  2:1012633/1‑150 (MQ=255)
                                                                                                                                 tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                     <  2:923304/150‑1 (MQ=255)
                                                                                                                                 tactTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCACCTGCATttgcg                     <  2:1383330/150‑3 (MQ=255)
                                                                                                                                 ggagatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                     >  2:1391454/6‑150 (MQ=255)
                                                                                                                                 ggagatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgtt                     >  2:1434806/6‑150 (MQ=255)
                                                                                                                                 aactTTCCTCCTGACGCTTCCGCAGGCGGATCTCCCCTCAGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGTCACTGGATGCGCAAATGCAGGTGCATttgtt                     <  1:1189666/149‑1 (MQ=255)
                                                                                                                                  gagatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATttgttg                    >  2:571813/5‑150 (MQ=255)
                                                                                                                                   agatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGc                   >  2:524563/4‑150 (MQ=255)
                                                                                                                                    tttCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCa                  >  2:958863/1‑150 (MQ=255)
                                                                                                                                    gatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATCCAGGTGCATTTGTTGCa                  >  1:444575/3‑150 (MQ=255)
                                                                                                                                    gatCCGCCTGCTGCAGCGTCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTAACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCa                  >  1:1095532/3‑150 (MQ=255)
                                                                                                                                       ccGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAat               <  1:399037/150‑1 (MQ=255)
                                                                                                                                        cGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAatt              >  2:1033967/1‑150 (MQ=255)
                                                                                                                                           cTGACGCTGCAGCAGGCGCATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAACCCAGCGCCAAACTCACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAattatt           <  2:1748618/150‑1 (MQ=255)
                                                                                                                                             gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcg         <  2:338973/150‑1 (MQ=255)
                                                                                                                                               cGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcgcg       >  2:258693/1‑150 (MQ=255)
                                                                                                                                                 cTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGaa     <  2:134427/150‑1 (MQ=255)
                                                                                                                                                  tGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAg    >  2:1805293/1‑150 (MQ=255)
                                                                                                                                                   gcagcaGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGc   <  2:418940/150‑1 (MQ=255)
                                                                                                                                                    cagcagGCGGATCTCCCCTCCGCATTGAGGGATATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGAGTGGCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCg  <  1:248476/150‑1 (MQ=255)
                                                                                                                                                    cagcagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCg  >  2:1005322/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCG  >  NC_000913/2586864‑2587161

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