Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,418,013 G→T 6.0% V158F (GTC→TTC)  purC → phosphoribosylaminoimidazole‑succinocarboxamide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,418,0130GT 6.0% 213.4 / 3.4 67V158F (GTC→TTC) purCphosphoribosylaminoimidazole‑succinocarboxamide synthase
Reads supporting (aligned to +/- strand):  ref base G (31/32);  new base T (2/2);  total (33/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 = 8.58e-01

GGCATCAAGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCC  >  NC_002947/1417869‑1418157
                                                                                                                                                |                                                                                                                                                
ggCATCAAGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCaa                                                                                                                                              >  6:5064/1‑149 (MQ=255)
  cATCAAGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAAc                                                                                                                                             >  2:65733/1‑148 (MQ=255)
  cATCAAGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAAc                                                                                                                                             <  1:65733/148‑1 (MQ=255)
      aaGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGACCTTCAACGacag                                                                                                                                        <  5:326017/149‑4 (MQ=255)
         cTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGACCTTCAACGaca                                                                                                                                         <  2:275081/145‑3 (MQ=255)
         cTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGACCTTCAACGaca                                                                                                                                         >  1:275081/1‑143 (MQ=255)
           ggAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCt                                                                                                                                     >  2:14350/1‑147 (MQ=255)
            gAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGa                                                                                                                                   <  5:19364/148‑1 (MQ=255)
              gCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCa                                                                                                                                <  6:240136/149‑1 (MQ=255)
               ccGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCaa                                                                                                                               <  5:5064/149‑1 (MQ=255)
                 gTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGc                                                                                                                             >  6:152255/1‑149 (MQ=255)
                  tCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCt                                                                                                                            >  6:351573/1‑149 (MQ=255)
                  tCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCt                                                                                                                            >  8:275601/1‑149 (MQ=255)
                    cACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGt                                                                                                                          <  5:276619/149‑1 (MQ=255)
                          cGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCGACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGAt                                                                                                                     <  7:322362/148‑1 (MQ=255)
                              cTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGc                                                                                                                <  7:275601/149‑1 (MQ=255)
                                 ttCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCg                                                                                                              <  5:351573/148‑1 (MQ=255)
                                  tCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGc                                                                                                            >  1:108408/1‑149 (MQ=255)
                                   ccTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGATGTGCTGAGCAAGCTGTTCGATGACGCCGGcc                                                                                                           >  2:171321/1‑149 (MQ=255)
                                    cTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCct                                                                                                          >  8:6104/1‑149 (MQ=255)
                                         gaaCGACGAGAAGGGCGACCCCTTCATCAACGAATCCCTCGTTGTGGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTTGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCctgctgc                                                                                                     <  2:64172/149‑1 (MQ=255)
                                           acgacgAGCAGGGCGACCCCTTCATCCACGAATCCCACGTGGTCGCGTTCGGCTGGGTCACCGCGTAGCAACTGTTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCctgctgctg                                                                                                   <  6:26915/149‑1 (MQ=255)
                                           acgacgAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCctgctgctg                                                                                                   <  1:171321/149‑1 (MQ=255)
                                            cgacgaGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAAc                                                                                                                                             >  8:121364/1‑106 (MQ=255)
                                             gacgaGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAg                                                                                                                              >  1:239274/1‑120 (MQ=255)
                                             gacgaGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAg                                                                                                                              <  2:239274/120‑1 (MQ=255)
                                              acgaGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCTCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTc                                                                                                <  5:152255/149‑1 (MQ=255)
                                                gagaAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGa                                                                                              >  8:12202/1‑149 (MQ=255)
                                                   aaGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACtt                                                                                           <  6:133583/149‑1 (MQ=255)
                                                           ccccTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGAACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTg                                                                                    >  6:55582/1‑148 (MQ=255)
                                                            cccTTCATCAACGAATCCCACGTTGTCGCGTTAGGCTGGGGCACAGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCTTTGAAGGTCAAAGAAGTGCTGAGAAAGCTGTTCGATGAAGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGa                                                                                  >  8:20075/1‑149 (MQ=255)
                                                              cTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGt                                                                                >  2:153204/1‑149 (MQ=255)
                                                              cTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGt                                                                                >  5:318272/1‑149 (MQ=255)
                                                               ttCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCg                                                                                                                       <  6:131881/109‑1 (MQ=255)
                                                               ttCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCg                                                                                                                       >  5:131881/1‑109 (MQ=255)
                                                                     aaCGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCg                                                                          >  3:354582/1‑148 (MQ=255)
                                                                       cGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTAt                                                                       >  1:299759/1‑149 (MQ=255)
                                                                             ccACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCCGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACg                                                                 <  8:230987/149‑1 (MQ=255)
                                                                              cACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGt                                                                                                 >  4:26969/1‑116 (MQ=255)
                                                                              cACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGt                                                                                                 <  3:26969/116‑1 (MQ=255)
                                                                              cACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTGGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTAGTACACgg                                                                >  6:260381/1‑149 (MQ=255)
                                                                                 gTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCg                                                                          >  4:139513/1‑136 (MQ=255)
                                                                                 gTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCg                                                                          <  3:139513/136‑1 (MQ=255)
                                                                                  ttGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGc                                                               >  1:348881/1‑146 (MQ=255)
                                                                                     tCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATc                                                         <  7:173812/149‑1 (MQ=255)
                                                                                       gcgTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGt                                                       >  3:246497/1‑149 (MQ=255)
                                                                                          ttCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGACGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCt                                                    >  1:74585/1‑149 (MQ=255)
                                                                                            cGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTg                                                   <  4:354582/148‑1 (MQ=255)
                                                                                             ggCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTgg                                                  <  7:12202/148‑1 (MQ=255)
                                                                                              gCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGc                                                <  1:14350/149‑1 (MQ=255)
                                                                                               cTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTggg                                                 <  2:23759/147‑1 (MQ=255)
                                                                                               cTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTAGGCGTATTCCACGGCCAGATCGTGCTggg                                                 >  1:23759/1‑147 (MQ=255)
                                                                                                          ccgAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGc                                    <  2:108408/149‑1 (MQ=255)
                                                                                                                    ggtggAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCt                          <  5:260381/149‑1 (MQ=255)
                                                                                                                      tggAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTgc                        <  2:299759/149‑1 (MQ=255)
                                                                                                                             gaagaaGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCACGCTGGAGTTCGGCGTATTCCACGGcc                                                              <  1:192472/104‑1 (MQ=255)
                                                                                                                             gaagaaGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGcc                                                              >  2:192472/1‑104 (MQ=255)
                                                                                                                             gaagaaGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCtgt                 >  3:302509/1‑149 (MQ=255)
                                                                                                                              aagaagCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCtgtg                <  5:55582/149‑1 (MQ=255)
                                                                                                                               agaagGTGTCGTTGAAGGTCAACGAAGTGCTGAGGATGCTGTTCGATGACGCCGGGCTGGTGCTGGTGGACTTCAAGCTGTAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCGtgtg                <  4:246497/148‑1 (MQ=255)
                                                                                                                               agaagCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTTGAGTTCGGCGTATTCCACGGCCATATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTgg               <  7:20075/149‑1 (MQ=255)
                                                                                                                               agaagCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTgg                                                  >  2:264724/1‑114 (MQ=255)
                                                                                                                               agaagCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTgg                                                  <  1:264724/114‑1 (MQ=255)
                                                                                                                                  gcacttCGTTGACCTTCAACGACGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCAGGAGTTCCGCGTATTCCACGGCCa                                                             >  3:269927/6‑100 (MQ=255)
                                                                                                                                         gTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGaaa     >  3:105652/1‑149 (MQ=255)
                                                                                                                                         gTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGaa      <  2:74585/148‑1 (MQ=255)
                                                                                                                                            gAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAAccc  >  5:30196/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GGCATCAAGCTGGAGCCGTCCACCTTCGAGCTGTTCCTGAAGAACGACGAGAAGGGCGACCCCTTCATCAACGAATCCCACGTTGTCGCGTTCGGCTGGGGCACCGCCGAGCAACTGGTGGAAATGAAGAAGCTGTCGTTGAAGGTCAACGAAGTGCTGAGCAAGCTGTTCGATGACGCCGGCCTGCTGCTGGTCGACTTCAAGCTGGAGTTCGGCGTATTCCACGGCCAGATCGTGCTGGGCGACGAGTTCAGCCCGGACGGCTGCCGCCTGTGGGACAAGGAAACCC  >  NC_002947/1417869‑1418157

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: