Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,415,725 C→A 5.9% L327M (CTG→ATG) ‡ pta → phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,415,7250CA 5.9% 649.2 / 22.5 185L327M (CTG→ATG) ‡ptaphosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (72/102);  new base A (8/3);  total (80/105)
Fisher's exact test for biased strand distribution p-value = 5.87e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.71e-01

Only 100 of 187 total aligned reads displayed.

GACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAT  >  NC_000913/2415580‑2415872
                                                                                                                                                 |                                                                                                                                                   
gACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGc                                                                                                                                                  <  4:155899/149‑1 (MQ=255)
  cTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCg                                                                                                                                                <  2:1408852/149‑1 (MQ=255)
  cTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGCCCTGCCg                                                                                                                                                <  2:807327/149‑1 (MQ=255)
   ttCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCgg                                                                                                                                               >  4:1402862/1‑149 (MQ=255)
     ccGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTa                                                                                                                                             <  3:1311246/149‑1 (MQ=255)
     ccGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCAGGCCGGTa                                                                                                                                             <  4:283051/149‑1 (MQ=255)
       gCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTAtt                                                                                                                                           <  3:1203963/149‑1 (MQ=255)
          gACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTAt                                                                                                                                        <  3:427188/149‑1 (MQ=255)
          gACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTAt                                                                                                                                        <  1:843524/149‑1 (MQ=255)
            ccGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATgg                                                                                                                                      <  1:788456/149‑1 (MQ=255)
              gTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTg                                                                                                                                    <  2:10478/149‑1 (MQ=255)
               tCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGa                                                                                                                                   >  3:492996/1‑149 (MQ=255)
                ccagACGTGCCGGTGGCCGCGTGCCTGGCAGCCATGAACGGCGTAGAACTCGGTGCCCTGCTTCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAGCGTGTTTTCGCTCCCGGCCTGCCGGTATTTATGGTGaa                                                                                                                                  <  2:1445336/146‑1 (MQ=255)
                ccTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaa                                                                                                                                  <  2:247389/149‑1 (MQ=255)
                  tGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaaca                                                                                                                                >  3:1422286/1‑149 (MQ=255)
                  tGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaaca                                                                                                                                <  1:1289872/149‑1 (MQ=255)
                    aCGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacacc                                                                                                                              >  4:1316217/1‑149 (MQ=255)
                      gTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacaccaa                                                                                                                            <  4:257200/149‑1 (MQ=255)
                           ggtggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCt                                                                                                                       <  1:95062/149‑1 (MQ=255)
                              ggCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGc                                                                                                                    >  4:1248396/1‑149 (MQ=255)
                                 cGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGa                                                                                                                 >  2:1130387/1‑149 (MQ=255)
                                  gCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGAc                                                                                                                <  2:811293/149‑1 (MQ=255)
                                  gCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGAc                                                                                                                >  1:1007443/1‑149 (MQ=255)
                                    ttGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACct                                                                                                              <  3:46937/149‑1 (MQ=255)
                                        cTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACctctct                                                                                                          >  1:1414840/1‑149 (MQ=255)
                                         tGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTg                                                                                                         >  1:337403/1‑149 (MQ=255)
                                          ggCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacac                                                                                                                               <  4:715878/126‑1 (MQ=255)
                                          ggCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGaacac                                                                                                                               >  3:715878/1‑126 (MQ=255)
                                           gCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAg                                                                                                       >  1:39526/1‑149 (MQ=255)
                                              gCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCt                                                                                                    <  4:396217/149‑1 (MQ=255)
                                                   gAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCaga                                                                                               >  4:1377276/1‑149 (MQ=255)
                                                     aCGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGc                                                                                             <  4:1456582/149‑1 (MQ=255)
                                                      cGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                            <  4:1196486/149‑1 (MQ=255)
                                                      cGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCt                                                                                            <  3:650658/149‑1 (MQ=255)
                                                          gTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCaa                                                                                        >  2:1363930/1‑149 (MQ=255)
                                                            aGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAAcc                                                                                      >  2:633277/1‑149 (MQ=255)
                                                                aTCGGTGCCCTGCTGCTGACTTGCGGTTACGAAATTGACGCGCGCATTTCTAAACTGTGCGCACGTGCCTTCACTACCGGCCTGCCGGTATTTATGGTGAAGACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGa                                                                                  >  2:492903/1‑149 (MQ=255)
                                                                 tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCACCCTGGaa                                                                                 <  4:736698/149‑1 (MQ=255)
                                                                 tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGaa                                                                                 <  1:633277/149‑1 (MQ=255)
                                                                 tCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGaa                                                                                 <  2:30177/149‑1 (MQ=255)
                                                                  cGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAg                                                                                <  4:708865/149‑1 (MQ=255)
                                                                     tGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTc                                                                             <  4:1422286/149‑1 (MQ=255)
                                                                     tGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTc                                                                             <  3:152116/149‑1 (MQ=255)
                                                                      gCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTcc                                                                            <  4:955969/149‑1 (MQ=255)
                                                                       cccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCg                                                                           >  3:965329/1‑149 (MQ=255)
                                                                       cccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCg                                                                           <  1:409526/149‑1 (MQ=255)
                                                                        ccTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCgg                                                                          >  4:78660/1‑149 (MQ=255)
                                                                         ctgctgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTcc                                                                            >  2:438811/1‑146 (MQ=255)
                                                                         ctgctgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGt                                                                         >  4:302159/1‑149 (MQ=255)
                                                                          tgctgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGtt                                                                        >  1:1031949/1‑149 (MQ=255)
                                                                             tgctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGCTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGAc                                                                     <  3:898506/149‑1 (MQ=255)
                                                                              gctgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGa                                                                                                        <  2:1238035/113‑1 (MQ=255)
                                                                                tgACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTTTGAGCCTGCAGAGCTACAACCTGGAAGTTCCGGTTGACGAt                                                                  <  3:951514/149‑1 (MQ=255)
                                                                                 gACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATc                                                                 <  3:1000966/149‑1 (MQ=255)
                                                                                     ggCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGa                                                             <  4:468418/149‑1 (MQ=255)
                                                                                      gCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                            <  4:967083/149‑1 (MQ=255)
                                                                                      gCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGaa                                                            >  2:1349657/1‑149 (MQ=255)
                                                                                       cGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCAGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAAc                                                           <  3:302159/149‑1 (MQ=255)
                                                                                         gTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGt                                                         >  1:170978/1‑149 (MQ=255)
                                                                                          ttACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTg                                                                                    >  2:35664/1‑121 (MQ=255)
                                                                                            aCGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGt                                                                         >  1:748823/1‑130 (MQ=255)
                                                                                              gAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGATTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTTCAGAGCTTCAACATGGAAGTTCCGGTTGACGATCACGAACGTATCga                                                    >  1:1213366/1‑149 (MQ=255)
                                                                                                 aTGGACGCGCGCATTTTTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaa                                                 >  3:1292846/1‑149 (MQ=255)
                                                                                                  tGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                <  2:72108/149‑1 (MQ=255)
                                                                                                  tGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                <  2:224390/149‑1 (MQ=255)
                                                                                                   ggACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAg                                               >  3:815290/1‑149 (MQ=255)
                                                                                                     aCGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGa                                                                      <  3:224442/124‑1 (MQ=255)
                                                                                                       gcgcgcATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                           >  4:512669/1‑149 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                 <  3:134946/149‑1 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                 <  4:800806/149‑1 (MQ=255)
                                                                                                                 cTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGtt                                 <  4:1002865/149‑1 (MQ=255)
                                                                                                                  tAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTg                                >  3:1443197/1‑149 (MQ=255)
                                                                                                                   aaaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                >  3:1353630/1‑132 (MQ=255)
                                                                                                                   aaaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGaaa                                                <  4:1353630/132‑1 (MQ=255)
                                                                                                                    aaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCt                              >  2:1259860/1‑149 (MQ=255)
                                                                                                                    aaCTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCt                              >  4:119311/1‑149 (MQ=255)
                                                                                                                       tgtgCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTTAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGGATACGTTGCTAAc                           >  4:641590/1‑149 (MQ=255)
                                                                                                                       tgtgCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAAc                           <  1:1051006/149‑1 (MQ=255)
                                                                                                                         tgCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGtt                                             <  2:1037379/129‑1 (MQ=255)
                                                                                                                           cGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACa                       <  3:934956/149‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATAc                                    <  4:494969/128‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAAc                  <  2:769161/146‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAAc                  >  1:769161/1‑146 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCt               <  4:832189/149‑1 (MQ=255)
                                                                                                                                   cTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCt               <  3:497986/149‑1 (MQ=255)
                                                                                                                                    tttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTg              >  2:270299/1‑149 (MQ=255)
                                                                                                                                     ttCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACCCCAACGCTGa             <  2:1069197/149‑1 (MQ=255)
                                                                                                                                      tCGCTACTGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGAc            <  4:909630/149‑1 (MQ=255)
                                                                                                                                       aaatACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCACCCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCa                                           <  3:1047876/114‑1 (MQ=255)
                                                                                                                                         atACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg         >  1:1472988/2‑149 (MQ=255)
                                                                                                                                         atACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg         >  3:528622/2‑149 (MQ=255)
                                                                                                                                         atACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg         >  2:962666/2‑149 (MQ=255)
                                                                                                                                         atACCGGCAGGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTgg         >  2:1475383/2‑149 (MQ=255)
                                                                                                                                          tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGGGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACTTCAACGCTGACTgca        >  2:430609/1‑147 (MQ=255)
                                                                                                                                          tACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGa        >  1:444115/1‑149 (MQ=255)
                                                                                                                                           aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGAt       <  2:1302181/149‑1 (MQ=255)
                                                                                                                                           aCCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCACTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGAt       <  3:119311/149‑1 (MQ=255)
                                                                                                                                              ggCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGa    >  1:242473/1‑149 (MQ=255)
                                                                                                                                               gCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGAt       <  2:556380/145‑1 (MQ=255)
                                                                                                                                               gCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGaa   <  1:519019/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACTTCCGCAGACCGTCCTGACGTGCTGGTGGCCGCTTGCCTGGCAGCCATGAACGGCGTAGAAATCGGTGCCCTGCTGCTGACTGGCGGTTACGAAATGGACGCGCGCATTTCTAAACTGTGCGAACGTGCTTTCGCTACCGGCCTGCCGGTATTTATGGTGAACACCAACACCTGGCAGACCTCTCTGAGCCTGCAGAGCTTCAACCTGGAAGTTCCGGTTGACGATCACGAACGTATCGAGAAAGTTCAGGAATACGTTGCTAACTACATCAACGCTGACTGGATCGAAT  >  NC_000913/2415580‑2415872

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