Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,585,998 T→C Y90H (TAT→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,585,9980TC99.1% 432.6 / ‑4.3 111Y90H (TAT→CAT) narQsensory histidine kinase in two‑component regulatory system with NarP
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (53/57);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

Only 100 of 111 total aligned reads displayed.

TGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTT  >  NC_000913/2585850‑2586145
                                                                                                                                                    |                                                                                                                                                   
tGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCAt                                                                                                                                                   >  1:62463/1‑151 (MQ=255)
tGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCAt                                                                                                                                                   >  2:374090/1‑151 (MQ=255)
 gAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCAt                                                                                                                                                   <  2:582108/150‑1 (MQ=255)
   ggCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCAtgtg                                                                                                                                                <  2:1292599/151‑1 (MQ=255)
   ggCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCAtgt                                                                                                                                                 <  2:489467/150‑1 (MQ=255)
         cAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGaa                                                                                                                                          <  1:1340201/151‑1 (MQ=255)
            tatTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGc                                                                                                                                        >  2:1020237/1‑150 (MQ=255)
            tatTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGc                                                                                                                                        >  1:1305381/1‑150 (MQ=255)
              ttGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGc                                                                                                                                        >  1:1060405/1‑148 (MQ=255)
              ttGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGc                                                                                                                                        <  2:1060487/148‑1 (MQ=255)
              ttGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGt                                                                                                                                     >  1:936399/1‑151 (MQ=255)
                gCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTa                                                                                                                                    >  2:853422/1‑150 (MQ=255)
                    gATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAAc                                                                                                                               <  1:482970/151‑1 (MQ=255)
                    gATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAAc                                                                                                                               <  2:809448/151‑1 (MQ=255)
                     aTCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCTCATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACt                                                                                                                              >  1:1372901/1‑151 (MQ=255)
                     aTCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACt                                                                                                                              >  1:1001539/1‑151 (MQ=255)
                      acGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTc                                                                                                                             >  1:1101220/2‑151 (MQ=255)
                       cGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCg                                                                                                                            >  1:1331508/1‑151 (MQ=255)
                       cGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCg                                                                                                                            <  2:1670439/151‑1 (MQ=255)
                       cGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCg                                                                                                                            <  1:1216113/151‑1 (MQ=255)
                         cTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCt                                                                                                                          >  2:155376/1‑151 (MQ=255)
                            cGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATg                                                                                                                       >  1:1017528/1‑151 (MQ=255)
                             gTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATgc                                                                                                                      <  2:1447459/151‑1 (MQ=255)
                             gTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATg                                                                                                                       >  1:185649/1‑150 (MQ=255)
                                 gCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCa                                                                                                                   <  2:300402/150‑1 (MQ=255)
                                 gCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCAt                                                                                                                  >  2:195721/1‑151 (MQ=255)
                                   agagTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCt                                                                                                                <  2:1419285/151‑1 (MQ=255)
                                       ttACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAAt                                                                                                            >  1:179416/1‑151 (MQ=255)
                                        tACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATg                                                                                                           >  2:1161998/1‑151 (MQ=255)
                                            gCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCaa                                                                                                       <  1:424565/151‑1 (MQ=255)
                                              cTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAAct                                                                                                     >  1:281760/1‑151 (MQ=255)
                                              cTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAAct                                                                                                     >  2:1311618/1‑151 (MQ=255)
                                                gggCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAActg                                                                                                    <  2:221338/150‑1 (MQ=255)
                                                  gCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAActggct                                                                                                 <  1:1436731/151‑1 (MQ=255)
                                                           tGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGa                                                                                         >  1:1059149/1‑150 (MQ=255)
                                                           tGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTAACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGa                                                                                         >  1:886772/1‑150 (MQ=255)
                                                            gCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGaat                                                                                       >  2:117922/1‑151 (MQ=255)
                                                            gCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGaat                                                                                       >  2:1225925/1‑151 (MQ=255)
                                                              aaaGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGaata                                                                                      >  2:1663564/1‑150 (MQ=255)
                                                              aaaGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGaata                                                                                      >  1:965847/1‑150 (MQ=255)
                                                               aaGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGaataa                                                                                     >  2:1127453/1‑150 (MQ=255)
                                                                 gTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCg                                                                                  <  1:1640275/151‑1 (MQ=255)
                                                                  tGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCgg                                                                                 >  1:1543726/1‑151 (MQ=255)
                                                                   ggCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGc                                                                                <  1:1308146/151‑1 (MQ=255)
                                                                      aGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTc                                                                              >  2:1592766/1‑150 (MQ=255)
                                                                       gTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAg                                                                            <  2:1469498/151‑1 (MQ=255)
                                                                        tCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGc                                                                           >  2:1614455/1‑151 (MQ=255)
                                                                         ccACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCa                                                                          <  1:1164429/151‑1 (MQ=255)
                                                                                tCAATGCACATAGCCCGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCg                                                                    >  1:1333252/1‑150 (MQ=255)
                                                                                   aTGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGAttt                                                                <  1:1305271/151‑1 (MQ=255)
                                                                                    tGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTg                                                               >  1:1089687/1‑151 (MQ=255)
                                                                                     gCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTg                                                               <  1:1154184/150‑1 (MQ=255)
                                                                                     gCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTg                                                               <  1:39196/150‑1 (MQ=255)
                                                                                        caTCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGt                                                           <  2:1338160/151‑1 (MQ=255)
                                                                                          tCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTgg                                                         >  2:1575835/1‑151 (MQ=255)
                                                                                           cGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGt                                                        >  2:1340057/1‑151 (MQ=255)
                                                                                            gCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGt                                                        >  2:377505/1‑150 (MQ=255)
                                                                                               aGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCa                                                    <  2:1113436/151‑1 (MQ=255)
                                                                                                gTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAg                                                   <  1:482808/151‑1 (MQ=255)
                                                                                                   aTTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGcc                                                >  1:234753/1‑151 (MQ=255)
                                                                                                   aTTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGc                                                 <  2:59568/150‑1 (MQ=255)
                                                                                                     ttCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCa                                               <  1:1050708/150‑1 (MQ=255)
                                                                                                       cagcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAata                                            <  2:1173912/151‑1 (MQ=255)
                                                                                                       cagcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAata                                            <  2:562587/151‑1 (MQ=255)
                                                                                                       cagcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAat                                             >  2:1306142/1‑150 (MQ=255)
                                                                                                        agcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAatat                                           <  2:1697332/151‑1 (MQ=255)
                                                                                                         gcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATAtt                                          <  1:553167/151‑1 (MQ=255)
                                                                                                         gcagGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATAtt                                          >  1:568548/1‑151 (MQ=255)
                                                                                                           agGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaa                                        <  2:1577138/151‑1 (MQ=255)
                                                                                                           agGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaa                                        <  2:1059231/151‑1 (MQ=255)
                                                                                                               aCTGCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATtaataa                                     >  1:772628/1‑150 (MQ=255)
                                                                                                                  gCATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTAt                                 <  1:168631/151‑1 (MQ=255)
                                                                                                                   cATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATg                                <  2:62467/151‑1 (MQ=255)
                                                                                                                   cATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATg                                <  2:1553050/151‑1 (MQ=255)
                                                                                                                   cATTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATg                                <  2:694179/151‑1 (MQ=255)
                                                                                                                    aTTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGt                               <  2:39989/151‑1 (MQ=255)
                                                                                                                    aTTCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGt                               <  2:1305475/151‑1 (MQ=255)
                                                                                                                     ttCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGtt                              >  1:848805/1‑151 (MQ=255)
                                                                                                                      tCACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGt                               <  1:374062/149‑1 (MQ=255)
                                                                                                                       cACCGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTa                             <  1:680529/150‑1 (MQ=255)
                                                                                                                         ccGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATc                          <  2:1101306/151‑1 (MQ=255)
                                                                                                                          cGGTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATc                          >  2:820839/1‑150 (MQ=255)
                                                                                                                           ggTATTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATc                          >  1:654072/1‑149 (MQ=255)
                                                                                                                              aTTAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATa                     <  2:1017605/151‑1 (MQ=255)
                                                                                                                               ttAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGAt                      >  1:240778/1‑149 (MQ=255)
                                                                                                                                tAACCAACCTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAg                    <  2:1449273/150‑1 (MQ=255)
                                                                                                                                 aaccaaccTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTCCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGAc                  >  2:1171232/1‑151 (MQ=255)
                                                                                                                                    caaccTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTg               <  1:1668239/151‑1 (MQ=255)
                                                                                                                                    caaccTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTg               >  1:1024628/1‑151 (MQ=255)
                                                                                                                                    caaccTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTg               <  1:6653/151‑1 (MQ=255)
                                                                                                                                      accTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGtt             <  1:1663456/151‑1 (MQ=255)
                                                                                                                                      accTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGtt             <  2:415069/151‑1 (MQ=255)
                                                                                                                                      accTCAACGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGtt             >  2:1125785/1‑151 (MQ=255)
                                                                                                                                           aaCGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTAc        >  1:809564/1‑151 (MQ=255)
                                                                                                                                            aCGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACt       <  1:1311524/151‑1 (MQ=255)
                                                                                                                                            aCGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACt       <  1:430240/151‑1 (MQ=255)
                                                                                                                                             cGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTg      <  1:520623/151‑1 (MQ=255)
                                                                                                                                             cGTCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAAACGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTg      >  1:736822/1‑151 (MQ=255)
                                                                                                                                               tCTGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGc    >  1:1543589/1‑151 (MQ=255)
                                                                                                                                                 tGGCATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCtt  >  2:873411/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGAGGCTATCAATATTGCCGGATCGCTGCGTATGCAGAGTTACCGCCTGGGCTACGACTTGCAAAGTGGCAGTCCACAACTCAATGCACATCGCCAGTTATTTCAGCAGGCACTGCATTCACCGGTATTAACCAACCTCAACGTCTGGTATGTGCCAGAAGCAGTAAAAACTCGCTATGCGCATCTGAATGCCAACTGGCTGGAGATGAATAATCGGCTCAGCAAGGGCGATTTGCCGTGGTATCAGGCCAATATTAATAATTATGTTAATCAGATAGACCTGTTCGTACTGGCTT  >  NC_000913/2585850‑2586145

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