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.

TTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  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)
        tagggttgaTAGAAAAATACCACACGATGTATTTTCATAGATTATTCACAACAAAACCGGATGAATTGTGATAAATAAATCAAATAAACAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑gg                                                                                                                                               <  1:1333817/144‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑ggg                                                                                                                                              <  2:1453911/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑ggg                                                                                                                                              <  1:1262369/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑ggg                                                                                                                                              <  2:147465/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGC‑ggg                                                                                                                                              <  1:794912/150‑1 (MQ=255)
         agTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGC‑ggg                                                                                                                                              <  2:1236762/150‑1 (MQ=255)
          gTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGC‑GGGt                                                                                                                                             <  1:1401325/150‑1 (MQ=255)
           tgattaGAGAATATTACCACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTa                                                                                                                                            <  2:825507/150‑1 (MQ=255)
           tgatgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTa                                                                                                                                            >  1:1022223/1‑150 (MQ=255)
           tgatgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAAATCATTCTGC‑GGGTa                                                                                                                                            <  1:986133/150‑1 (MQ=255)
              tgaGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATTATTTCTGC‑GGGTAAGt                                                                                                                                         >  2:937728/1‑150 (MQ=255)
                agagaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GTGTAAGTAt                                                                                                                                       >  2:1424682/1‑150 (MQ=255)
                 gagaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATATCTGC‑GGGTAAGTAtt                                                                                                                                      >  2:1687322/1‑150 (MQ=255)
                  agaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTa                                                                                                                                     <  1:188229/150‑1 (MQ=255)
                   gaATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAg                                                                                                                                    <  1:615594/150‑1 (MQ=255)
                       atTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTa                                                                                                                                >  1:421566/1‑150 (MQ=255)
                            aacaCGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTAttttt                                                                                                                           >  1:1684542/1‑150 (MQ=255)
                             acacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAAACAGGATTAATCTGTATTAATTAATAAGAAAGCAACTTAATACCCGCATAATGATATCAGC‑GGGTAAGTATTATCTTAtatttt                                                                                                                          >  2:1025782/1‑150 (MQ=255)
                              cacGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCTGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTc                                                                                                                         <  1:1449605/150‑1 (MQ=255)
                                cGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGATTATTTTTTcta                                                                                                                       >  1:370597/1‑148 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             >  2:1103612/1‑150 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATcc                                                                                                             >  1:517184/1‑150 (MQ=255)
                                          gCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGATTATTTTTTCGAGCATTAATcc                                                                                                             >  2:86556/1‑150 (MQ=255)
                                           cAGAGATTATGAAGAACTATACCGGATGATTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTATTTTATAAGAAAGCAACTTAAAACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATccc                                                                                                            >  2:110385/1‑150 (MQ=255)
                                           cAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTATGC‑TTGTAAGTATTAGATTATTTTTTCGAGCATTAATccc                                                                                                            >  2:617858/1‑150 (MQ=255)
                                            agagATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcg                                                                                                           <  1:325997/150‑1 (MQ=255)
                                            agagATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCcg                                                                                                           <  2:78224/150‑1 (MQ=255)
                                                 ttATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATTCCGCGCGt                                                                                                      >  1:1658726/1‑150 (MQ=255)
                                                           cTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACg                                                                                            <  1:1350620/150‑1 (MQ=255)
                                                            tataCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGc                                                                                           >  1:837739/1‑150 (MQ=255)
                                                               aCCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAc                                                                                        >  1:1357130/1‑150 (MQ=255)
                                                                ccGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTAcc                                                                                       <  2:19410/150‑1 (MQ=255)
                                                                   gATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAAc                                                                                    <  1:1128188/150‑1 (MQ=255)
                                                                     tGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACgg                                                                                  <  2:865819/150‑1 (MQ=255)
                                                                         ttgTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTa                                                                              <  2:1510084/148‑1 (MQ=255)
                                                                         tGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTa                                                                              >  2:831956/1‑150 (MQ=255)
                                                                          ggTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACACTACCAACGGATTAg                                                                             >  2:694437/1‑150 (MQ=255)
                                                                            tGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAggg                                                                           <  1:287851/150‑1 (MQ=255)
                                                                               taaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCTCGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                        >  2:344554/1‑150 (MQ=255)
                                                                               taaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAAt                                                                        >  1:1219448/1‑150 (MQ=255)
                                                                                aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCTAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  2:1657129/1‑150 (MQ=255)
                                                                                aaataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATa                                                                       >  2:980580/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGATCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      <  1:706938/150‑1 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      <  2:46209/150‑1 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  1:1710100/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  2:1249216/1‑150 (MQ=255)
                                                                                 aataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATaa                                                                      >  1:722715/1‑150 (MQ=255)
                                                                                  ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAATCAACTTAATACCCGCAGAATGATTTCTGC‑TGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAAAGGAGTAGGGAATAac                                                                     >  2:152343/1‑150 (MQ=255)
                                                                                  ataaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAac                                                                     <  2:722715/150‑1 (MQ=255)
                                                                                    aaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAacac                                                                   >  1:1711778/1‑150 (MQ=255)
                                                                                    aaaGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAacac                                                                   <  2:761991/150‑1 (MQ=255)
                                                                                      aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCACCGCTACCAACGGATTAGGGAATAACACAt                                                                 <  2:699727/150‑1 (MQ=255)
                                                                                      aGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACAt                                                                 <  1:971494/150‑1 (MQ=255)
                                                                                       tcAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGCGGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACAt                                                                 <  1:466377/149‑1 (MQ=255)
                                                                                       gCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCAATAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATa                                                                >  2:1342459/1‑150 (MQ=255)
                                                                                        cAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTACTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATAc                                                               >  2:435428/1‑150 (MQ=255)
                                                                                          aaTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActc                                                             >  2:790804/1‑150 (MQ=255)
                                                                                           aTAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATActct                                                            <  2:196425/150‑1 (MQ=255)
                                                                                            tAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTa                                                           >  1:879170/1‑150 (MQ=255)
                                                                                             aaCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAAACCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAc                                                          <  2:621685/150‑1 (MQ=255)
                                                                                              aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         >  1:592035/1‑150 (MQ=255)
                                                                                              aCCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTAca                                                         <  1:1616001/150‑1 (MQ=255)
                                                                                                      tAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGt                                                 <  1:1424682/150‑1 (MQ=255)
                                                                                                       taTCTGTATTAATTTTATCCAAATCAACATAATACCCGCCGAATGCATTATGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTCATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTTCATCATGGGtt                                                <  2:329715/149‑1 (MQ=255)
                                                                                                       aaTTTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGtt                                                <  1:594122/150‑1 (MQ=255)
                                                                                                         tCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑TGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGTGAATAACACATACTCTACATCATGGGTTAc                                              >  1:1631589/1‑150 (MQ=255)
                                                                                                           tGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACgg                                            >  2:1139348/1‑150 (MQ=255)
                                                                                                             tATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTa                                          >  1:875047/1‑150 (MQ=255)
                                                                                                               ttAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTaaa                                        <  2:1711778/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      <  1:1049107/150‑1 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCAAAACGCTACCAATGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      >  1:278686/1‑150 (MQ=255)
                                                                                                                 aaTTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGAATTAATCCCGCGCGTAATCCCAACGCTACAAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAAc                                      >  1:1025734/1‑150 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:462137/150‑1 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 <  2:837739/150‑1 (MQ=255)
                                                                                                                      ataAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTcc                                 >  2:778113/1‑150 (MQ=255)
                                                                                                                        aaGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctc                               <  1:831956/150‑1 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGTGTTACGGTAAAACGTTCctct                              >  2:673553/1‑150 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCctct                              <  1:538845/150‑1 (MQ=255)
                                                                                                                         aGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GAGTAAGTATAAGCTAATTTTTTCGAGCATTAATACCTCGCGTAATACAAACGCTAACAACGGATTAGGGAATAACACAAACTCTACATCATGGATTACGATAAAACGTTCatct                              >  2:927647/1‑150 (MQ=255)
                                                                                                                           aaaGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTcc                            >  2:174293/1‑150 (MQ=255)
                                                                                                                            aatcAACTTAATTCCCTCATAATGATTTCTGC‑GTGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATTGGTTACGGTAAAACGTTCCTCTCCg                           <  2:370597/147‑1 (MQ=255)
                                                                                                                            aaGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACAAATACTCTACATAATGGGTTACGGTAAAACGTTCCTCTCCg                           >  2:1143321/1‑150 (MQ=255)
                                                                                                                                 aCTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTg                      <  1:527646/150‑1 (MQ=255)
                                                                                                                                  cTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTACTCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgc                     >  1:117696/1‑150 (MQ=255)
                                                                                                                                    acttACCCGCAGAAATCATTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTgcgc                   >  2:1444967/4‑150 (MQ=255)
                                                                                                                                        aCCCGCAGAAATCATTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACTGATTAGGGAATAACACATACTCTACATCATCGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGc               >  1:790528/1‑150 (MQ=255)
                                                                                                                                         cccGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              >  1:511627/1‑150 (MQ=255)
                                                                                                                                         cccGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCa              <  1:937728/150‑1 (MQ=255)
                                                                                                                                          ccGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCaa             >  2:382695/1‑150 (MQ=255)
                                                                                                                                           cGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAAt            >  2:314850/1‑150 (MQ=255)
                                                                                                                                            gCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg           <  1:1575980/150‑1 (MQ=255)
                                                                                                                                            gCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATg           <  2:643201/150‑1 (MQ=255)
                                                                                                                                              aGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGtt         >  1:283655/1‑150 (MQ=255)
                                                                                                                                                 aTGATTTCTGC‑GGTTAAGTATTAGCTTATTTTTTCGAGCATTAAGCCAGCGAGTAATCCCAACTCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCt      <  1:110385/150‑1 (MQ=255)
                                                                                                                                                  tGATTTCTGC‑GTGTAAGTATTAGCTTATTTTATCGAGCATTAATCCCGCGCGTAATCCCAACGCGACCAACGGATTAGGGAATAACACATACTCTACATCATGGGGTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCtt     <  1:86426/150‑1 (MQ=255)
                                                                                                                                                     tttCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCACCGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  2:852498/150‑1 (MQ=255)
                                                                                                                                                     tttCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTtct  <  2:596792/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTGATGGAGAGTGATGAGAGAATATTACAACACGATGATTTTGCAGAGATTATGAAGAACTATACCGGATGACTGGTGATAAATAAAGCAAATAACCAGGATTAATCTGTATTAATTTATAAGAAAGCAACTTAATACCCGCAGAATGATTTCTGC‑GGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCT  >  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: ‑