Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,074,559 G→A 6.4% T44M (ACG→ATG)  yggC ← putative PanK family P‑loop kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,074,5590GA 6.4% 249.8 / 5.0 78T44M (ACG→ATG) yggCputative PanK family P‑loop kinase
Reads supporting (aligned to +/- strand):  ref base G (36/37);  new base A (3/2);  total (39/39)
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.95e-01

ATCCAGCCAGCTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCTTTAAG  >  NC_000913/3074414‑3074707
                                                                                                                                                 |                                                                                                                                                    
aTCCAGCCAGCTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTAc                                                                                                                                                   >  3:194447/1‑149 (MQ=255)
  ccagccagCTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTAcgc                                                                                                                                                 >  7:45045/1‑149 (MQ=255)
   cagccagcTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTAcgcc                                                                                                                                                <  7:220221/149‑1 (MQ=255)
   cagccagcTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTAcgcc                                                                                                                                                >  7:348291/1‑149 (MQ=255)
    agccagcTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTAcgccg                                                                                                                                               >  4:72273/1‑149 (MQ=255)
        agcTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGc                                                                                                                                           <  5:15182/149‑1 (MQ=255)
           tATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGa                                                                                                                                        >  3:312711/1‑149 (MQ=255)
              tGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGAtt                                                                                                                                      >  5:372736/1‑148 (MQ=255)
                tAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTAc                                                                                                                                   >  4:177914/1‑149 (MQ=255)
                  aTGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCg                                                                                                                                 >  4:267258/1‑149 (MQ=255)
                   tGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGt                                                                                                                                <  4:312711/149‑1 (MQ=255)
                   tGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGt                                                                                                                                <  3:177914/149‑1 (MQ=255)
                        aaaaCCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCt                                                                                                                                             <  4:100755/131‑1 (MQ=255)
                        aaaaCCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCt                                                                                                                                             >  3:100755/1‑131 (MQ=255)
                              atccatcGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGGTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGgccg                                                                                                                     >  6:89199/1‑149 (MQ=255)
                               tccatcGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGgccgc                                                                                                                    <  8:45045/149‑1 (MQ=255)
                                 catcGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCACGTCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGTATTTACCGTCTGTAAGGCCGCCa                                                                                                                  <  4:314927/149‑1 (MQ=255)
                                     gggAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACa                                                                                                               <  8:116253/148‑1 (MQ=255)
                                      ggAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACAtg                                                                                                             <  6:335498/149‑1 (MQ=255)
                                       gAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACAtgt                                                                                                            <  4:194447/149‑1 (MQ=255)
                                          cttCTGAATAGCGGGCAGTTCCGGGTCTTGTTTCGCGAGATATTCCCATAAGGTTGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCTGATTTAACGTCTGTAATGCCGCCAACATGTgc                                                                                                          <  8:348291/146‑1 (MQ=255)
                                            tCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTgcag                                                                                                        >  3:362835/1‑148 (MQ=255)
                                            gcTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTgcag                                                                                                        >  7:289294/2‑148 (MQ=255)
                                              tGAATAGCGGGCAGTTCCGGGTCTTGTTGCGAGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCaa                                                                                                     >  8:245205/1‑149 (MQ=255)
                                                     cGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTa                                                                                              <  3:75225/149‑1 (MQ=255)
                                                     cGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTa                                                                                              >  4:75225/1‑149 (MQ=255)
                                                     cGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTa                                                                                              <  6:372736/149‑1 (MQ=255)
                                                               cGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACAtttt                                                                                    >  3:280923/1‑149 (MQ=255)
                                                                gggTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTc                                                                                   <  1:37155/149‑1 (MQ=255)
                                                                 ggTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCa                                                                                  <  1:75474/149‑1 (MQ=255)
                                                                       gttgCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATctct                                                                            >  6:181965/1‑149 (MQ=255)
                                                                       gttgCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGACTGTAAGGCCGCCACCATGTGCAGCAACGGTTTATGTACATTTTCAATctct                                                                            >  5:226972/1‑149 (MQ=255)
                                                                          gcgcGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTc                                                                          <  7:245205/148‑1 (MQ=255)
                                                                          gcgcGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTc                                                                          <  3:267258/148‑1 (MQ=255)
                                                                          gcgcGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTc                                                                          <  7:153419/148‑1 (MQ=255)
                                                                           cgcgAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGt                                                                        <  4:280923/149‑1 (MQ=255)
                                                                              gagaTATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCt                                                                     >  2:185863/1‑149 (MQ=255)
                                                                                gaTATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAAGCACCGTACGCCGCTGCTGATTTACCTTCTGTAAGGCCTCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTg                                                                    <  7:308609/148‑1 (MQ=255)
                                                                                 atatTCCCAGAAGGTGGTTAGGGGTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGt                                                                  <  5:89199/149‑1 (MQ=255)
                                                                                   atTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTTTGTACATTTTCAATCTCTTc                                                                          <  1:185872/139‑1 (MQ=255)
                                                                                   atTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTc                                                                          >  2:185872/1‑139 (MQ=255)
                                                                                   atTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTATTCGTTCTGGTAc                                                                >  7:219337/1‑149 (MQ=255)
                                                                                    ttCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACt                                                               <  5:342624/149‑1 (MQ=255)
                                                                                    ttCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACt                                                               >  5:385541/1‑149 (MQ=255)
                                                                                       ccAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACCTGTGCAGCAACGGTTTATGTACATTTTCACTCTCTTCGTACTGGTACtgtg                                                            >  1:229342/1‑149 (MQ=255)
                                                                                        cAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGc                                                           <  8:86672/149‑1 (MQ=255)
                                                                                            aGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGt                                                       <  3:72273/149‑1 (MQ=255)
                                                                                             ggaggtTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCATACGCCGCTGCGGATTTACCGTCTGTATGGCCGCCAACATGTGCAGCAAAGGTTTATGTCAAATTTCAATCTCTTCGTTCTGGAACTGTGACTCTa                                                      >  2:291639/4‑149 (MQ=255)
                                                                                                       gTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTg                                                                                                                             >  1:312877/1‑68 (MQ=255)
                                                                                                       gTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTg                                                                                                                             <  2:312877/68‑1 (MQ=255)
                                                                                                          gATTTACCCGTCCCCGGCGGCGCGCAAAGAAAACCCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCAt                                          >  5:258729/1‑148 (MQ=255)
                                                                                                               aCCCGTCCCCGGAGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTACTGGTACTGTGCCTGTATCTTCAGCCCATTCACCg                                    >  6:64374/1‑149 (MQ=255)
                                                                                                                 ccGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGt                                   >  1:381624/1‑148 (MQ=255)
                                                                                                                     ccccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCtt                                                  >  8:379808/1‑129 (MQ=255)
                                                                                                                     ccccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCtt                                                  <  7:379808/129‑1 (MQ=255)
                                                                                                                     ccccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAAtt                              >  6:366418/1‑149 (MQ=255)
                                                                                                                      cccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCa                                           <  5:387542/135‑1 (MQ=255)
                                                                                                                      cccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCa                                           >  6:387542/1‑135 (MQ=255)
                                                                                                                      cccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAAtt                              >  7:341298/1‑148 (MQ=255)
                                                                                                                      cccGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAAtt                              >  2:175793/1‑148 (MQ=255)
                                                                                                                        cggcggAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCaa                           <  4:362835/149‑1 (MQ=255)
                                                                                                                          gcggAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAAtt                         <  1:353055/149‑1 (MQ=255)
                                                                                                                              aGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCa                     <  1:185863/149‑1 (MQ=255)
                                                                                                                                  cAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGc                 >  4:16101/1‑149 (MQ=255)
                                                                                                                                  cAAAGAAAAACCACCATACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTc                                                                                   <  8:279775/83‑1 (MQ=255)
                                                                                                                                  cAAAGAAAAACCACCATACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTc                                                                                   >  7:279775/1‑83 (MQ=255)
                                                                                                                                  cAAAGAAAAACCACCATACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTc                                                                                   >  5:54797/1‑83 (MQ=255)
                                                                                                                                  cAAAGAAAAACCACCATACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTc                                                                                   <  6:54797/83‑1 (MQ=255)
                                                                                                                                   aaaGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCg                <  1:175793/149‑1 (MQ=255)
                                                                                                                                   aaaGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCg                >  4:141446/1‑149 (MQ=255)
                                                                                                                                      gAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGtt              >  6:67893/1‑148 (MQ=255)
                                                                                                                                       aaaaaCCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTaa            <  8:289294/149‑1 (MQ=255)
                                                                                                                                         aaaCCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAAtt          <  8:341298/149‑1 (MQ=255)
                                                                                                                                             caccGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCtt      <  6:385541/149‑1 (MQ=255)
                                                                                                                                             caccGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCtt      >  3:67486/1‑149 (MQ=255)
                                                                                                                                             caccGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGCGCCTGGATCTTAAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCtt      >  6:5493/1‑149 (MQ=255)
                                                                                                                                              accGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCttt     <  5:181965/149‑1 (MQ=255)
                                                                                                                                               ccGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGAAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCTTTa    <  6:258729/149‑1 (MQ=255)
                                                                                                                                                 gTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCTTTAAg  <  4:67486/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
ATCCAGCCAGCTATTGTAATGGTGAAAACCATCCATCGGGAGCGTCTGAATAGCGGGCAGTTCCGGGTCTTGTTGCGCGAGATATTCCCAGAAGGTGGTTAGGGTTGATTTACCCGTCCCCGGCGGAGCGCAAAGAAAAACCACCGTACGCCGCTGCGGATTTACCGTCTGTAAGGCCGCCAACATGTGCAGCAACGGTTTATGTACATTTTCAATCTCTTCGTTCTGGTACTGTGCCTGTATCTTCAGCCCATTCACCGTTAATTCAATTTTCACGGCGTTAATTCCTTTAAG  >  NC_000913/3074414‑3074707

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