Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,110,400 A→T 5.4% noncoding (35/76 nt) pheV → tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,110,4000AT 5.4% 493.4 / 10.9 161noncoding (35/76 nt)pheVtRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (79/72);  new base T (3/6);  total (82/79)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 165 total aligned reads displayed.

TGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACG  >  NC_000913/3110254‑3110549
                                                                                                                                                  |                                                                                                                                                     
 gCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAAt                                                                                                                                                   >  2:2097671/1‑150 (MQ=255)
  ccTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATc                                                                                                                                                  >  1:854602/1‑150 (MQ=255)
  ccTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATc                                                                                                                                                  >  1:1843786/1‑150 (MQ=255)
  ccTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATc                                                                                                                                                  >  2:1349475/1‑150 (MQ=255)
   cTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcc                                                                                                                                                 <  2:95523/150‑1 (MQ=255)
   cTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATcc                                                                                                                                                 <  2:802827/150‑1 (MQ=255)
    tGACAATGCGTGCAATATCGGCAAAGTGTTGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCTCCGCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATccc                                                                                                                                                <  1:1616385/150‑1 (MQ=255)
      ccAATGCGTGCAATATCGGCAAAGTGCTGATAGATTGTGCAGTCTGCAGGAAATTGAAGAAATTTGCTTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGAGAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCg                                                                                                                                              <  2:1251493/149‑1 (MQ=255)
      aCAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCg                                                                                                                                              <  1:1920400/150‑1 (MQ=255)
       cAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctg                                                                                                                                             <  1:1722017/150‑3 (MQ=255)
         aTGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTTCAATcccctgct                                                                                                                                           <  1:1214889/150‑5 (MQ=255)
          tGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTc                                                                                                                                          >  2:1678517/1‑150 (MQ=255)
             gTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCtt                                                                                                                                       >  1:990854/1‑150 (MQ=255)
               gCAATATCGGCAAAGTGATGATAGATTGTGTAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTgg                                                                                                                                     >  2:1599705/1‑150 (MQ=255)
               gCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTgg                                                                                                                                     <  1:1510644/150‑1 (MQ=255)
               gCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTgg                                                                                                                                     <  2:1090836/150‑1 (MQ=255)
                  atatCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                  <  1:1663948/150‑1 (MQ=255)
                  atatCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                  <  2:1519710/150‑1 (MQ=255)
                  atatCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATATGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTc                                                                                                                                  <  2:205376/150‑1 (MQ=255)
                       ggCAAAGTGATGATAGATGGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTc                                                                                                                             <  1:1693402/150‑1 (MQ=255)
                         cAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCg                                                                                                                           <  2:790285/150‑1 (MQ=255)
                           aaGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAg                                                                                                                         <  1:1183758/150‑1 (MQ=255)
                                 tgatAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGAAAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCggg                                                                                                                   >  2:95063/1‑150 (MQ=255)
                                  tatAGATTGTGCAGTCTGCTGTAAATTGAAGAATTGTGATTGACGAGCCGAGGCGACTCAGGTTTAATGCGCCGCGTTGCCTGGATAGCTCAGTCGGTAGAGCAGGGGATTGAATATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGc                                                                                                                  <  1:2148292/149‑1 (MQ=18)
                                   atAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGAAGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGTATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCACGTGTCCTTGGTTCGATTCCGAGTCAGGGca                                                                                                                 >  1:46253/1‑150 (MQ=255)
                                          gtgCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTaa                                                                                                          >  2:439149/1‑150 (MQ=255)
                                           tgCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATTCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGTTTCGATTCCGAGTCCGGGCACCACTAAt                                                                                                         <  2:37252/150‑1 (MQ=255)
                                              aGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCt                                                                                                      <  1:973541/150‑1 (MQ=255)
                                              aGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTATAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCt                                                                                                      >  2:113160/1‑150 (MQ=255)
                                               gTCTGCAGTAAATTGAAGAAATATGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCAGGGCACCACTAATTCtt                                                                                                     >  1:1372139/1‑150 (MQ=255)
                                                 cTGCAGTAAATTGACGAACTTTGATTGACGACACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCACAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTGCGATTCCGAGTCCGGGCACCACTAATTCTTaa                                                                                                   <  1:1144670/150‑1 (MQ=21)
                                                 cTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGCTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCATaa                                                                                                   <  2:1126578/150‑1 (MQ=255)
                                                  ggCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAg                                                                                                  <  1:649909/149‑1 (MQ=255)
                                                      gTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcc                                                                                              >  1:641382/1‑150 (MQ=255)
                                                       tAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAccc                                                                                             >  1:942946/1‑150 (MQ=255)
                                                        aaaTTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAAcccc                                                                                            >  2:2040555/1‑149 (MQ=255)
                                                            tGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                                                        >  2:625670/1‑150 (MQ=255)
                                                            tGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCca                                                                                        >  2:1994356/1‑150 (MQ=255)
                                                             gaagaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                                                       >  2:1091153/1‑150 (MQ=255)
                                                             gaagaaATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCcac                                                                                       >  2:1843079/1‑150 (MQ=255)
                                                                      tgattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                                                              <  1:595207/150‑1 (MQ=255)
                                                                      tgattgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCAACACTAATTCTTAAGAACCCGCCCACAAGGCGGGt                                                                              >  1:2090646/1‑150 (MQ=255)
                                                                        attgaCGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGttt                                                                            <  2:1499284/150‑1 (MQ=255)
                                                                           gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                         >  2:69983/1‑150 (MQ=255)
                                                                           gacgagacgaGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTg                                                                         >  2:1235827/1‑150 (MQ=255)
                                                                            acgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGc                                                                        >  2:973498/1‑150 (MQ=255)
                                                                             cgagacgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCt                                                                       <  1:1448239/150‑1 (MQ=255)
                                                                                 ccgagGCGAATCAGGTTTACTGCGCGCCGTTGCCCGGATAGCGCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                                                                   <  2:1475206/149‑1 (MQ=255)
                                                                                 acgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTg                                                                   <  1:1897215/150‑1 (MQ=255)
                                                                                  cgagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCAACAAGTCGGGAGTTTGGTTTTgg                                                                  >  2:320020/1‑150 (MQ=18)
                                                                                   gagGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGa                                                                 >  1:1893904/1‑150 (MQ=255)
                                                                                     ggCGACTCAGGTTTAATGCGGCCCGGTGCCCGGATACCTCAGTCGGTAGAGCAGGGGATTGACAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                               <  1:1677209/150‑1 (MQ=21)
                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                               >  1:577278/1‑150 (MQ=255)
                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                               >  1:1636670/1‑150 (MQ=255)
                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                               <  1:1703065/150‑1 (MQ=255)
                                                                                     ggCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATc                                                               >  2:2072491/1‑150 (MQ=255)
                                                                                      gcttaTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCAGAGTCCGGTCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                              <  2:377680/146‑1 (MQ=255)
                                                                                      gCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                              <  2:864247/150‑1 (MQ=255)
                                                                                      gCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCt                                                              <  1:2097671/150‑1 (MQ=255)
                                                                                       ggAATCAGGTTTAATGCTCCCCGTTGCGCGTATAGCTCAGTCGGTAGAGCAGGGGTTTGAAAATCCCCGTGTCGTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                                                             <  2:649725/149‑1 (MQ=21)
                                                                                       cGAATCAGGTTTACTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCATGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTg                                                             <  1:2124700/150‑1 (MQ=255)
                                                                                          aTCAGGTTTAATGTGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                                                          <  2:6112/150‑1 (MQ=255)
                                                                                          aTCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACa                                                          >  1:533616/1‑150 (MQ=255)
                                                                                              ggTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCATTTGGATCTGAAAATaa                                                      >  2:914735/1‑150 (MQ=255)
                                                                                               gTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAAc                                                     >  2:1201829/1‑150 (MQ=255)
                                                                                                   aaTGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTc                                                 <  2:456166/150‑1 (MQ=255)
                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                               >  2:1037445/1‑150 (MQ=255)
                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                               <  1:1128553/150‑1 (MQ=255)
                                                                                                     tGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAc                                               <  1:439149/150‑1 (MQ=255)
                                                                                                          cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCAAGaaaa                                          >  1:410616/1‑150 (MQ=255)
                                                                                                          cccGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACACGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGAAAATAACCTTCAAGaaaa                                          >  2:231790/1‑150 (MQ=255)
                                                                                                               tGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTTATAATCCCCGTATCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGATCCCGCCCACACGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAAtt                                     <  2:1465104/150‑1 (MQ=25)
                                                                                                                  ccGGATTGATCAGTCGTTAGAGCAGTGGTTTGAAAATCCCCGTGTCCTTGTTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCAACAAGGCGGGTTTTTGCTTTTGGATCTTACAATAACCTTCACGAAAAAAATTAGc                                  <  1:1270708/150‑1 (MQ=37)
                                                                                                                    ggATAGGACAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCtt                                <  1:1002225/150‑1 (MQ=255)
                                                                                                                         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTATGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGTTCTGACAATAACCTTCACGAAACAAATTAGCTTATaaa                           >  1:1946270/1‑150 (MQ=255)
                                                                                                                         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATaaa                           >  2:1431261/1‑150 (MQ=255)
                                                                                                                         gCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGCTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATaaa                           >  2:127590/1‑150 (MQ=255)
                                                                                                                          cTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAg                          <  1:1843079/150‑1 (MQ=255)
                                                                                                                               tCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTgg                     <  2:1843786/150‑1 (MQ=255)
                                                                                                                                cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggg                    >  1:588656/1‑150 (MQ=255)
                                                                                                                                cGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggg                    >  1:2059598/1‑150 (MQ=255)
                                                                                                                                 ggTAGAGCAGGTGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGCCTgggg                   <  2:429598/150‑1 (MQ=255)
                                                                                                                                 ggTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTgggg                   <  1:1033070/150‑1 (MQ=255)
                                                                                                                                  gTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTggggg                  <  1:1349475/150‑1 (MQ=255)
                                                                                                                                     gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAt               >  1:908109/1‑150 (MQ=255)
                                                                                                                                     gagCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAt               <  2:1904327/150‑1 (MQ=255)
                                                                                                                                     gacacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAt               >  2:683299/6‑150 (MQ=255)
                                                                                                                                      acacggggATTTTCAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAAtt              >  2:2155430/5‑150 (MQ=255)
                                                                                                                                       gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTa             <  2:1944406/150‑1 (MQ=255)
                                                                                                                                       gCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTa             >  2:1980345/1‑150 (MQ=255)
                                                                                                                                          ggggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTActc          >  1:720445/1‑150 (MQ=255)
                                                                                                                                            ggATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACtctc        >  2:2061586/1‑150 (MQ=255)
                                                                                                                                             gATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCTCAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCg       >  2:1148192/1‑150 (MQ=255)
                                                                                                                                               ttGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCCCCACTAATTCTTAAGAACCCGCCCACAAGGCGGGGTTTTGCTTTTGGAGCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGTGAATTACTCTCgac     >  1:1129739/1‑148 (MQ=255)
                                                                                                                                                tGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGGCCGTGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCa    <  1:148868/150‑1 (MQ=255)
                                                                                                                                                 gAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc   >  2:560201/1‑150 (MQ=255)
                                                                                                                                                 gAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCAc   >  1:906507/1‑150 (MQ=255)
                                                                                                                                                  aaaaTCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACg  >  1:710310/1‑150 (MQ=255)
                                                                                                                                                  aaaaTCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACg  <  1:470519/150‑1 (MQ=255)
                                                                                                                                                  aaaaTCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAAATCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCATATAAAGTCTGGGGGAATTACTCTCGCCACg  >  1:2009569/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGCCTGACAATGCGTGCAATATCGGCAAAGTGATGATAGATTGTGCAGTCTGCAGTAAATTGAAGAAATTTGATTGACGAGACGAGGCGAATCAGGTTTAATGCGCCCCGTTGCCCGGATAGCTCAGTCGGTAGAGCAGGGGATTGAAAATCCCCGTGTCCTTGGTTCGATTCCGAGTCCGGGCACCACTAATTCTTAAGAACCCGCCCACAAGGCGGGTTTTTGCTTTTGGATCTGACAATAACCTTCACGAAAAAAATTAGCTTATAAAGTCTGGGGGAATTACTCTCGCCACG  >  NC_000913/3110254‑3110549

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