Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,186,448 A→C K395Q (AAA→CAA)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,186,4480AC98.5% 209.1 / NA 66K395Q (AAA→CAA) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (38/27);  total (39/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGA  >  NC_000913_3_pae_tpiA/4186303‑4186593
                                                                                                                                                 |                                                                                                                                                 
tgtTGACTACTCCGGTCGTTCTGTAATCAACGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCACACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAg                                                                                                                                                <  1:220366/149‑1 (MQ=255)
    aactactACGGTCGTTATGTAATCACCGTAGGTCCATATCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGc                                                                                                                                            <  1:215760/148‑1 (MQ=255)
        actCCGGTCGTTCTGTAAGCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCgaag                                                                                                                                        <  1:7317/149‑1 (MQ=255)
                 gTTTTGTTATCACCGTCGGTCCATACCTGCGTCTGCACCAGTGCGGTCTGCCGAAGAAAATGGCACTGTAGCTGTTCAAACCGTTCTTCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGtt                                                                                                                               <  2:981824/149‑1 (MQ=255)
                 gTTCTGTAATCACCGTAGGTCCATCCCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGtt                                                                                                                               <  1:757598/149‑1 (MQ=255)
                 gTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGtt                                                                                                                               <  1:539810/149‑1 (MQ=255)
                 gTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGACACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGtt                                                                                                                               >  1:265375/1‑149 (MQ=255)
                    cTGTGATCACCGTAGGTCCATTCCGGCGTCAGCAACAGTGCGGTAGCCGGAAGAACATAGCACTGGAGCTGTTACAACCGTTCATCTACGGCAAGCTGGACCTGCGTGGTCTTACTACCACCATCAAAGCGCAGAAGGAAATggcggtg                                                                                                                            >  2:502408/1‑144 (MQ=255)
                     tGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAgc                                                                                                                           >  2:910174/1‑149 (MQ=255)
                      gTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCTTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAGTGGTTGAgcg                                                                                                                          >  2:239211/1‑149 (MQ=255)
                       tAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAgcgc                                                                                                                         >  1:473037/1‑149 (MQ=255)
                         aTCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCga                                                                                                                       >  2:152430/1‑149 (MQ=255)
                            aCCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCGGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACAACCATTCAAGCTGCGAAGAAAATGGTTGAACGCgaaga                                                                                                                    >  2:825070/1‑149 (MQ=255)
                             ccGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTCTTCAAACCGTTTATGTACGGCAAGCTGGAACTACGTGGTCTTGCTAACACCATTCAAGCTGCCATGAAATTAATGGAGTGGgaagaa                                                                                                                   >  2:322660/1‑149 (MQ=255)
                                       aTACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAAGTTCGTGGTCTTGCTACTACCATTCAAGCTGCGAAAAAAATGGTTGAGCGTGAAGAAGCTGTCGttt                                                                                                         >  1:537075/1‑149 (MQ=255)
                                       aTACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCtatagtat                                                                                                         >  2:364916/1‑142 (MQ=255)
                                       aTACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGttt                                                                                                         <  2:473037/149‑1 (MQ=255)
                                        tACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTg                                                                                                        <  2:722834/149‑1 (MQ=255)
                                          cccgcctCTGCACCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCCGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTggg                                                                                                      <  2:555300/143‑1 (MQ=255)
                                          ccTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTggg                                                                                                      >  2:275594/1‑149 (MQ=255)
                                           cTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCCCTGGAGCTGTTCAAACCGATCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGAGAAGAAAATGGTTGAGCGGAAAGACGCTGTCGTTTggaa                                                                                                     >  1:25853/1‑147 (MQ=255)
                                             gCGTATACATAAGTGAGGTCTGCAGAAGAAAATGGCACTGGAGCTGTTAAAACAGTTCATCTACAGCAAACTGGCAATCCGGGGTCTTGCTACCACCATTCCAGCTGCGAAGAAAAGGGTTGAGCGCGAACAGGACGTTGTTTggagta                                                                                                   >  2:15520/1‑145 (MQ=255)
                                                  tGCATCAGGGCGGTCTCCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTg                                                                                              <  1:445869/149‑1 (MQ=255)
                                                    cATCAGTGCGGTCTGCGGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGa                                                                                            <  1:1041019/149‑1 (MQ=255)
                                                         gTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAg                                                                                       <  1:275594/149‑1 (MQ=255)
                                                                   ccGAAGAAATTGGCACTGTAGCTGTTCAAACCGTTCATCTACGTCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGa                                                                             <  1:799484/149‑1 (MQ=255)
                                                                     gaagaaAATGGCACTGGAGCTGTTCAAACCGTTCATATACAGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGAGATCCCGGACGAAGTTATCCGCGAac                                                                           >  2:321029/1‑149 (MQ=255)
                                                                      aagaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAaca                                                                          >  2:595521/1‑149 (MQ=255)
                                                                        gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTg                                                                                              >  1:409951/1‑127 (MQ=255)
                                                                        gaaAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTg                                                                                              <  2:409951/127‑1 (MQ=255)
                                                                           aaTGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCgg                                                                     >  1:726897/1‑149 (MQ=255)
                                                                               gCACCGGTGCTTTTCAATCCGTTTTTCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGCCGAAGTTATCCGCGAACACCCGGTAct                                                                 <  1:358351/149‑1 (MQ=255)
                                                                                 aCTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAAAACCCGGTActgc                                                               >  2:147721/1‑149 (MQ=255)
                                                                                     gAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAATTATCAGCGAACACACGGTACTGCTGaa                                                           >  1:370737/1‑149 (MQ=255)
                                                                                         tGTTTAATCCGTTCATCTACGGCAAGTTGGATCTGCGTGGTGTTGCGACCATCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGt                                                       <  1:454459/149‑1 (MQ=255)
                                                                                          gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGTTATCCTGGTTGAAGTTTTCTGCGAAAAACCGGTACTGCTGAACtgtg                                                      >  1:125966/1‑149 (MQ=255)
                                                                                          gTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGTACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTg                                                      >  1:415724/1‑149 (MQ=255)
                                                                                           ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGCAGCTGTCGTTTGGGATATCCTGGACGAAATTATCCGCGAACACCCGGTACTGCTGAACCGTGc                                                     >  1:930167/1‑149 (MQ=255)
                                                                                           ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGc                                                     >  2:336448/1‑149 (MQ=255)
                                                                                           ttCAAACCGTTCATCTACGGCAAGCTGGAACTGCGCGGTCTTGCTACCAAACCTCAAGCCGAGAACAACATGGTTGTGCGGACCGAACCCGACGCTTTGGTTCGCCTTAACGAAGGTTTCCACGAACACCCGGTACTGCTGAACCGTGc                                                     >  2:770103/1‑149 (MQ=255)
                                                                                             cAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAc                                                   >  1:1002931/1‑149 (MQ=255)
                                                                                              aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                  >  1:658963/1‑149 (MQ=255)
                                                                                              aaaCCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCAcc                                                  >  1:626682/1‑149 (MQ=255)
                                                                                                aCCGTTCATCTACGGCTAGCTGGAATTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGa                                                <  1:379308/149‑1 (MQ=255)
                                                                                                    gtCACCTACGGCCAGCTGGAACTGGCTACCCCTGCTACTGCCATTCAAGCTAGAATGAAGATGGTTGAGGTTGAAGAAGCTTGCGTTTGGGATATCCTGGTCGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGActct                                            <  2:989095/148‑1 (MQ=255)
                                                                                                          tACGGCAAGCTTGAACTGCGTGGTCTTGCTAACACCATACAAGATGCGAATAAAATGGTTAGGCGGGAAGAAGGTGTCGTTTGGGATATCCAGGGCGAAGGTATCCGCAAAGAACCGGGATTGATGGCCGTTGCACCGACGCCGCACCg                                      >  2:315038/1‑149 (MQ=255)
                                                                                                             ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAct                                              >  2:26834/1‑138 (MQ=255)
                                                                                                             ggCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAct                                              <  1:26834/138‑1 (MQ=255)
                                                                                                              cagaGCGGGTACTGCGTGTTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTg                                  <  2:822277/146‑1 (MQ=255)
                                                                                                               cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                 <  1:1012740/149‑1 (MQ=255)
                                                                                                               cAAGCTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACACGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTgg                                 >  1:824224/1‑149 (MQ=255)
                                                                                                                   cTGGAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCCCCGAATCTGCACCATCTGGGTAt                             >  1:397525/1‑149 (MQ=255)
                                                                                                                    tGGAACTGCGTGGTCTTGTTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATc                            <  1:712376/149‑1 (MQ=255)
                                                                                                                      gAACTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCa                          <  2:839902/149‑1 (MQ=255)
                                                                                                                       aaCTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAACATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                         >  1:276316/1‑149 (MQ=255)
                                                                                                                       aaCTGCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAg                         <  2:658963/149‑1 (MQ=255)
                                                                                                                           gCGTGGTCTTGCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCAt                     <  1:595521/149‑1 (MQ=255)
                                                                                                                                gTCTTGCTACCACCATTCAAGATGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGaa                >  2:722732/1‑149 (MQ=255)
                                                                                                                                   ttGCTACCACCATTCACGCTGCGAGGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATACTGGACGCAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGAGTCTGCACCGTATGGGTATCCAGGCATTTGACCCg             >  2:907741/1‑149 (MQ=255)
                                                                                                                                     gCTACCACCATTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGt           <  1:152430/149‑1 (MQ=255)
                                                                                                                                       tgacaccaTTCAAGTCGCGAAGAAAATGGTTGAGGGGGAAGAAGCTGGCGTTTGGGATATCCTGGACGAAGTTATTCTCAAACCCCCGGCACTGTTTAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTAc         <  2:432617/146‑1 (MQ=255)
                                                                                                                                          caccaTTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTGCTGa      >  2:457680/1‑149 (MQ=255)
                                                                                                                                          caccaTTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGa      >  1:535586/1‑149 (MQ=255)
                                                                                                                                          caccaTTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGAGTATCCAGGCATTTGAACCGGTACTGa      >  1:570781/1‑149 (MQ=255)
                                                                                                                                              aTTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGa  >  1:756214/1‑149 (MQ=255)
                                                                                                                                              aTTCAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATACTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGCATCCAGGCATTGGAAACAGTAGTGATCGa  >  2:499916/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TGTTGACTACTCCGGTCGTTCTGTAATCACCGTAGGTCCATACCTGCGTCTGCATCAGTGCGGTCTGCCGAAGAAAATGGCACTGGAGCTGTTCAAACCGTTCATCTACGGCAAGCTGGAACTGCGTGGTCTTGCTACCACCATTAAAGCTGCGAAGAAAATGGTTGAGCGCGAAGAAGCTGTCGTTTGGGATATCCTGGACGAAGTTATCCGCGAACACCCGGTACTGCTGAACCGTGCACCGACTCTGCACCGTCTGGGTATCCAGGCATTTGAACCGGTACTGATCGA  >  NC_000913_3_pae_tpiA/4186303‑4186593

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑