Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,775,272 C→T 7.4% A196V (GCC→GTC)  aroD → 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,775,2720CT 7.4% 209.0 / 3.8 68A196V (GCC→GTC) aroD3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base T (2/3);  total (34/34)
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.00e-01

AACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTA  >  NC_000913/1775126‑1775420
                                                                                                                                                  |                                                                                                                                                    
aaCGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGcc                                                                                                                                                     >  7:201072/1‑148 (MQ=255)
  cGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGAt                                                                                                                                                  <  4:186938/149‑1 (MQ=255)
   gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                                                                                                                 <  7:2117/149‑1 (MQ=255)
   gCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATc                                                                                                                                                 <  3:175117/149‑1 (MQ=255)
     cGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGt                                                                                                                                               >  3:334358/1‑149 (MQ=255)
      ggAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTc                                                                                                                                              <  8:201072/149‑1 (MQ=255)
        aaGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTcc                                                                                                                                             >  5:72872/1‑148 (MQ=255)
             gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                                                                                                                       >  4:316283/1‑149 (MQ=255)
             gaagaaATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTAt                                                                                                                                       >  6:248465/1‑149 (MQ=255)
                      aTTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGTCGATCGTCCAATTATCACGATGt                                                                                                                               <  2:421958/148‑1 (MQ=255)
                      aTTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTc                                                                                                                              >  6:250321/1‑149 (MQ=255)
                          cccGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATg                                                                                                                          >  6:177552/1‑149 (MQ=255)
                                  cgcAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAAc                                                                                                                  >  2:300514/1‑149 (MQ=255)
                                     aaaaTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGCCCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTgg                                                                                                               >  6:63010/1‑149 (MQ=255)
                                     aaaaTGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTgg                                                                                                               >  8:41099/1‑149 (MQ=255)
                                      aaaTGCAATCCGTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGGGCTGAGGTTGCTTGCCGCGACCCGGGAGATGCAGGAGCAGTATGCCGTTCGTCCAATTATCACGATGTCGATGGCAAAAACTGGc                                                                                                              <  6:112998/149‑1 (MQ=255)
                                           cAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAt                                                                                                         <  5:308067/149‑1 (MQ=255)
                                             aTCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAAttt                                                                                                       >  6:77383/1‑149 (MQ=255)
                                              tCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTtc                                                                                                      <  8:203883/149‑1 (MQ=255)
                                              tCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTtc                                                                                                      <  4:334358/149‑1 (MQ=255)
                                               ccTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTtct                                                                                                     >  6:343673/1‑149 (MQ=255)
                                                        ccGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggct                                                                                            >  8:58697/1‑149 (MQ=255)
                                                         cGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTctggctg                                                                                           >  6:85224/1‑149 (MQ=255)
                                                            tatTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTg                                                                                        <  8:359609/149‑1 (MQ=255)
                                                               tCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAg                                                                                     >  1:398875/1‑149 (MQ=255)
                                                                 cTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGTCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGaa                                                                                      >  7:674/1‑146 (MQ=255)
                                                                 cTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGTCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGaa                                                                                      <  4:244844/146‑1 (MQ=255)
                                                                 cTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGTCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGaa                                                                                      <  8:674/146‑1 (MQ=255)
                                                                 cTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGTCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGaa                                                                                      >  3:244844/1‑146 (MQ=255)
                                                                       ttGCGCTGCTGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGc                                                                             <  4:236813/149‑1 (MQ=255)
                                                                             tGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTcggcg                                                                       >  8:136255/1‑149 (MQ=255)
                                                                               aTGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTcggcg                                                                       <  7:58697/147‑1 (MQ=255)
                                                                                   cGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACt                                                                 <  3:316283/149‑1 (MQ=255)
                                                                                   cGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTCGCGTAATTTCTCGTCTGGCTGGTGAAGTATATGGCTCGGCGGCAACt                                                                 >  1:414500/1‑149 (MQ=255)
                                                                                     cAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACttt                                                               <  5:250321/149‑1 (MQ=255)
                                                                                        aGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTgg                                                            <  4:52012/149‑1 (MQ=255)
                                                                                          tACCAGCGATGTGCTGACGTTGATTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTACAATTATCACGATGTCGATGGCAAAAACTGGc                                                                                                              <  8:104216/97‑1 (MQ=255)
                                                                                            ccAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGc                                                         <  2:398875/148‑1 (MQ=255)
                                                                                             cAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCgg                                                       >  8:50507/1‑149 (MQ=255)
                                                                                               gCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCTAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGt                                                      >  3:255937/1‑148 (MQ=255)
                                                                                                     tgCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTaaaaaaa                                               <  7:50507/149‑1 (MQ=255)
                                                                                                       cTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCATTTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGc                                             <  5:372146/149‑1 (MQ=255)
                                                                                                        tGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTg                                                          <  6:15594/135‑1 (MQ=255)
                                                                                                          aCGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGt                                           <  5:63010/148‑1 (MQ=255)
                                                                                                                 ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                                                                                                           <  3:29497/61‑1 (MQ=255)
                                                                                                                 ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGAt                                                                                                                           >  4:29497/1‑61 (MQ=255)
                                                                                                                 ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTgcgc                                     >  8:263509/1‑147 (MQ=255)
                                                                                                                 ttgcCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAg                                  >  4:166548/1‑149 (MQ=255)
                                                                                                                     cgcgACCCTGGAGATGCAGGAGCAGTATGCCGATCGGCCAATTATCACGATGTCGATGTCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGc                               <  5:343673/149‑1 (MQ=255)
                                                                                                                           ccTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtct                         <  5:77383/149‑1 (MQ=255)
                                                                                                                           ccTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtct                         <  1:300514/149‑1 (MQ=255)
                                                                                                                           ccTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtctc                        >  7:392243/1‑149 (MQ=255)
                                                                                                                            cTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtctc                        >  7:38279/1‑149 (MQ=255)
                                                                                                                            cTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAAtctc                        >  3:323679/1‑149 (MQ=255)
                                                                                                                             tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCg                       >  2:224370/1‑149 (MQ=255)
                                                                                                                             tGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCg                       >  2:234514/1‑149 (MQ=255)
                                                                                                                               gagaTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCg                       <  5:85224/147‑1 (MQ=255)
                                                                                                                                agaTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGt                     >  4:260164/1‑148 (MQ=255)
                                                                                                                                  aTGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTaaa                  >  7:46969/1‑149 (MQ=255)
                                                                                                                                   tGCAGGATCAGTATTACGATCCTCCAATTATAACGATGTCGATGGCAAAAACTGGAGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAATAAAATCTTCTGAGACCGGGCAAATATCGGTAAAt                 >  4:317158/1‑149 (MQ=255)
                                                                                                                                        gAGCAGTATGCCGAGCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTAGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGAttt            <  5:178337/149‑1 (MQ=255)
                                                                                                                                         aGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTg           <  6:72872/149‑1 (MQ=255)
                                                                                                                                         aGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTg           <  5:248465/149‑1 (MQ=255)
                                                                                                                                          gCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTgc          <  1:234514/149‑1 (MQ=255)
                                                                                                                                           cAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTgcg         <  3:116997/149‑1 (MQ=255)
                                                                                                                                             gTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTcgg                                                                         >  2:32139/1‑83 (MQ=255)
                                                                                                                                             gTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTcgg                                                                         <  1:32139/83‑1 (MQ=255)
                                                                                                                                             gTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCa       >  8:260745/1‑149 (MQ=255)
                                                                                                                                              tATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCAc      >  1:150309/1‑149 (MQ=255)
                                                                                                                                                  ccGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTa  >  1:368606/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AACGCCGGAAGCCGAAGAAATCATTGCCCGTCTGCGCAAAATGCAATCCTTCGACGCCGATATTCCTAAGATTGCGCTGATGCCGCAAAGTACCAGCGATGTGCTGACGTTGCTTGCCGCGACCCTGGAGATGCAGGAGCAGTATGCCGATCGTCCAATTATCACGATGTCGATGGCAAAAACTGGCGTAATTTCTCGTCTGGCTGGTGAAGTATTTGGCTCGGCGGCAACTTTTGGTGCGGTAAAAAAAGCGTCTGCGCCAGGGCAAATCTCGGTAAATGATTTGCGCACGGTA  >  NC_000913/1775126‑1775420

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: