Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,825,939 T→A 7.1% intergenic (‑314/+16) spy ← / ← astE periplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,825,9390TA 7.1% 312.0 / 9.7 110intergenic (‑314/+16)spy/astEperiplasmic ATP‑independent protein refolding chaperone, stress‑induced/succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base T (53/49);  new base A (3/5);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 114 total aligned reads displayed.

TTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  NC_000913/1825790‑1826088
                                                                                                                                                     |                                                                                                                                                      
ttGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAt                                                                                                                                                        <  1:1362280/150‑1 (MQ=255)
 tGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGAtt                                                                                                                                                       >  2:1350620/1‑150 (MQ=255)
     ggAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTg                                                                                                                                                   <  2:1451706/150‑1 (MQ=255)
        tagggttgaTAGAAAAATACCACACGATGTATTTTCATAGATTATTCACAACAAAACCGGATGAATTGTGATAAATAAATCAAATAAACAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCgg                                                                                                                                               <  1:1333817/144‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                              <  2:1453911/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                              <  1:1262369/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCggg                                                                                                                                              <  2:147465/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                              <  1:794912/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCggg                                                                                                                                              <  2:1236762/150‑1 (MQ=255)
          gTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCGGGt                                                                                                                                             <  1:1401325/150‑1 (MQ=255)
           tgattaGAGAATATTACCACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTa                                                                                                                                            <  2:825507/150‑1 (MQ=255)
           tgatgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTa                                                                                                                                            >  1:1022223/1‑150 (MQ=255)
           tgatgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGCGGGTa                                                                                                                                            <  1:986133/150‑1 (MQ=255)
              tgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATTATTTCTGCGGGTAAGt                                                                                                                                         >  2:937728/1‑150 (MQ=255)
                agagaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGTGTAAGTAt                                                                                                                                       >  2:1424682/1‑150 (MQ=255)
                 gagaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATATCTGCGGGTAAGTAtt                                                                                                                                      >  2:1687322/1‑150 (MQ=255)
                  agaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTa                                                                                                                                     <  1:188229/150‑1 (MQ=255)
                   gaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAg                                                                                                                                    <  1:615594/150‑1 (MQ=255)
                       atTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTa                                                                                                                                >  1:421566/1‑150 (MQ=255)
                            aacaCGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTAttttt                                                                                                                           >  1:1684542/1‑150 (MQ=255)
                             acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAAACAGGATTAATCTGTATTAATTAATAAGAAAGCAACTTAATACCCGCATAATGATATCAGCGGGTAAGTATTATCTTAtatttt                                                                                                                          >  2:1025782/1‑150 (MQ=255)
                              cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCTGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTc                                                                                                                         <  1:1449605/150‑1 (MQ=255)
                                cGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGATTATTTTTTcta                                                                                                                       >  1:370597/1‑148 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             >  2:1103612/1‑150 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             >  1:517184/1‑150 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGATTATTTTTTCGAGCATTAATcc                                                                                                             >  2:86556/1‑150 (MQ=255)
                                           cAGAGATTATGAAGAACTATACCGGATGATTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTATTTTATAAGAAAGCAACTTAAAACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATccc                                                                                                            >  2:110385/1‑150 (MQ=255)
                                           cAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTATGCTTGTAAGTATTAGATTATTTTTTCGAGCATTAATccc                                                                                                            >  2:617858/1‑150 (MQ=255)
                                            agagATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcg                                                                                                           <  1:325997/150‑1 (MQ=255)
                                            agagATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcg                                                                                                           <  2:78224/150‑1 (MQ=255)
                                                 ttATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATTCCGCGCGt                                                                                                      >  1:1658726/1‑150 (MQ=255)
                                                           cTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACg                                                                                            <  1:1350620/150‑1 (MQ=255)
                                                            tataCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGc                                                                                           >  1:837739/1‑150 (MQ=255)
                                                               aCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAc                                                                                        >  1:1357130/1‑150 (MQ=255)
                                                                ccGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  2:19410/150‑1 (MQ=255)
                                                                   gATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAAc                                                                                    <  1:1128188/150‑1 (MQ=255)
                                                                     tGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACgg                                                                                  <  2:865819/150‑1 (MQ=255)
                                                                         ttgTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTa                                                                              <  2:1510084/148‑1 (MQ=255)
                                                                         tGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTa                                                                              >  2:831956/1‑150 (MQ=255)
                                                                          ggTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACACTACCAACGGATTAg                                                                             >  2:694437/1‑150 (MQ=255)
                                                                            tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                           <  1:287851/150‑1 (MQ=255)
                                                                               taaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                        >  2:344554/1‑150 (MQ=255)
                                                                               taaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                        >  1:1219448/1‑150 (MQ=255)
                                                                                aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCTAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  2:1657129/1‑150 (MQ=255)
                                                                                aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  2:980580/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGATCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      <  1:706938/150‑1 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      <  2:46209/150‑1 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  1:1710100/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  2:1249216/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  1:722715/1‑150 (MQ=255)
                                                                                  ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAATCAACTTAATACCCGCAGAATGATTTCTGCTGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAAAGGAGTAGGGAATAac                                                                     >  2:152343/1‑150 (MQ=255)
                                                                                  ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                     <  2:722715/150‑1 (MQ=255)
                                                                                    aaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAacac                                                                   >  1:1711778/1‑150 (MQ=255)
                                                                                    aaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAacac                                                                   <  2:761991/150‑1 (MQ=255)
                                                                                      aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCACCGCTACCAACGGATTAGGGAATAACACAt                                                                 <  2:699727/150‑1 (MQ=255)
                                                                                      aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACAt                                                                 <  1:971494/150‑1 (MQ=255)
                                                                                       tcAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACAt                                                                 <  1:466377/149‑1 (MQ=255)
                                                                                       gCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCAATAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                >  2:1342459/1‑150 (MQ=255)
                                                                                        cAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTACTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATAc                                                               >  2:435428/1‑150 (MQ=255)
                                                                                          aaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActc                                                             >  2:790804/1‑150 (MQ=255)
                                                                                           aTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActct                                                            <  2:196425/150‑1 (MQ=255)
                                                                                            tAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTa                                                           >  1:879170/1‑150 (MQ=255)
                                                                                             aaCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAc                                                          <  2:621685/150‑1 (MQ=255)
                                                                                              aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         >  1:592035/1‑150 (MQ=255)
                                                                                              aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         <  1:1616001/150‑1 (MQ=255)
                                                                                                      tAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGt                                                 <  1:1424682/150‑1 (MQ=255)
                                                                                                       taTCTGTATTAATTTTATCCAAATCAACATAATACCCGCCGAATGCATTATGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTCATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTTCATCATGGGtt                                                <  2:329715/149‑1 (MQ=255)
                                                                                                       aaTTTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGtt                                                <  1:594122/150‑1 (MQ=255)
                                                                                                         tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCTGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGTGAATAACACATACTCTACATCATGGGTTAc                                              >  1:1631589/1‑150 (MQ=255)
                                                                                                           tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                            >  2:1139348/1‑150 (MQ=255)
                                                                                                             tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTa                                          >  1:875047/1‑150 (MQ=255)
                                                                                                               ttAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaa                                        <  2:1711778/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      <  1:1049107/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCAAAACGCTACCAATGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      >  1:278686/1‑150 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGAATTAATCCCGCGCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      >  1:1025734/1‑150 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:462137/150‑1 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:837739/150‑1 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 >  2:778113/1‑150 (MQ=255)
                                                                                                                        aaGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctc                               <  1:831956/150‑1 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGTGTTACGGTAAAACGTTCctct                              >  2:673553/1‑150 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctct                              <  1:538845/150‑1 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGAGTAAGTATAAGCTAATTTTTTCGAGCATTAATACCTCGCGTAATACAAACGCTAACAACGGATTAGGGAATAACACAAACTCTACATCATGGATTACGATAAAACGTTCatct                              >  2:927647/1‑150 (MQ=255)
                                                                                                                           aaaGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTcc                            >  2:174293/1‑150 (MQ=255)
                                                                                                                            aatcAACTTAATTCCCTCATAATGATTTCTGCGTGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATTGGTTACGGTAAAACGTTCCTCTCCg                           <  2:370597/147‑1 (MQ=255)
                                                                                                                            aaGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACAAATACTCTACATAATGGGTTACGGTAAAACGTTCCTCTCCg                           >  2:1143321/1‑150 (MQ=255)
                                                                                                                                 aCTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTg                      <  1:527646/150‑1 (MQ=255)
                                                                                                                                  cTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTACTCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                     >  1:117696/1‑150 (MQ=255)
                                                                                                                                    acttACCCGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcgc                   >  2:1444967/4‑150 (MQ=255)
                                                                                                                                        aCCCGCAGAAATCATTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACTGATTAGGGAATAACACATACTCTACATCATCGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGc               >  1:790528/1‑150 (MQ=255)
                                                                                                                                         cccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              >  1:511627/1‑150 (MQ=255)
                                                                                                                                         cccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              <  1:937728/150‑1 (MQ=255)
                                                                                                                                          ccGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCaa             >  2:382695/1‑150 (MQ=255)
                                                                                                                                           cGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAAt            >  2:314850/1‑150 (MQ=255)
                                                                                                                                            gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg           <  1:1575980/150‑1 (MQ=255)
                                                                                                                                            gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg           <  2:643201/150‑1 (MQ=255)
                                                                                                                                              aGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGtt         >  1:283655/1‑150 (MQ=255)
                                                                                                                                                 aTGATTTCTGCGGTTAAGTATTAGCTTATTTTTTCGAGCATTAAGCCAGCGAGTAATCCCAACTCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      <  1:110385/150‑1 (MQ=255)
                                                                                                                                                  tGATTTCTGCGTGTAAGTATTAGCTTATTTTATCGAGCATTAATCCCGCGCGTAATCCCAACGCGACCAACGGATTAGGGAATAACACATACTCTACATCATGGGGTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCtt     <  1:86426/150‑1 (MQ=255)
                                                                                                                                                     tttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCACCGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  2:852498/150‑1 (MQ=255)
                                                                                                                                                     tttCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  2:596792/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  NC_000913/1825790‑1826088

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