Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,478,806 T→A 5.2% N329I (AAC→ATC) ‡ proA ← gamma‑glutamyl phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,478,8060TA 5.2% 1126.8 / 32.9 328N329I (AAC→ATC) ‡proAgamma‑glutamyl phosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (167/143);  new base A (6/11);  total (173/155)
Fisher's exact test for biased strand distribution p-value = 1.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 332 total aligned reads displayed.

GAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATT  >  NC_002947/5478658‑5478955
                                                                                                                                                    |                                                                                                                                                     
gAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGtt                                                                                                                                                      >  2:1329302/1‑150 (MQ=255)
  aGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCa                                                                                                                                                    <  2:2082482/150‑1 (MQ=255)
  aGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCa                                                                                                                                                    >  2:572571/1‑150 (MQ=255)
   gCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAg                                                                                                                                                   <  2:3284332/150‑1 (MQ=255)
   gCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAg                                                                                                                                                   <  2:6139550/150‑1 (MQ=255)
   gCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAg                                                                                                                                                   <  1:3007915/150‑1 (MQ=255)
    ccATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAgg                                                                                                                                                  <  1:3071040/150‑1 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                >  1:6412754/1‑150 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                >  2:3786531/1‑150 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                >  2:6234903/1‑150 (MQ=255)
      aTCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGCGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGcc                                                                                                                                                >  2:102101/1‑150 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGAACCAGGCCATc                                                                                                                                             <  2:4599755/150‑1 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGAACCAGGCCATc                                                                                                                                             <  2:3607255/150‑1 (MQ=255)
         ggCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGAACCAGGCCATc                                                                                                                                             <  2:2428934/150‑1 (MQ=255)
          gCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCATGGTTCAGGCCGTCg                                                                                                                                            >  2:3663020/1‑150 (MQ=255)
           cAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGa                                                                                                                                           <  2:5061391/150‑1 (MQ=255)
           cAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGa                                                                                                                                           <  1:892260/150‑1 (MQ=255)
           cAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGACGGCCTGAACCAGGCCATCGa                                                                                                                                           <  1:4879083/150‑1 (MQ=255)
            aaaGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGAc                                                                                                                                          >  1:2036734/1‑150 (MQ=255)
              aGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCa                                                                                                                                        <  2:3147133/150‑1 (MQ=255)
                 aGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGGGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCAcgc                                                                                                                                     <  2:883942/150‑1 (MQ=255)
                      ggggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATc                                                                                                                                >  2:4522057/1‑150 (MQ=255)
                        ggTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGa                                                                                                                              <  2:2989479/150‑1 (MQ=255)
                          tgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                            >  1:2056844/1‑150 (MQ=255)
                          tgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                            <  2:5639237/150‑1 (MQ=255)
                          tgtTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACa                                                                                                                            <  1:541092/150‑1 (MQ=255)
                              gAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGAt                                                                                                                        >  2:1495825/1‑150 (MQ=255)
                              gAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGAt                                                                                                                        <  2:5002931/150‑1 (MQ=255)
                                      ccGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGa                                                                                                                <  2:1749507/150‑1 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                               >  1:6775717/1‑150 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                               <  1:6564835/150‑1 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAg                                                                                                               >  1:1554412/1‑150 (MQ=255)
                                       cGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGACGATCGACACGGTGTGGCGCGAGCCAGAGTCGTTGATGTGCTCGATGGCCTGGTTCAGGCCCTCGACCACGCGCCACGACAGGATCGCAACGAg                                                                                                               <  2:2924436/150‑1 (MQ=255)
                                              gAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTcc                                                                                                        <  1:4074548/150‑1 (MQ=255)
                                              gAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTcc                                                                                                        >  2:5396621/1‑150 (MQ=255)
                                              gAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTcc                                                                                                        >  2:4224798/1‑150 (MQ=255)
                                               aaTCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCg                                                                                                       >  2:1451581/1‑150 (MQ=255)
                                                 tCGACCTCAGCCATGAACTGACGGGCTTCTCCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCgtg                                                                                                     <  1:3611445/150‑1 (MQ=255)
                                                  cGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCgtgt                                                                                                    <  1:6132422/150‑1 (MQ=255)
                                                    aCCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                  <  1:2429793/150‑1 (MQ=255)
                                                    aCCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGc                                                                                                  <  1:1839938/150‑1 (MQ=255)
                                                      cTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCa                                                                                                >  1:730983/1‑150 (MQ=255)
                                                           ccATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCgg                                                                                           <  1:117739/150‑1 (MQ=255)
                                                             aTGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCt                                                                                         >  1:58486/1‑150 (MQ=255)
                                                             aTGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCt                                                                                         >  2:1958731/1‑150 (MQ=255)
                                                               gAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTc                                                                                       <  1:1769522/150‑1 (MQ=255)
                                                                  cTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGt                                                                                    <  1:3445405/150‑1 (MQ=255)
                                                                   tGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                   <  1:3219818/150‑1 (MQ=255)
                                                                   tGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTg                                                                                   >  1:3517103/1‑150 (MQ=255)
                                                                    gACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTgg                                                                                  <  2:2568928/150‑1 (MQ=255)
                                                                     aCGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGc                                                                                 <  2:5242517/150‑1 (MQ=255)
                                                                     aCGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGc                                                                                 >  2:6274147/1‑150 (MQ=255)
                                                                       gggCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGa                                                                                                                                           >  2:4693963/1‑90 (MQ=255)
                                                                       gggCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGa                                                                                                                                           <  1:4693963/90‑1 (MQ=255)
                                                                          cTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGtt                                                                            <  1:1457236/150‑1 (MQ=255)
                                                                           ttCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGttt                                                                           >  2:3739401/1‑150 (MQ=255)
                                                                           ttCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGttt                                                                           >  2:2612327/1‑150 (MQ=255)
                                                                             cACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTgg                                                                         >  1:3311705/1‑150 (MQ=255)
                                                                               cccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCgg                                                                                     <  1:3762725/136‑1 (MQ=255)
                                                                                ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcgc                                                                      <  1:4437898/150‑1 (MQ=255)
                                                                                ccTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGGTTGgcgc                                                                      >  2:6818375/1‑150 (MQ=255)
                                                                                  tGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtg                                                                    >  1:4058804/1‑150 (MQ=255)
                                                                                   ggTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtga                                                                   <  1:1964605/150‑1 (MQ=255)
                                                                                      gtTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCtgatga                                                                <  2:3545918/150‑1 (MQ=255)
                                                                                       ttCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTgatgat                                                               <  1:2837881/150‑1 (MQ=255)
                                                                                        tCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATa                                                              <  2:4071439/150‑1 (MQ=255)
                                                                                        tCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATa                                                              <  2:3334009/150‑1 (MQ=255)
                                                                                          gCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGc                                                            <  1:1509972/150‑1 (MQ=255)
                                                                                            tgatgaTCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCt                                                          >  2:5128039/1‑150 (MQ=255)
                                                                                            tgatgaTCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCt                                                          >  1:2322390/1‑150 (MQ=255)
                                                                                                               gtgCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACcgcg                                       >  1:520500/1‑150 (MQ=255)
                                                                                                                tgCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCgcgc                                      >  2:1569471/1‑150 (MQ=255)
                                                                                                                 gCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCa                                     <  1:1002103/150‑1 (MQ=255)
                                                                                                                        aTAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCAc                              >  1:2651338/1‑150 (MQ=255)
                                                                                                                          aGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTAc                                                    <  1:636897/126‑1 (MQ=255)
                                                                                                                          aGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGc                            >  2:6517152/1‑150 (MQ=255)
                                                                                                                            tGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGccc                          <  1:319942/150‑1 (MQ=255)
                                                                                                                            tGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGccc                          <  2:4733751/150‑1 (MQ=255)
                                                                                                                             ggTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCt                         >  1:2544505/1‑150 (MQ=255)
                                                                                                                               ttGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtc                       <  2:642128/150‑1 (MQ=255)
                                                                                                                               ttGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtc                       <  1:572598/150‑1 (MQ=255)
                                                                                                                                tGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtct                      <  2:3437674/150‑1 (MQ=255)
                                                                                                                                 gATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTtctc                     >  1:4272019/1‑150 (MQ=255)
                                                                                                                                   tgtgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                   >  2:5912104/1‑150 (MQ=255)
                                                                                                                                   cgtgGTCGACGGCCTGAACCAGGCCATCGAGCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCa                   >  1:6675453/2‑150 (MQ=255)
                                                                                                                                     tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGgcg                                                                       <  1:3899481/96‑1 (MQ=255)
                                                                                                                                     tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACg                 >  1:5139317/1‑150 (MQ=255)
                                                                                                                                     tgCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACg                 >  1:5479885/1‑150 (MQ=255)
                                                                                                                                      gCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGa                >  1:1924272/1‑150 (MQ=255)
                                                                                                                                       cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGaa               >  1:5656156/1‑150 (MQ=255)
                                                                                                                                       cTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGaa               <  1:2305670/150‑1 (MQ=255)
                                                                                                                                         cGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAgc             <  1:1709604/150‑1 (MQ=255)
                                                                                                                                               ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGATGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg       >  1:2668828/1‑150 (MQ=255)
                                                                                                                                               ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg       >  1:4979583/1‑150 (MQ=255)
                                                                                                                                               ccTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCggg       <  1:6067484/150‑1 (MQ=255)
                                                                                                                                                cTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGc      >  1:4611879/1‑150 (MQ=255)
                                                                                                                                                cTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGc      >  2:6644566/1‑150 (MQ=255)
                                                                                                                                                 tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGcc     >  1:1619096/1‑150 (MQ=255)
                                                                                                                                                 tGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGcc     >  2:6785018/1‑150 (MQ=255)
                                                                                                                                                    ttCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCAtt  >  1:6206020/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAAGCCATCGGCAAAGCAGGTCGGGGTGTTGAGCATGACCGACGCCGAATCGACCTCAGCCATGAACTGACGGGCTTCACCCTGGTGTTCGCTGATGATCGAGTCGGTGTGGTGCGAGCCATAGTGGTTGATGTGCTCGATGGCCTGGTTCAGGCCGTCGACCACGCGAATCGACAGGATCGCATCGAGGTATTCCGTGTGCCAGTCGGCTTCGGTGGCCGGTTTGGCGCTGATGATAGCCTGGGTACGCTCGCAACCGCGCAGCTCCACGCCCTTCTCCACGAAGCGGCGGGCCATT  >  NC_002947/5478658‑5478955

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