Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 351,452 C→G 6.4% Q80E (CAG→GAG)  prpD → 2‑methylcitrate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913351,4520CG 6.4% 256.2 / 5.0 76Q80E (CAG→GAG) prpD2‑methylcitrate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (33/38);  new base G (2/3);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTCGCGC  >  NC_000913/351303‑351601
                                                                                                                                                     |                                                                                                                                                     
ctAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTc                                                                                                                                                       <  2:126682/150‑1 (MQ=255)
    aGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGtt                                                                                                                                                   <  1:952223/150‑1 (MQ=255)
     gTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTc                                                                                                                                                  <  2:1310352/150‑1 (MQ=255)
      gggCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTcc                                                                                                                                                 <  2:409915/149‑1 (MQ=255)
        gCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTGAGTTccgg                                                                                                                                               >  2:177467/1‑148 (MQ=255)
        gCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAg                                                                                                                                               <  1:502554/150‑1 (MQ=255)
         ccTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTGAGTTccggg                                                                                                                                              <  1:177467/150‑4 (MQ=255)
           tACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTc                                                                                                                                            <  1:473811/150‑1 (MQ=255)
           tACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTc                                                                                                                                            >  2:417842/1‑150 (MQ=255)
            aCGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTGAGTTccggggac                                                                                                                                           <  2:577987/150‑7 (MQ=255)
                cacCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGAccc                                                                                                                                       <  1:332429/150‑1 (MQ=255)
                  ccGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAGAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGCCCCCg                                                                                                                                     >  2:865361/1‑150 (MQ=255)
                  ccGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCg                                                                                                                                     >  2:1311645/1‑150 (MQ=255)
                    gCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTc                                                                                                                                   <  2:144087/150‑1 (MQ=255)
                     cacaTTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTcc                                                                                                                                  <  2:48410/150‑1 (MQ=255)
                        aTTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAgg                                                                                                                               >  1:279455/1‑150 (MQ=255)
                                ctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGt                                                                                                                                    >  2:365877/1‑137 (MQ=255)
                                ctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGt                                                                                                                                    <  1:365877/137‑1 (MQ=255)
                                ctgcTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCAttt                                                                                                                       <  2:9733/150‑1 (MQ=255)
                                    tCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACa                                                                                                                   <  1:1485210/150‑1 (MQ=255)
                                     cGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACAt                                                                                                                  >  1:801064/1‑150 (MQ=255)
                                        cacGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCgg                                                                                                               <  1:257994/150‑1 (MQ=255)
                                          cGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGgcg                                                                                                             >  2:131599/1‑150 (MQ=255)
                                                   gCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATcc                                                                                                    <  2:716336/150‑1 (MQ=255)
                                                          ggAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCt                                                                                             >  2:950862/1‑150 (MQ=255)
                                                           gAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTc                                                                                            <  1:417842/150‑1 (MQ=255)
                                                           gAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTc                                                                                            >  2:689626/1‑150 (MQ=255)
                                                            aaGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCg                                                                                           <  1:219091/150‑1 (MQ=255)
                                                             aGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGGTCCCGGCACCGGCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGa                                                                                          <  1:1435064/150‑1 (MQ=255)
                                                                tctcGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGAttt                                                                                       >  2:547494/1‑150 (MQ=255)
                                                                 ctcGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTc                                                                                      <  1:689626/150‑1 (MQ=255)
                                                                 ctcGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTc                                                                                      <  1:602416/150‑1 (MQ=255)
                                                                   cGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCaa                                                                                    <  1:1366121/150‑1 (MQ=255)
                                                                   cGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCATACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCaa                                                                                    <  1:627594/150‑1 (MQ=255)
                                                                        aCCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATa                                                                               >  1:787372/1‑150 (MQ=255)
                                                                         cccGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATAc                                                                              >  1:1194021/1‑150 (MQ=255)
                                                                            ggCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACctg                                                                           >  1:1419824/1‑150 (MQ=255)
                                                                                 gTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACctggctgg                                                                      >  1:865494/1‑150 (MQ=255)
                                                                                  tAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCtggctggc                                                                     >  1:1320169/1‑150 (MQ=255)
                                                                                    aaaaaCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCGGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTggcgg                                                                   >  1:1457641/1‑150 (MQ=255)
                                                                                     aaaaCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTggcggc                                                                  >  1:385433/1‑150 (MQ=255)
                                                                                      aaaCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGgcggcg                                                                 <  1:561666/150‑1 (MQ=255)
                                                                                      aaaCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGgcggcg                                                                 <  2:27738/150‑1 (MQ=255)
                                                                                      aaaCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATAGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGgcggcg                                                                 >  2:1280937/1‑150 (MQ=255)
                                                                                        aCTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGa                                                               <  1:1311645/150‑1 (MQ=255)
                                                                                         ctgctgGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAg                                                              <  2:245163/150‑1 (MQ=255)
                                                                                           gctgGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTg                                                            <  1:363961/150‑1 (MQ=255)
                                                                                            gtgGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGGGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTgg                                                           <  1:687056/149‑1 (MQ=255)
                                                                                            ctgGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTgg                                                           >  2:1123468/1‑150 (MQ=255)
                                                                                              ggggCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTgggg                                                         <  1:1199616/150‑1 (MQ=255)
                                                                                               gggCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGc                                                        <  1:563529/150‑1 (MQ=255)
                                                                                                 gCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCa                                                      >  1:964100/1‑150 (MQ=255)
                                                                                                   cAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCTTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATc                                                    <  1:925663/150‑1 (MQ=255)
                                                                                                         ttCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCCACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCg                                              >  1:1358911/1‑150 (MQ=255)
                                                                                                         ttCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCg                                              <  2:1350869/150‑1 (MQ=255)
                                                                                                            ccGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACa                                           <  1:1381982/150‑1 (MQ=255)
                                                                                                             cGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCCGCGCGATGATCCGCTGGCTCGATTTCAACGATAACTGGCTGGCGGCGGAGTGGGGCAATCCTTACGACaa                                          >  2:1125185/1‑150 (MQ=255)
                                                                                                              ggCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAAc                                         >  1:1358403/1‑150 (MQ=255)
                                                                                                               gCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAAcc                                        >  1:857420/1‑150 (MQ=255)
                                                                                                                 aCCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTc                                      <  1:547494/150‑1 (MQ=255)
                                                                                                                   cgtcgtACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCTGAGTGGGGCCATCCTTCCGACAACCTcgg                                    <  2:385433/150‑1 (MQ=255)
                                                                                                                     gcgtACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTcggcg                                  <  1:685145/149‑1 (MQ=255)
                                                                                                                      cgtACCCAACGGCGTGCGCGTCCCCGGAACTGAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTcggcgg                                 <  2:286829/150‑1 (MQ=255)
                                                                                                                      cgtACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTcggcgg                                 >  1:319048/1‑150 (MQ=255)
                                                                                                                         aCCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCAt                              >  2:91444/1‑150 (MQ=255)
                                                                                                                           ccAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTc                            >  1:1369281/1‑150 (MQ=255)
                                                                                                                                        cgTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGa               >  2:562464/1‑150 (MQ=255)
                                                                                                                                            gctggAACTGAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTgg           >  1:981050/4‑150 (MQ=255)
                                                                                                                                             ccGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGc          >  2:303105/1‑150 (MQ=255)
                                                                                                                                             ccGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGc          >  2:972990/1‑150 (MQ=255)
                                                                                                                                              cGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCt         <  1:1125185/150‑1 (MQ=255)
                                                                                                                                               ggAACTCGGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCtt        >  2:559241/1‑150 (MQ=255)
                                                                                                                                                gAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCttt       <  1:464219/150‑1 (MQ=255)
                                                                                                                                                  aCTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGAGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTcg     >  1:1005323/1‑150 (MQ=255)
                                                                                                                                                   cTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTcgc    <  2:633489/150‑1 (MQ=255)
                                                                                                                                                   cTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTcgc    >  1:536892/1‑150 (MQ=255)
                                                                                                                                                    tCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTcgcg   <  1:131599/150‑1 (MQ=255)
                                                                                                                                                     cAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTCgcgc  >  1:824136/1‑150 (MQ=255)
                                                                                                                                                     cAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTCgcgc  >  2:268954/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTAAAGTGGCCTACGACACCGCACATTACTGCCTGCTCGACACGCTCGGCTGCGGTCTGGAAGCTCTCGAATACCCGGCCTGTAAAAAACTGCTGGGGCCAATTGTTCCCGGCACCGTCGTACCCAACGGCGTGCGCGTCCCCGGAACTCAGTTCCAGCTCGACCCCGTCCAGGCGGCATTTAACATCGGCGCGATGATCCGCTGGCTCGATTTCAACGATACCTGGCTGGCGGCGGAGTGGGGCCATCCTTCCGACAACCTCGGCGGCATTCTGGCAACGGCGGACTGGCTTTCGCGC  >  NC_000913/351303‑351601

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