Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,528,007 T→A 5.6% F351L (TTT→TTA) ‡ purB → adenylosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,528,0070TA 5.6% 220.3 / 2.5 71F351L (TTT→TTA) ‡purBadenylosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (2/2);  total (34/37)
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.41e-01

AGGTCAACCCGATCGACTTCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAA  >  NC_002947/4527859‑4528154
                                                                                                                                                    |                                                                                                                                                   
aGGTCAACCCGATCGACTTCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCttt                                                                                                                                                     <  7:666/149‑1 (MQ=255)
       cccGATCGACTTCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCca                                                                                                                                                >  2:130311/1‑147 (MQ=255)
        ccGATCGACTTCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCTGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAg                                                                                                                                             >  2:52349/1‑149 (MQ=255)
                 ttCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGc                                                                                                                                    >  4:223469/1‑149 (MQ=255)
                    gAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCtt                                                                                                                                 >  3:96727/1‑149 (MQ=255)
                          tCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGc                                                                                                                           >  7:301213/1‑149 (MQ=255)
                          tCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGc                                                                                                                           >  7:228212/1‑149 (MQ=255)
                          tCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCAAAGCCcacgccc                                                                                                                                           >  8:1022/1‑129 (MQ=255)
                          tCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCAAAGCCcacgccc                                                                                                                                           <  7:1022/133‑5 (MQ=255)
                            cGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCa                                                                                                                         <  6:9571/149‑1 (MQ=255)
                              aaGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGc                                                                                                                       >  1:251803/1‑149 (MQ=255)
                                  cAACCTGGGTATCGCCAACGCGCTGTTCCTGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGa                                                                                                                   <  1:52349/149‑1 (MQ=255)
                                    aCCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCCTCGAGGCCAGCCTGa                                                                                                                   >  8:226622/1‑147 (MQ=255)
                                         ggTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCAt                                                                                                            <  8:90358/149‑1 (MQ=255)
                                          gTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATc                                                                                                           >  3:211772/1‑149 (MQ=255)
                                               gCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGaaa                                                                                                                 >  1:124621/1‑138 (MQ=255)
                                               gCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGaaa                                                                                                                 <  2:124621/138‑1 (MQ=255)
                                               gCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGc                                                                                                        >  2:236869/1‑147 (MQ=255)
                                               gCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCCTCGAGGCCAGCCTGAAAGGCATCGGCaa                                                                                                      <  7:226622/149‑1 (MQ=255)
                                                    cgcgCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAg                                                                                                                <  7:269165/134‑1 (MQ=255)
                                                    cgcgCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAg                                                                                                                >  8:269165/1‑134 (MQ=255)
                                                      cgcTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCTTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGaa                                                                                               <  8:301213/149‑1 (MQ=255)
                                                      cgcTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGaa                                                                                               <  5:44770/149‑1 (MQ=255)
                                                      cgcTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGaa                                                                                               <  2:251803/149‑1 (MQ=255)
                                                       gcTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAg                                                                                              >  8:44105/1‑149 (MQ=255)
                                                           ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCaa                                                                                          >  8:318971/1‑149 (MQ=255)
                                                           ttCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCaa                                                                                          <  3:223469/149‑1 (MQ=255)
                                                             ccAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCAAAGCCcacgc                                                                                                                                             <  3:165818/96‑3 (MQ=255)
                                                             ccAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCAAAGCCcacgc                                                                                                                                             >  4:165818/1‑94 (MQ=255)
                                                              cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGc                                                                                                    <  3:209028/136‑1 (MQ=255)
                                                              cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGc                                                                                                    >  4:209028/1‑136 (MQ=255)
                                                              cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACg                                                                                        >  7:208616/1‑148 (MQ=255)
                                                              cAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGa                                                                                       <  1:130311/149‑1 (MQ=255)
                                                               aGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACg                                                                                        >  2:44990/1‑147 (MQ=255)
                                                                  aCCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGc                                                                                                                                                   >  1:100497/1‑85 (MQ=255)
                                                                  aCCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGc                                                                                                                                                   <  2:100497/85‑1 (MQ=255)
                                                                             aaGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATcgccgc                                                                        <  7:318971/149‑1 (MQ=255)
                                                                               gCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCgccgcc                                                                       >  1:303347/1‑148 (MQ=255)
                                                                               gCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCgccgcc                                                                       >  1:95137/1‑148 (MQ=255)
                                                                                  gCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGAcc                                                                   >  3:49448/1‑149 (MQ=255)
                                                                                     gATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTg                                                                 >  4:284248/1‑148 (MQ=255)
                                                                                           gcgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACg                                                             <  7:44105/146‑1 (MQ=255)
                                                                                           gcgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGc                                                            <  8:208616/147‑1 (MQ=255)
                                                                                            cgcTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCctg                                                         <  4:96727/149‑1 (MQ=255)
                                                                                                 gCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCctgctg                                                      <  4:49448/147‑1 (MQ=255)
                                                                                                 gCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCctgctg                                                      <  4:968/147‑1 (MQ=255)
                                                                                                 gCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCctgctg                                                      <  1:236869/147‑1 (MQ=255)
                                                                                                 gCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTgg                                                     <  2:95137/148‑1 (MQ=255)
                                                                                                   aGCGTTACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTTGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGGAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGaa                                                  <  7:244886/149‑1 (MQ=255)
                                                                                                   aGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATcg                                                                            >  5:147139/1‑123 (MQ=255)
                                                                                                   aGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATcg                                                                            <  6:147139/123‑1 (MQ=255)
                                                                                                   aGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTggg                                                    >  5:13340/1‑147 (MQ=255)
                                                                                                   aGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGaa                                                  <  8:228212/149‑1 (MQ=255)
                                                                                                        gacctgaccGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCt                                                                                                                     >  1:21986/1‑77 (MQ=255)
                                                                                                        gacctgaccGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCt                                                                                                                     <  2:21986/77‑1 (MQ=255)
                                                                                                           ctgaccGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGgc                                          <  7:71750/149‑1 (MQ=255)
                                                                                                              accGACTCCACCGTGCTGCGCAACCTGGGCGTTGGCTTTGCCCACAGCGTCATCGCCTTCGATGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCTCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGgcc                                         <  7:257775/147‑1 (MQ=255)
                                                                                                                cGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGCAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGgccgagc                                     <  4:211772/149‑1 (MQ=255)
                                                                                                                  aCTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTgg                                                     >  1:213988/1‑131 (MQ=255)
                                                                                                                  aCTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTgg                                                     <  2:213988/131‑1 (MQ=255)
                                                                                                                      cACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCCCAGCGTCATCGCCTTCGCGGCCAGCCAGAAAGGCCTCGGCAAGCGGGAACTCAACGAAGCTAGCAAACCCGCCGCCCAGGACGCCAGCTGGGAAGTGCTGGGCGCGCCGATc                                >  2:327002/1‑148 (MQ=255)
                                                                                                                      cACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATcc                               <  3:284248/149‑1 (MQ=255)
                                                                                                                         cGTGCTGCGCAACCTGGGCGTGGGCTTTTCCCACAGCGTCATCGCCTTCCAAGCCAGCCTGAAACGCATCGGCAAGCTGGAAGTCAACGAAGCTCACAAAGCCCCCGACCTGGACGCCTGCTGGGAACTGCTGGCCCAGCCGATGCAGa                            >  6:32631/1‑149 (MQ=255)
                                                                                                                              tgcGCAACCTGGGCGTGGGGTTTTGCCACAGCGTCATCGCCGTCGAGGCCAGCCTGAAAGGGATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGCCCGGGCCGCCGGGGGGGAAGTGCTGGCCGGCCCGAGCCAGACCgg                        >  4:26350/1‑147 (MQ=255)
                                                                                                                              tgcGCAACCTGCGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCAGCCGACCTGGCCGCCTGCTGGGAAGTGCTGGCCGAGCGGATCCAGACAGTg                       >  4:247974/1‑149 (MQ=255)
                                                                                                                                   aaCCTGGGCGTGGGCTTTGGCCACAGCGTCATCGCCTTCGAGGCCAGCCGGAAAGGCATCGGCAACCTGGAAGTCAACGAAGCTCGCATCGCCGCCGCCCGGGACGCCTGCTGGGCAGTGCTGGCCGCGCCGCTCCCGACCGTGATgcg                  >  2:143979/1‑149 (MQ=255)
                                                                                                                                   aaCCTGGGCGTGGGCTTTGCCCACAGCGTCATCTCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGa                      <  7:114425/145‑1 (MQ=255)
                                                                                                                                     ccTGGGCGTGGGCTTTGGCCACAGCGTGATCGCCTTAGAGGCCAGCCGGAAATGCATCGGCAAGCTGGAAGTCAACGAAGCGAGCATAGCCGCCGACCTGGACGCCGGCTGGGAAGTGCTGGCCGCGCCGCTCCAGAGCGTGAATcgcc                >  4:19132/1‑149 (MQ=255)
                                                                                                                                           cGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCa          <  3:19132/149‑1 (MQ=255)
                                                                                                                                           cGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCa          <  2:274189/149‑1 (MQ=255)
                                                                                                                                            gTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCaa         <  2:303347/149‑1 (MQ=255)
                                                                                                                                                gCTTTGCCCAGAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATc     <  7:88988/149‑1 (MQ=255)
                                                                                                                                                    tGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGaa  <  1:44990/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGTCAACCCGATCGACTTCGAAAACTCCGAAGGCAACCTGGGTATCGCCAACGCGCTGTTCCAGCACCTGGCCAGCAAGCTGCCGATCTCGCGCTGGCAGCGTGACCTGACCGACTCCACCGTGCTGCGCAACCTGGGCGTGGGCTTTGCCCACAGCGTCATCGCCTTCGAGGCCAGCCTGAAAGGCATCGGCAAGCTGGAAGTCAACGAAGCTCGCATCGCCGCCGACCTGGACGCCTGCTGGGAAGTGCTGGCCGAGCCGATCCAGACCGTGATGCGCCGCTTCAACATCGAA  >  NC_002947/4527859‑4528154

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