Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,903,078 C→A 7.9% intergenic (+43/+72) PP_3445 → / ← ilvA‑I hypothetical protein/threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,903,0780CA 7.9% 232.1 / 8.1 76intergenic (+43/+72)PP_3445/ilvA‑Ihypothetical protein/threonine deaminase
Reads supporting (aligned to +/- strand):  ref base C (32/38);  new base A (3/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 = 8.44e-01

GGCCCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902934‑3903224
                                                                                                                                                |                                                                                                                                                   
ggCCCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTcgggc                                                                                                                                                 <  4:72483/149‑1 (MQ=255)
ggCCCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTcgggc                                                                                                                                                 >  5:176002/1‑149 (MQ=255)
   cccGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTcgggcggg                                                                                                                                              <  8:16448/149‑1 (MQ=255)
    ccGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCgggcgggg                                                                                                                                             >  3:142798/1‑149 (MQ=255)
       aGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg                                                                                                                                          <  2:128235/149‑1 (MQ=255)
           ggATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        >  5:317013/1‑147 (MQ=255)
           ggATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        <  6:317013/147‑1 (MQ=255)
            gATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                                     >  4:42958/1‑149 (MQ=255)
                ccGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  6:193957/148‑1 (MQ=255)
                ccGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  4:199251/148‑1 (MQ=255)
                 cGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGc                                                                                                                                 >  4:275544/1‑148 (MQ=255)
                   aggaCTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGaa                                                                                                                              <  4:205306/146‑1 (MQ=255)
                       cTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAAcaca                                                                                                                          <  7:47352/149‑1 (MQ=255)
                       cTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        >  1:245314/1‑135 (MQ=255)
                       cTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        <  2:245314/135‑1 (MQ=255)
                        ttCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         >  6:360824/1‑149 (MQ=255)
                        ttCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         >  8:178472/1‑149 (MQ=255)
                        ttCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         >  2:284754/1‑149 (MQ=255)
                          cAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         <  3:192438/147‑1 (MQ=255)
                            aGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcgg                                                                                                                    >  3:36129/1‑149 (MQ=255)
                                cTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                                     >  5:160388/1‑129 (MQ=255)
                                cTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                                     <  6:160388/129‑1 (MQ=255)
                                   cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcg                                                                                                              <  3:275544/148‑1 (MQ=255)
                                    gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCa                                                                                                                                       <  8:287680/123‑1 (MQ=255)
                                    gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCa                                                                                                                                       >  7:287680/1‑123 (MQ=255)
                                    gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcag                                                                                                                 >  3:182819/1‑144 (MQ=255)
                                    gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcag                                                                                                                 <  4:182819/144‑1 (MQ=255)
                                    gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACaggcggcaggcggc                                                                                                            >  2:286553/1‑149 (MQ=255)
                                            cAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAgg                                                                                                    >  6:272613/1‑149 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCTCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAgg                                                                                                    <  8:95486/148‑1 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAaca                                                                                                                            >  5:193019/1‑125 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAaca                                                                                                                            <  6:193019/125‑1 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  3:89335/1‑149 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  5:245295/1‑149 (MQ=255)
                                             aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   <  8:218462/149‑1 (MQ=255)
                                                      gAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTgg                                                                                          >  3:26883/1‑149 (MQ=255)
                                                       aaaCGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGc                                                                                         <  1:65205/149‑1 (MQ=255)
                                                       aaaCGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGc                                                                                         <  6:176002/149‑1 (MQ=255)
                                                           gTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGca                                                                                     >  1:271730/1‑149 (MQ=255)
                                                                             cGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  7:245574/87‑1 (MQ=255)
                                                                             cGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  >  8:245574/1‑87 (MQ=255)
                                                                             cGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAg                                                                    >  1:343930/1‑148 (MQ=255)
                                                                              ggCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCTGAAGCCCAgg                                                                   >  1:123826/1‑148 (MQ=255)
                                                                                cTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAAc                                                                <  4:142798/149‑1 (MQ=255)
                                                                                    ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  1:284754/148‑1 (MQ=255)
                                                                                    ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  3:42958/148‑1 (MQ=255)
                                                                                    ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  7:178472/148‑1 (MQ=255)
                                                                                      cATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCtt                                                           <  4:36129/148‑1 (MQ=255)
                                                                                      cATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCtt                                                           <  5:360824/148‑1 (MQ=255)
                                                                                          cgcgAGCGTTTAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGCCTGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTAc                                                       <  1:286553/149‑1 (MQ=255)
                                                                                           gcgAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGAcgg                                                                                                     >  7:166808/1‑99 (MQ=255)
                                                                                           gcgAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGAcgg                                                                                                     <  8:166808/101‑3 (MQ=255)
                                                                                           gcgAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGc                                                     >  3:54579/1‑149 (MQ=255)
                                                                                             gAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCCTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAg                                                                         >  2:200895/1‑127 (MQ=255)
                                                                                             gAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCCTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAg                                                                         <  1:200895/127‑1 (MQ=255)
                                                                                                cGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGt                                                <  4:54579/149‑1 (MQ=255)
                                                                                                 gTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               >  6:140450/1‑149 (MQ=255)
                                                                                                 gTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               <  6:245295/149‑1 (MQ=255)
                                                                                                 gTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               >  1:14816/1‑149 (MQ=255)
                                                                                                    gAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGt                                             <  2:14816/148‑1 (MQ=255)
                                                                                                       aTTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGt                                         >  7:58157/1‑149 (MQ=255)
                                                                                                            cTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGt                                    <  5:272613/149‑1 (MQ=255)
                                                                                                                    gTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTAc                            <  2:271730/149‑1 (MQ=255)
                                                                                                                        ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTTAGCCCAGGAACAGCTTGTACGCCGGGTTGt                                             <  5:104744/128‑1 (MQ=255)
                                                                                                                        ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTTAGCCCAGGAACAGCTTGTACGCCGGGTTGt                                             >  6:104744/1‑128 (MQ=255)
                                                                                                                        ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCgg                         >  3:142260/1‑148 (MQ=255)
                                                                                                                        ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCgg                         <  4:26883/148‑1 (MQ=255)
                                                                                                                             aGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTAcccc                   <  5:140450/149‑1 (MQ=255)
                                                                                                                              ggAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCa                  <  2:343930/149‑1 (MQ=255)
                                                                                                                               gagaTCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCAt                 <  8:58157/149‑1 (MQ=255)
                                                                                                                                 gaccgccccgcccGAAGGGCAGGGCGATCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGt                                                         >  1:65364/4‑107 (MQ=38)
                                                                                                                                 gaccgccccgcccGAAGGGCAGGGCGATCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGt                                                         <  2:65364/104‑1 (MQ=38)
                                                                                                                                      ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               >  2:134182/1‑112 (MQ=255)
                                                                                                                                      ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               <  1:134182/112‑1 (MQ=255)
                                                                                                                                       cctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  7:40949/1‑59 (MQ=255)
                                                                                                                                       cctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   <  8:40949/59‑1 (MQ=255)
                                                                                                                                              ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc                                  <  2:120667/117‑1 (MQ=255)
                                                                                                                                              ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc                                  >  1:120667/1‑117 (MQ=255)
                                                                                                                                              ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCCGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGACCAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGa  >  7:119121/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GGCCCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902934‑3903224

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: