Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,689,897 A→T 6.1% V401D (GTC→GAC)  glrK ← sensory histidine kinase GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,689,8970AT 6.1% 260.5 / 5.0 82V401D (GTC→GAC) glrKsensory histidine kinase GlrK
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base T (3/2);  total (40/42)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

CCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGA  >  NC_000913/2689749‑2690041
                                                                                                                                                    |                                                                                                                                                
cccTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGAc                                                                                                                                                 <  2:1580883/150‑1 (MQ=255)
   ttGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATc                                                                                                                                              <  1:1789637/150‑1 (MQ=255)
          cGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGTCATCAatacag                                                                                                                                       <  1:172877/150‑4 (MQ=255)
          cGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGTCATCAatacagg                                                                                                                                      <  1:1600690/150‑5 (MQ=255)
             cGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAcc                                                                                                                                    <  1:434439/150‑1 (MQ=255)
              gAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAAccc                                                                                                                                   >  1:1048595/1‑150 (MQ=255)
                 tACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCcgcg                                                                                                                                >  1:1245418/1‑150 (MQ=255)
                   cAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATTTAAACCCGCGCa                                                                                                                              >  1:1579021/1‑150 (MQ=255)
                   cAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCa                                                                                                                              >  2:129529/1‑150 (MQ=255)
                    aaTCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCAc                                                                                                                             >  2:397481/1‑150 (MQ=255)
                      tCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCa                                                                                                                           <  1:212164/150‑1 (MQ=255)
                        ccTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATg                                                                                                                         <  2:1048595/150‑1 (MQ=255)
                        ccTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGGGCCTGTATTGATGACATCAATATAAACCCGCGCACCATg                                                                                                                         >  2:1649032/1‑150 (MQ=255)
                        ccTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCACCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATg                                                                                                                         <  1:320338/150‑1 (MQ=255)
                         cTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCACGt                                                                                                                        <  2:324421/150‑1 (MQ=255)
                          tGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTa                                                                                                                       >  1:300985/1‑150 (MQ=255)
                                 gCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAActgct                                                                                                                >  1:1253646/1‑150 (MQ=255)
                                   ttAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACtgctgc                                                                                                              <  1:131375/150‑1 (MQ=255)
                                    tAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCg                                                                                                             >  1:276947/1‑150 (MQ=255)
                                      aTCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGat                                                                                                           <  1:67726/150‑2 (MQ=255)
                                      aTCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGaa                                                                                                           <  1:67723/150‑1 (MQ=255)
                                       tCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAg                                                                                                          >  1:281283/1‑150 (MQ=255)
                                        cccAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAgg                                                                                                         <  2:1224139/150‑1 (MQ=255)
                                         ccAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGc                                                                                                        >  2:1503580/1‑150 (MQ=255)
                                            gACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCaaa                                                                                                     >  2:321726/1‑150 (MQ=255)
                                              ccGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATg                                                                                                   <  2:1479282/150‑1 (MQ=255)
                                                 cTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTa                                                                                                >  1:1722246/1‑150 (MQ=255)
                                                   gCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTAcc                                                                                              <  1:1503580/150‑1 (MQ=255)
                                                    cccTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCg                                                                                             <  1:174583/150‑1 (MQ=255)
                                                         cACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTc                                                                                        <  2:1316828/150‑1 (MQ=255)
                                                         cACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTc                                                                                        >  1:1778021/1‑150 (MQ=255)
                                                           ccGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCac                                                                                      >  2:393311/1‑149 (MQ=255)
                                                           ccGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAg                                                                                      >  2:807867/1‑150 (MQ=255)
                                                            cGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGc                                                                                     <  1:76721/150‑1 (MQ=255)
                                                             gCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGcc                                                                                    <  2:32865/150‑1 (MQ=255)
                                                             gCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCACATGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGcc                                                                                    >  2:1679463/1‑150 (MQ=255)
                                                                ccTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCg                                                                                 <  1:1647415/150‑1 (MQ=255)
                                                                ccTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCg                                                                                 >  1:207088/1‑150 (MQ=255)
                                                                     tCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTg                                                                            <  1:397481/150‑1 (MQ=255)
                                                                          ggtggCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCg                                                                       <  2:276947/150‑1 (MQ=255)
                                                                           gtggCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGc                                                                      >  2:1443577/1‑150 (MQ=255)
                                                                            tggCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCa                                                                     >  2:1878528/1‑150 (MQ=255)
                                                                             ggCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCg                                                                                             <  1:436222/125‑1 (MQ=255)
                                                                             ggCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCg                                                                                             >  2:436222/1‑125 (MQ=255)
                                                                              gCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCAtt                                                                   <  1:1388198/150‑1 (MQ=255)
                                                                               cTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTg                                                                  >  2:968728/1‑150 (MQ=255)
                                                                               cTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTg                                                                  <  2:974138/150‑1 (MQ=255)
                                                                                 tCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCa                                                                          >  2:675562/1‑140 (MQ=255)
                                                                                 tCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCa                                                                          <  1:675562/140‑1 (MQ=255)
                                                                                    cTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTa                                                             >  2:1128934/1‑150 (MQ=255)
                                                                                    cTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTa                                                             >  1:1845629/1‑150 (MQ=255)
                                                                                    cTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTa                                                             >  1:56605/1‑150 (MQ=255)
                                                                                    cTGAAAAAAGGGTTCGAAGATCATGGCCCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTa                                                             >  2:1402572/1‑150 (MQ=255)
                                                                                      gAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTg                                                                  >  1:1830410/1‑143 (MQ=255)
                                                                                      gAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTg                                                                  <  2:1830410/143‑1 (MQ=255)
                                                                                        aaaaaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGa                                                         >  1:660521/1‑150 (MQ=255)
                                                                                         aaacGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCGGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGAt                                                        <  2:1139970/150‑1 (MQ=255)
                                                                                         aaaaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATTTAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGAt                                                        <  2:1579021/150‑1 (MQ=255)
                                                                                         aaaaGGGATCCAAGATCATGCCCCGATACTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGAt                                                        >  2:1911977/1‑150 (MQ=255)
                                                                                          aaaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGAtt                                                       >  2:1452445/1‑150 (MQ=255)
                                                                                           aaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTa                                                      >  1:1035272/1‑150 (MQ=255)
                                                                                           aaGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTa                                                      <  1:916012/150‑1 (MQ=255)
                                                                                              ggTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCGGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATcc                                                   >  1:1990870/1‑150 (MQ=255)
                                                                                              ggTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATcc                                                   <  1:532671/150‑1 (MQ=255)
                                                                                                   gAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTAc                                              >  2:453554/1‑150 (MQ=255)
                                                                                                          aTGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCtcatc                                       <  2:1845629/150‑1 (MQ=255)
                                                                                                          aTGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCtcatc                                       <  1:129529/150‑1 (MQ=255)
                                                                                                             gcgcGTTCCTCTTGCGGAATGGGCGTGCCTGTAATGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATcagc                                    <  2:37325/150‑1 (MQ=255)
                                                                                                                      tctTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGc                           <  1:453554/150‑1 (MQ=255)
                                                                                                                         tGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTg                                                                            <  2:795043/98‑1 (MQ=255)
                                                                                                                         tGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTg                                                                            >  1:795043/1‑98 (MQ=255)
                                                                                                                          gCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCg                       <  2:660521/150‑1 (MQ=255)
                                                                                                                               aTGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCGCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAgg                  <  2:300985/150‑1 (MQ=255)
                                                                                                                                      ggtttatATTGATGTCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTg           >  1:423564/7‑150 (MQ=255)
                                                                                                                                       gCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGt          <  2:836422/150‑1 (MQ=255)
                                                                                                                                        tttatATTGATGTCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGtt         >  2:1613022/5‑150 (MQ=255)
                                                                                                                                        tttatATTGATGTCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGtt         >  2:380677/5‑150 (MQ=255)
                                                                                                                                          tGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGc       <  1:807867/150‑1 (MQ=255)
                                                                                                                                           gTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCt      >  2:15917/1‑150 (MQ=255)
                                                                                                                                             aTTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCttt    >  1:1189208/1‑150 (MQ=255)
                                                                                                                                              ttGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTg   <  2:1778021/150‑1 (MQ=255)
                                                                                                                                              ttGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTg   <  1:1128934/150‑1 (MQ=255)
                                                                                                                                               tgatgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTga  >  2:413893/1‑150 (MQ=255)
                                                                                                                                               tgatgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTga  <  1:15917/150‑1 (MQ=255)
                                                                                                                                               tgatgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTga  >  1:30943/1‑150 (MQ=255)
                                                                                                                                               tgatgaCATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGCGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTgc  <  1:1501071/150‑2 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCCTTGCATACGGCGAATACAATCCCTGGCAATGCTTAATCCCAGACCGCTGCCCTTCACCGCCCCTTTTCGCTGGTGGCTTCCCTGAAAAAAGGGTTCGAAGATCATGGCGCGTTCCTCTTGCGGAATGGGCGTGCCTGTATTGATGACATCAATATAAACCCGCGCACCATGTAAACTGCTGCGAAGGCAAATGTTACCGGATTCAGCCCCGTAGTGCACCGCATTGGAGTAAAGATTATCCAGTACGCTCATCAGCAGCATTGGCTCCGCCAGGCAAGCTGTTGCTTTGA  >  NC_000913/2689749‑2690041

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