Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,632,774 C→G 5.9% intergenic (‑184/+78) yhiJ ← / ← yhiL DUF4049 family protein/pseudogene, DUF4049 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,632,7740CG 5.9% 449.8 / 14.4 135intergenic (‑184/+78)yhiJ/yhiLDUF4049 family protein/pseudogene, DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base C (61/64);  new base G (4/4);  total (67/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 135 total aligned reads displayed.

GTTAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCA  >  NC_000913/3632627‑3632922
                                                                                                                                                   |                                                                                                                                                    
gTTAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACAt                                                                                                                                                    >  2:400555/1‑150 (MQ=255)
 ttAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATg                                                                                                                                                   <  1:1639061/150‑1 (MQ=255)
  tAACATTTTATCTTATGACATTCACCTCATGTTCGATGGGAAAATAATTATATTTTTCACAAATAATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCATAGAAAGTCACTAAAGATTTTATAGGCGTAATCAGAAAAtta                                                                                                                                                  >  1:35268/1‑148 (MQ=255)
  tAACATTTTATCTTATGACATACCCCTCATGTGCGCTGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGc                                                                                                                                                  >  1:335024/1‑150 (MQ=255)
  tAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGc                                                                                                                                                  >  2:686213/1‑150 (MQ=255)
  tAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGc                                                                                                                                                  >  2:192996/1‑150 (MQ=255)
   aaCATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGcc                                                                                                                                                 <  1:973324/150‑1 (MQ=255)
    aCATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCct                                                                                                                                                >  2:1420109/1‑150 (MQ=255)
     cATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCctc                                                                                                                                               <  2:1433689/150‑1 (MQ=255)
       ttatagcTTATGACATACACCCCATGATCGATGGGAAAATAATTATATTTTGCCCAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCATAGAAAGTCCCTAAATATTGTAGATTCGTCAGCATCACATGCCTCTa                                                                                                                                             >  2:1130708/7‑150 (MQ=255)
        tttATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCATGTTCTGATGACGCCTCTAc                                                                                                                                            <  1:580458/150‑1 (MQ=255)
         ttATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCAACa                                                                                                                                           >  2:1445739/1‑150 (MQ=255)
          tATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCATGTTCTGATGACGCCTCTACaa                                                                                                                                          <  2:1576861/150‑1 (MQ=255)
            tCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATc                                                                                                                                        <  2:551836/150‑1 (MQ=255)
               tATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGAt                                                                                                                                     >  1:966446/1‑150 (MQ=255)
                        cacCTCATGTTCGATGTGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAATCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAAc                                                                                                                            <  2:1663369/150‑1 (MQ=255)
                        cacCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAAc                                                                                                                            >  1:968549/1‑150 (MQ=255)
                          ccTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCt                                                                                                                          >  2:1221505/1‑150 (MQ=255)
                               tGTTCGATGGGAAAATAATTATATTTTGAACAAATCATTCAGCCTCAAAGCATAATATGAAATTTATTATCTAATAGCAGAGAAAGTCCCTAAAGATTGTATAGGCGTCATCAGAACATGCCTCTACAATAGATTAAGAAAACCTATGaa                                                                                                                     >  1:331801/1‑150 (MQ=255)
                                   cGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATAc                                                                                                                 >  1:1292406/1‑150 (MQ=255)
                                       gggAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGACATACTGGc                                                                                                             <  1:1076056/150‑1 (MQ=255)
                                        ggAAAATAATATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAg                                                                                                           >  1:372160/1‑150 (MQ=255)
                                         gAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAg                                                                                                           >  1:1027405/1‑150 (MQ=255)
                                         gAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAg                                                                                                           <  1:1420109/150‑1 (MQ=255)
                                            aataatTATATTTTTCACAAATCATTCAGGCTAAAAGCATAATATGAAATTTCTTATCTCATATTAGTGAAAGTTCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGAt                                                                                                        <  2:3228/150‑1 (MQ=255)
                                            aataatTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATATCAGAGAAAGTCCCTAAAGATTTTAGAGGCGTCATCAGAACATGCCTCTACAATAGATTAAGATAACCAATGAAATACTGTCAGGAt                                                                                                        >  1:299587/1‑150 (MQ=255)
                                               aatTATATTTTGCACAAATCATTCACCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaa                                                                                                     <  1:272578/150‑1 (MQ=255)
                                                atTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaaa                                                                                                    <  2:384457/150‑1 (MQ=255)
                                                atTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaaa                                                                                                    <  2:646660/150‑1 (MQ=255)
                                                 ttATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaaaa                                                                                                   >  2:1345341/1‑150 (MQ=255)
                                                 ttATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaaaa                                                                                                   >  2:1352100/1‑150 (MQ=255)
                                                  tataTTTTGCACAAATCATGCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTaaaaa                                                                                                  <  1:471636/150‑1 (MQ=255)
                                                    tatTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGa                                                                                                <  1:1653844/150‑1 (MQ=255)
                                                    tatTTTGCACAAATCATTCAGCATCAAAGCATAATATGATATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCTTCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGa                                                                                                <  2:331801/150‑1 (MQ=255)
                                                        ttGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATa                                                                                            <  1:1097306/150‑1 (MQ=255)
                                                          gCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATacc                                                                                          <  1:1221505/150‑1 (MQ=255)
                                                           cacaAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCAAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGAAGGATTAAAAAGACATacca                                                                                         >  1:794360/1‑150 (MQ=255)
                                                            acaAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATaccac                                                                                        >  1:1292376/1‑150 (MQ=255)
                                                             caAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATaccacc                                                                                       >  2:1000672/1‑150 (MQ=255)
                                                              aaaTCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTa                                                                                                      >  2:739453/1‑134 (MQ=255)
                                                              aaaTCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTa                                                                                                      <  1:739453/134‑1 (MQ=255)
                                                                       aGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATgg                                                                             <  2:394406/150‑1 (MQ=255)
                                                                        gCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGTTTGTAGAGGCGTCATCAGAACATGCCTCTACAAACGTTTAAGAAAACCTATGAAATACTGGAATGATTAAAAAGACATACCACAGTATTCATggt                                                                            >  2:1271734/1‑149 (MQ=255)
                                                                        gCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGa                                                                            >  1:1349486/1‑150 (MQ=255)
                                                                           tCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGtt                                                                         >  1:1238330/1‑150 (MQ=255)
                                                                                gCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTGTACAATCGATTAAGAAAACCTATGAAATCCGGGCAGGATTAAAACGACATACCACCGTATTGAGGGAGGTTTAAt                                                                    >  2:829551/1‑150 (MQ=255)
                                                                                  ataataTGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGt                                                                  <  2:335024/150‑1 (MQ=255)
                                                                                  ataataTGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGt                                                                  <  2:1292406/150‑1 (MQ=255)
                                                                                   taataTGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTa                                                                 >  2:157065/1‑150 (MQ=255)
                                                                                    aataTGAAATTTCTTATCTCATAGCAGAGAAAGTCACTAAAGATTGTAGAGGCTTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAg                                                                <  1:1529509/150‑1 (MQ=255)
                                                                                      tatGAAATTTCTTATCTCATTGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAAAATGCCTCTACAATCGAATAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGtt                                                              >  2:539688/1‑150 (MQ=255)
                                                                                       atGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTg                                                             >  2:103627/1‑150 (MQ=255)
                                                                                         taaaTTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGaa                                                           <  2:1272227/149‑1 (MQ=255)
                                                                                             tttCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAaaca                                                       >  1:693348/1‑150 (MQ=255)
                                                                                             tttCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAaaca                                                       <  2:827352/150‑1 (MQ=255)
                                                                                             tttCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAaaca                                                       <  2:693348/150‑1 (MQ=255)
                                                                                             tttCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAaaca                                                       <  1:573937/150‑1 (MQ=255)
                                                                                               tCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAaacaac                                                     <  2:123385/150‑1 (MQ=255)
                                                                                                   atatcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAg                                                 <  2:1077201/147‑1 (MQ=255)
                                                                                                   aTCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAg                                                 <  2:664759/150‑1 (MQ=255)
                                                                                                   aTCTCATAGAAGAGAAAGTACCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTTGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAg                                                 <  2:1083067/150‑1 (MQ=255)
                                                                                                    tctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGc                                                <  1:1302131/150‑1 (MQ=255)
                                                                                                    tctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGc                                                <  2:968549/150‑1 (MQ=255)
                                                                                                    tctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGc                                                <  1:1345341/150‑1 (MQ=255)
                                                                                                    tctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGc                                                <  1:1352100/150‑1 (MQ=255)
                                                                                                    tctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGAAGTTGAAAACAACTAAGc                                                <  2:966446/150‑1 (MQ=255)
                                                                                                     ctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGcc                                               >  2:1054736/1‑150 (MQ=255)
                                                                                                     ctcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGcc                                               >  1:247861/1‑150 (MQ=255)
                                                                                                      tcATTGCAGAGAAAGTCCATAAAGATTGTAGAGGCGTCAGCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGTCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGccc                                              <  1:1130708/150‑1 (MQ=255)
                                                                                                      tcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGccc                                              >  2:54896/1‑150 (MQ=255)
                                                                                                      tcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGccc                                              >  1:1375010/1‑150 (MQ=255)
                                                                                                      tcATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGccc                                              <  1:1000672/150‑1 (MQ=255)
                                                                                                       cATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCt                                             >  2:663930/1‑150 (MQ=255)
                                                                                                         tAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGa                                           >  1:1544593/1‑150 (MQ=255)
                                                                                                            cAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGa                                        >  1:996126/1‑150 (MQ=255)
                                                                                                             agagAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGaa                                       >  1:1448105/1‑150 (MQ=255)
                                                                                                             agagAAAGTCCCTAAAGATTGTAGAGGAGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGaa                                       <  1:229957/150‑1 (MQ=255)
                                                                                                              gagaAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAAt                                      >  2:832955/1‑150 (MQ=255)
                                                                                                                   aGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAatt                                 >  1:174873/1‑150 (MQ=255)
                                                                                                                   aGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAatt                                 <  2:1349486/150‑1 (MQ=255)
                                                                                                                       ccTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattt                             >  1:504229/1‑150 (MQ=255)
                                                                                                                        taatcgATTGTAGAGGCATGTTCTGATGACGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattatttt                            >  2:202653/6‑150 (MQ=38)
                                                                                                                        cTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGTATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattatttt                            <  1:686213/150‑1 (MQ=255)
                                                                                                                         tAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattatttta                           >  1:846153/1‑150 (MQ=255)
                                                                                                                          aaaGATTGTCGAGGCGGCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAACTACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattttat                          <  2:228901/150‑1 (MQ=255)
                                                                                                                            aGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattttatta                        >  1:819590/1‑150 (MQ=255)
                                                                                                                               ttGTAGAGGCATGTTCTGATGACGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattttattattt                     <  1:224318/150‑1 (MQ=255)
                                                                                                                                    gagGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATc                >  1:1610273/1‑150 (MQ=255)
                                                                                                                                    gagGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATc                >  1:1494365/1‑150 (MQ=255)
                                                                                                                                      ggCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattttattatttt                    >  2:15476/1‑144 (MQ=255)
                                                                                                                                      ggCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAattattttattatttt                    <  1:15476/144‑1 (MQ=255)
                                                                                                                                      ggCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCAc              <  1:461371/150‑1 (MQ=255)
                                                                                                                                        cGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCAATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGt            >  1:460862/1‑150 (MQ=255)
                                                                                                                                           catcaGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGct         >  1:1088610/1‑150 (MQ=255)
                                                                                                                                              caGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGcttcc      <  2:247861/150‑1 (MQ=255)
                                                                                                                                               aGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGcttcct     <  1:103627/150‑1 (MQ=255)
                                                                                                                                               aGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGcttcct     <  1:754594/150‑1 (MQ=255)
                                                                                                                                               aGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGcttcct     <  1:550396/150‑1 (MQ=255)
                                                                                                                                                 aaCATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGcttccttc   <  2:1049302/150‑1 (MQ=255)
                                                                                                                                                  aCATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCa  <  1:1494320/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTAACATTTTATCTTATGACATACACCTCATGTTCGATGGGAAAATAATTATATTTTGCACAAATCATTCAGCCTCAAAGCATAATATGAAATTTCTTATCTCATAGCAGAGAAAGTCCCTAAAGATTGTAGAGGCGTCATCAGAACATGCCTCTACAATCGATTAAGAAAACCTATGAAATACTGGCAGGATTAAAAAGACATACCACCGTATTCATGGAGTTTTAATGTAGTTGAAAACAACTAAGCCCTGAGGAATAAATTATTTTATTATTTTTATCACGTGCTTCCTTCA  >  NC_000913/3632627‑3632922

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