Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,415,726 T→G 7.0% L327R (CTG→CGG) ‡ pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,415,7260TG 7.0% 461.1 / 20.7 129L327R (CTG→CGG) ‡ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (58/62);  new base G (4/5);  total (62/67)
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.83e-01

Only 100 of 131 total aligned reads displayed.

TGACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATC  >  NC_000913/2415579‑2415873
                                                                                                                                                   |                                                                                                                                                   
tGACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCt                                                                                                                                                     >  3:487126/1‑148 (MQ=255)
 gACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTg                                                                                                                                                    <  4:869490/148‑1 (MQ=255)
  aCTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   <  3:499142/148‑1 (MQ=255)
  aCTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   <  4:487126/148‑1 (MQ=255)
             ccGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATg                                                                                                                                        <  1:544684/148‑1 (MQ=255)
                      cGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTa                                                                                                                                          >  1:680311/1‑137 (MQ=255)
                       gTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacacca                                                                                                                              <  4:372341/148‑1 (MQ=255)
                         gCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   >  1:557927/1‑125 (MQ=255)
                         gCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                   <  2:557927/125‑1 (MQ=255)
                               ggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTgg                                                                                                                      >  4:832859/1‑148 (MQ=255)
                                   gCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                  >  3:831581/1‑148 (MQ=255)
                                         cTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctctc                                                                                                            <  3:832859/148‑1 (MQ=255)
                                          tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              >  1:412669/1‑113 (MQ=255)
                                          tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              <  2:412669/113‑1 (MQ=255)
                                          tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGa                                                                                                                                    <  4:691782/123‑1 (MQ=255)
                                          tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGa                                                                                                                                    >  3:691782/1‑123 (MQ=255)
                                           ggCAGCCATGAACGGCGTAGAAATCGGTGCCCGGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              <  2:775378/112‑1 (MQ=255)
                                                  aTGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              <  1:354550/105‑1 (MQ=255)
                                                  aTGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              >  2:354550/1‑105 (MQ=255)
                                                  aTGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAg                                                                                                        >  3:754342/1‑143 (MQ=255)
                                                  aTGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAg                                                                                                        <  4:754342/143‑1 (MQ=255)
                                                  aTGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGc                                                                                                   >  1:560889/1‑148 (MQ=255)
                                                            tAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              <  4:29868/95‑1 (MQ=255)
                                                            tAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                              >  3:29868/1‑95 (MQ=255)
                                                            tAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATg                                                                                                                                        <  1:996046/101‑1 (MQ=255)
                                                             aGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGt                                                                                                                                               >  1:19616/1‑93 (MQ=255)
                                                             aGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGt                                                                                                                                               <  2:19616/93‑1 (MQ=255)
                                                               aaaTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                         >  4:32647/1‑129 (MQ=255)
                                                               aaaTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                         <  3:32647/129‑1 (MQ=255)
                                                                 aTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaa                                                                                                                                   <  3:196354/101‑1 (MQ=255)
                                                                 aTCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaa                                                                                                                                   >  4:196354/1‑101 (MQ=255)
                                                                  tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTa                                                                                                                                          >  1:752808/1‑93 (MQ=255)
                                                                  tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTa                                                                                                                                          <  2:752808/93‑1 (MQ=255)
                                                                  tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAg                                                                                                        >  1:856327/1‑127 (MQ=255)
                                                                  tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAg                                                                                                        <  2:856327/127‑1 (MQ=255)
                                                                  tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCACGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCt                                                                                                                        >  2:859215/1‑111 (MQ=255)
                                                                       gCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTc                                                                              <  2:560889/148‑1 (MQ=255)
                                                                        cccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                  <  2:274640/111‑1 (MQ=255)
                                                                        cccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                  >  1:274640/1‑111 (MQ=255)
                                                                               gctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGAc                                                                      >  4:372638/1‑148 (MQ=255)
                                                                                   aCTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       >  2:40893/1‑79 (MQ=255)
                                                                                   aCTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                       <  1:40893/79‑1 (MQ=255)
                                                                                         ggTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGa                                                                                   >  4:174453/1‑125 (MQ=255)
                                                                                         ggTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGa                                                                                   <  3:174453/125‑1 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGc                                                                                                   >  2:100148/1‑106 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGc                                                                                                   <  1:100148/106‑1 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTa                                                         <  4:831581/148‑1 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCGGAACTTccagg                                                                          >  4:523699/1‑127 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCGGAACTTccagg                                                                          <  3:523699/130‑4 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCGGAACTTccagg                                                                          <  1:1109120/130‑4 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCGGAACTTccagg                                                                          >  2:1109120/1‑127 (MQ=255)
                                                                                            tACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACAGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                             >  4:598741/1‑144 (MQ=255)
                                                                                              cGAAATGGACGCGCGCATTTCTATACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGAt                                                                   >  1:955952/1‑136 (MQ=255)
                                                                                              cGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacac                                                                                                                                <  1:886968/75‑1 (MQ=255)
                                                                                              cGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACTTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGAt                                                                   <  2:955952/136‑1 (MQ=255)
                                                                                                aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctc                                                                                                              <  1:1090880/91‑1 (MQ=255)
                                                                                                aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctc                                                                                                              >  2:1090880/1‑91 (MQ=255)
                                                                                                aaaTGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCga                                                     >  4:775850/1‑148 (MQ=255)
                                                                                                    ggACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAAcc                                                                                       <  1:191500/110‑1 (MQ=255)
                                                                                                    ggACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAAcc                                                                                       >  2:191500/1‑110 (MQ=255)
                                                                                                           cgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTg                                                                        <  3:982321/118‑1 (MQ=255)
                                                                                                           cgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                             >  3:599516/1‑129 (MQ=255)
                                                                                                           cgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAgg                                          >  2:55203/1‑148 (MQ=255)
                                                                                                             cATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaa                                                  >  4:888670/1‑138 (MQ=255)
                                                                                                             cATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaa                                                  <  3:888670/138‑1 (MQ=255)
                                                                                                              aTTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAAt                                       <  1:55203/148‑1 (MQ=255)
                                                                                                                  cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGt                                   >  2:236756/1‑148 (MQ=255)
                                                                                                                  cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGt                                   <  3:1045286/148‑1 (MQ=255)
                                                                                                                   tAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGc                                                                                              >  1:996807/1‑88 (MQ=255)
                                                                                                                   tAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGc                                                                                              <  2:996807/88‑1 (MQ=255)
                                                                                                                    aaaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                 <  1:735013/132‑1 (MQ=255)
                                                                                                                    aaaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                 >  2:735013/1‑132 (MQ=255)
                                                                                                                     aaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCg                                                                            >  1:888878/1‑104 (MQ=255)
                                                                                                                     aaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCg                                                                            <  2:888878/104‑1 (MQ=255)
                                                                                                                           gCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTAt                                                        >  4:934453/1‑118 (MQ=255)
                                                                                                                           gCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTAt                                                        <  3:934453/118‑1 (MQ=255)
                                                                                                                           gCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAAt                                       <  4:1113387/135‑1 (MQ=255)
                                                                                                                           gCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAAt                                       >  3:1113387/1‑135 (MQ=255)
                                                                                                                              aaCGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGCTCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTaa                             >  4:857564/1‑142 (MQ=255)
                                                                                                                              aaCGTGCTTTCGCTACCGGCCTGCCGGGATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTaa                             <  3:857564/142‑1 (MQ=255)
                                                                                                                               aCGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                             <  4:1139248/109‑1 (MQ=255)
                                                                                                                               aCGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                             >  3:1139248/1‑109 (MQ=255)
                                                                                                                                  tGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATc                                                                  >  2:669432/1‑101 (MQ=255)
                                                                                                                                     ttgcGCTACCCGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGCGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCt                >  1:914813/4‑148 (MQ=255)
                                                                                                                                      ttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAg                                           >  1:285644/1‑120 (MQ=255)
                                                                                                                                      ttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAg                                           <  2:285644/120‑1 (MQ=255)
                                                                                                                                      ttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTg               <  3:372638/148‑1 (MQ=255)
                                                                                                                                       taaatACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACg                                                                     <  2:921385/89‑1 (MQ=255)
                                                                                                                                        aaatACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGt                                                                          <  1:849182/84‑1 (MQ=255)
                                                                                                                                        aaatACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGt                                                                          >  2:849182/4‑87 (MQ=255)
                                                                                                                                           tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGa              <  1:1079600/144‑1 (MQ=255)
                                                                                                                                           tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGa              >  2:1079600/1‑144 (MQ=255)
                                                                                                                                           tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg          >  4:117230/1‑148 (MQ=255)
                                                                                                                                           tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg          <  1:236756/148‑1 (MQ=255)
                                                                                                                                           tACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGAt                                                                   <  3:831429/91‑1 (MQ=255)
                                                                                                                                            aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg          <  2:263620/147‑1 (MQ=255)
                                                                                                                                            aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg          >  1:263620/1‑147 (MQ=255)
                                                                                                                                            aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGa         <  4:1144580/148‑1 (MQ=255)
                                                                                                                                                 ccTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACg                  >  1:828491/1‑134 (MQ=255)
                                                                                                                                                   tGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAtc  <  3:117230/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAATC  >  NC_000913/2415579‑2415873

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