Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,903,083 G→A 8.1% intergenic (+48/+67) 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,0830GA 8.1% 207.3 / 5.9 74intergenic (+48/+67)PP_3445/ilvA‑Ihypothetical protein/threonine deaminase
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base A (3/3);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

CCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902937‑3903224
                                                                                                                                                  |                                                                                                                                              
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                                                                                                                                     <  6:160388/129‑1 (MQ=255)
                             cTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                                     >  5:160388/1‑129 (MQ=255)
                                cGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcaggcg                                                                                                              <  3:275544/148‑1 (MQ=255)
                                 gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCa                                                                                                                                       >  7:287680/1‑123 (MQ=255)
                                 gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCa                                                                                                                                       <  8:287680/123‑1 (MQ=255)
                                 gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcag                                                                                                                 <  4:182819/144‑1 (MQ=255)
                                 gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcaggcggcag                                                                                                                 >  3:182819/1‑144 (MQ=255)
                                 gCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACaggcggcaggcggc                                                                                                            >  2:286553/1‑149 (MQ=255)
                                         cAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAgg                                                                                                    >  6:272613/1‑149 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCTCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAgg                                                                                                    <  8:95486/148‑1 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAaca                                                                                                                            <  6:193019/125‑1 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAaca                                                                                                                            >  5:193019/1‑125 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   <  8:218462/149‑1 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  3:89335/1‑149 (MQ=255)
                                          aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  5:245295/1‑149 (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                                                                                                                                  >  8:245574/1‑87 (MQ=255)
                                                                          cGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  7:245574/87‑1 (MQ=255)
                                                                          cGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAg                                                                    >  1:343930/1‑148 (MQ=255)
                                                                           ggCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCTGAAGCCCAgg                                                                   >  1:123826/1‑148 (MQ=255)
                                                                             cTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAAc                                                                <  4:142798/149‑1 (MQ=255)
                                                                                 ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  3:42958/148‑1 (MQ=255)
                                                                                 ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  7:178472/148‑1 (MQ=255)
                                                                                 ttCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGc                                                             <  1:284754/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                                               >  1:14816/1‑149 (MQ=255)
                                                                                              gTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               <  6:245295/149‑1 (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                                             >  6:104744/1‑128 (MQ=255)
                                                                                                                     ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTTAGCCCAGGAACAGCTTGTACGCCGGGTTGt                                             <  5:104744/128‑1 (MQ=255)
                                                                                                                     ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCgg                         <  4:26883/148‑1 (MQ=255)
                                                                                                                     ccTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCgg                         >  3:142260/1‑148 (MQ=255)
                                                                                                                          aGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTAcccc                   <  5:140450/149‑1 (MQ=255)
                                                                                                                           ggAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCa                  <  2:343930/149‑1 (MQ=255)
                                                                                                                            gagaTCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCAt                 <  8:58157/149‑1 (MQ=255)
                                                                                                                              gaccgccccgcccGAAGGGCAGGGCGATCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGt                                                         <  2:65364/104‑1 (MQ=38)
                                                                                                                              gaccgccccgcccGAAGGGCAGGGCGATCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGt                                                         >  1:65364/4‑107 (MQ=38)
                                                                                                                                   ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               <  1:134182/112‑1 (MQ=255)
                                                                                                                                   ccctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGtt                                               >  2:134182/1‑112 (MQ=255)
                                                                                                                                    cctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  7:40949/1‑59 (MQ=255)
                                                                                                                                    cctgccctTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   <  8:40949/59‑1 (MQ=255)
                                                                                                                                           ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc                                  >  1:120667/1‑117 (MQ=255)
                                                                                                                                           ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTcc                                  <  2:120667/117‑1 (MQ=255)
                                                                                                                                           ttCGGGCGGGGCGGTCACGCAGCGAACACAGGCCGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGACCAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGa  >  7:119121/1‑149 (MQ=255)
                                                                                                                                                gcggggcggTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGtctc                                      <  4:331193/108‑1 (MQ=255)
                                                                                                                                                gcggggcggTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGtctc                                      >  3:331193/1‑108 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CCCGAGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902937‑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: