Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,587,014 G→T 10.0% Q428H (CAG→CAT) ‡ 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,0140GT10.0% 275.8 / 17.9 105Q428H (CAG→CAT) ‡narQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base G (37/57);  new base T (7/4);  total (44/61)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.09e-01

Only 100 of 114 total aligned reads displayed.

TACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGT  >  NC_000913/2586865‑2587163
                                                                                                                                                     |                                                                                                                                                     
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)
              gTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCAGATCTccc                                                                                                                                         <  2:1737859/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                                                                                                           >  2:392198/1‑150 (MQ=255)
                                            aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                           <  1:890805/150‑1 (MQ=255)
                                            aGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGtt                                                                                                           <  2:1003009/150‑1 (MQ=255)
                                                  cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaaa                                                                                                     <  2:1573269/150‑1 (MQ=255)
                                                  cGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACaaaa                                                                                                     >  1:1498081/1‑150 (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)
                                                                            aTTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGATGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCg                                                                           >  2:234367/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)
                                                                                              atgatgCTTATCGGCAGTTACGCGAGCTGGTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAAccc                                                         <  2:405678/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)
                                                                                                     ttATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACt                                                  <  2:1495353/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)
                                                                                                                    cgAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCa                                   <  2:138618/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: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)
                                                                                                                                            gACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTcg           <  2:338973/150‑1 (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)
                                                                                                                                                     gcagGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGt  >  2:455655/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACTTACGTATCCAGTTGACGTTACTGAAGCGTTCGATACCGGAAGATAACGCCACCGCACAAAGTATCATGGCCGATTTTTCCCAGGCGTTGAATGATGCTTATCGGCAGTTACGCGAGCTGTTGACTACTTTCCGCCTGACGCTGCAGCAGGCGGATCTCCCCTCCGCATTGAGGGAAATGCTGGATACGTTACAAAATCAAACCAGCGCCAAACTGACCCTCGACTGCCGTCTGCCAACCCTGGCACTGGATGCGCAAATGCAGGTGCATTTGTTGCAAATTATTCGCGAAGCGGT  >  NC_000913/2586865‑2587163

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: