Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,448 G→A 5.5% A33V (GCC→GTC)  yfeX ← porphyrinogen oxidase, cytoplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,4480GA 5.5% 236.1 / 2.3 73A33V (GCC→GTC) yfeXporphyrinogen oxidase, cytoplasmic
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (2/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

GTAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAG  >  NC_000913/2550300‑2550593
                                                                                                                                                    |                                                                                                                                                 
gTAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTgg                                                                                                                                                   <  7:155048/149‑1 (MQ=255)
gTAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTgg                                                                                                                                                   <  3:208156/149‑1 (MQ=255)
gTAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTgg                                                                                                                                                   <  4:343277/149‑1 (MQ=255)
 tAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGc                                                                                                                                                  <  1:399224/149‑1 (MQ=255)
    cccGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCcgc                                                                                                                                               >  1:97512/1‑149 (MQ=255)
     ccGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCcgca                                                                                                                                              >  1:197456/1‑149 (MQ=255)
      cGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGc                                                                                                                                                  <  2:303837/144‑1 (MQ=255)
      cGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGc                                                                                                                                                  >  1:303837/1‑144 (MQ=255)
      cGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCcgcac                                                                                                                                             >  2:179365/1‑149 (MQ=255)
           aaTCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAgg                                                                                                                                        <  8:84999/149‑1 (MQ=255)
                 tCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGc                                                                                                                                      <  2:413224/145‑1 (MQ=255)
                 tCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGc                                                                                                                                      <  5:53667/145‑1 (MQ=255)
                 tCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTc                                                                                                                                   <  1:33798/148‑1 (MQ=255)
                 tCAGATCTTCTGCGCCAACGACGCAGCTCAGAGCTCGCCAGGTGTTGTTAACAAAGGCAACCACCGCACCAAGATGGGCGTCCGGTAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCTCGCAGGGc                                                                                                                                      <  6:401753/145‑1 (MQ=255)
                   aGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCACCCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCacc                                                                                                                                >  4:201663/1‑147 (MQ=255)
                        ttctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTccc                                                                                                                           >  8:11351/1‑149 (MQ=255)
                        ttctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTc                                                                                                                             <  2:97512/147‑1 (MQ=255)
                          ctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTcccc                                                                                                                          <  5:208926/148‑1 (MQ=255)
                          ctGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTcc                                                                                                                            <  7:82896/146‑1 (MQ=255)
                              cccAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCa                                                                                                                     >  2:43542/1‑149 (MQ=255)
                                  aCGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGtt                                                                                                                 >  6:113557/1‑149 (MQ=255)
                                              agagCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGc                                                                                                                                                  <  5:99939/104‑1 (MQ=255)
                                              agagCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGc                                                                                                                                                  >  6:99939/1‑104 (MQ=255)
                                               gagCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCaa                                                                                                    <  2:15734/149‑1 (MQ=255)
                                                    cgcCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATcgc                                                                                               <  5:113557/149‑1 (MQ=255)
                                                     gcCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGACCGGGAATTTCGCTTCAAAAGTTCCCAGATTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCATTCACGTTCCCTTCCAACCAAATcgcc                                                                                              >  4:173455/1‑149 (MQ=255)
                                                       cAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATcgccgc                                                                                            >  7:230461/1‑149 (MQ=255)
                                                        aGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCcgcg                                                                                           >  8:268662/1‑149 (MQ=255)
                                                         ggTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCcgcgc                                                                                          <  1:43542/149‑1 (MQ=255)
                                                           tgttgtTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCgg                                                                                        <  5:120603/149‑1 (MQ=255)
                                                                  aCCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGt                                                                                  >  8:15585/1‑148 (MQ=255)
                                                                   ccAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTc                                                                                >  3:347597/1‑149 (MQ=255)
                                                                   ccAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGGCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGtt                                                                                 <  7:312705/148‑1 (MQ=255)
                                                                     aaaGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTg                                                                              <  8:151/149‑1 (MQ=255)
                                                                           aaCCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCaaaa                                                                        >  6:35663/1‑149 (MQ=255)
                                                                            accaccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAAt                                                                       <  1:174685/149‑1 (MQ=255)
                                                                             ccaccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGAt                                                                                                        <  6:278248/115‑1 (MQ=255)
                                                                             ccaccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGAt                                                                                                        >  5:278248/1‑115 (MQ=255)
                                                                              caccGCACCAAGATGCGGGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGc                                                                     >  6:349110/1‑149 (MQ=255)
                                                                              caccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGc                                                                     >  8:399012/1‑149 (MQ=255)
                                                                              caccGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGc                                                                     <  1:179365/149‑1 (MQ=255)
                                                                                    aCCAAGATGCGGGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCActct                                                               <  5:349110/149‑1 (MQ=255)
                                                                                          aTGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCt                                                         >  2:271232/1‑149 (MQ=255)
                                                                                            gcgcGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGAAACTGTTTTACTGGCCGCACGCAGGGCGTAAACTTCCCCTTTCACGTTGGCGTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCCAAATGCCACTCTGAACCt                                                         >  7:286959/1‑147 (MQ=255)
                                                                                             cgcgTCCGGGAATTTCGCTTCAATTGTTGGGAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTga                                                       <  3:173455/148‑1 (MQ=255)
                                                                                              gcgTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCTCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTgag                                                      >  4:385314/1‑148 (MQ=255)
                                                                                                 tCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACAt                                                  <  5:35663/149‑1 (MQ=255)
                                                                                                   cGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCa                                                                   <  3:350884/130‑1 (MQ=255)
                                                                                                   cGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCa                                                                   >  4:350884/1‑130 (MQ=255)
                                                                                                   cGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTg                                                <  7:268662/149‑1 (MQ=255)
                                                                                                    gggAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGAAAAATGCCACTCTGAACCTGAGACATTGt                                               >  7:291830/1‑149 (MQ=255)
                                                                                                      gAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTtc                                             <  2:144553/149‑1 (MQ=255)
                                                                                                         tttCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTtcctc                                          <  7:15585/149‑1 (MQ=255)
                                                                                                              cTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTATGGCAAAATGCCCCTCTGAACCTGAGACATTGTTCCTCATGa                                      >  6:117017/1‑148 (MQ=255)
                                                                                                                tCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCATGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaaaa                                   >  7:62741/1‑149 (MQ=255)
                                                                                                                  aaaaGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaat                                  >  7:390232/1‑148 (MQ=255)
                                                                                                                   aaaGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGTCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataa                                <  4:49982/149‑1 (MQ=255)
                                                                                                                   aaaGTTGCCAGTCTGTCGTCAAATGTTTTACTGGCCGCACGCAGGGCGGCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGAAATGTTCTTGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataa                                <  8:286959/149‑1 (MQ=255)
                                                                                                                    aaGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataa                                >  3:120279/1‑148 (MQ=255)
                                                                                                                      gTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTg                                                >  4:80593/1‑130 (MQ=255)
                                                                                                                      gTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTg                                                <  3:80593/130‑1 (MQ=255)
                                                                                                                      gTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataataa                             <  1:18466/149‑1 (MQ=255)
                                                                                                                       ttGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAaataataat                            <  8:390232/149‑1 (MQ=255)
                                                                                                                       ttGCCAGTTTGTCGGCAAATGTTTTACTGACCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaaaa                                   >  5:197431/1‑142 (MQ=255)
                                                                                                                       ttGCCAGTTTGTCGGCAAATGTTTTACTGACCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaaaa                                   <  6:191518/142‑1 (MQ=255)
                                                                                                                       ttGCCAGTTTGTCGGCAAATGTTTTACTGACCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaaaa                                   >  5:191518/1‑142 (MQ=255)
                                                                                                                       ttGCCAGTTTGTCGGCAAATGTTTTACTGACCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGaaaa                                   <  6:197431/142‑1 (MQ=255)
                                                                                                                            aGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCAc                       <  1:271232/149‑1 (MQ=255)
                                                                                                                              tttGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACg                      >  6:112040/1‑148 (MQ=255)
                                                                                                                              tttGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGc                     <  4:347597/149‑1 (MQ=255)
                                                                                                                                     gCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGt              >  5:305122/1‑149 (MQ=255)
                                                                                                                                       aaaTGTTTTACTGTCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGAGCGGCAATGTTCTGGCAAAATGCCACTCTGAACATGAGACATTGTTCCTCCTGAAAATAATAATGACACGCAATCCGTgg            >  1:24274/1‑149 (MQ=255)
                                                                                                                                       aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg            <  8:230461/149‑1 (MQ=255)
                                                                                                                                          tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATc                                                                                                       >  7:329255/1‑55 (MQ=255)
                                                                                                                                          tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATc                                                                                                       <  8:329255/55‑1 (MQ=255)
                                                                                                                                                 cTGGCCGCACGCAGGGCGTCAACTTACCCTTTCACGTTGGGGTCGATACCAATCGCCGCGCGGGAATGTTATGGCAAAATGCCACTCTGAACCAGAGACATTGGTCCTCCTGAAAATAATAATTCCACGCAATCCGGGGCATTGTTGAg  >  2:255593/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTAGCCCGGAAAATCTTTCAGCTCTTCTGCCCCAACGCCGCCGCTCAGAGCGCGCCAGGTGTTGTTACCAAAGGCAACCACCGCACCAAGATGCGCGTCCGGGAATTTCGCTTCAAAAGTTGCCAGTTTGTCGGCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAG  >  NC_000913/2550300‑2550593

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