Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,903,089 G→A 8.1% intergenic (+54/+61) 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,0890GA 8.1% 215.2 / 6.5 74intergenic (+54/+61)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 = 8.90e-01

AGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902941‑3903224
                                                                                                                                                    |                                                                                                                                        
aGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGggcggggcgg                                                                                                                                          <  2:128235/149‑1 (MQ=255)
    ggATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        <  6:317013/147‑1 (MQ=255)
    ggATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGACCGCCCCGCCCGAAGGGCAGGGCGATc                                                                                                                                        >  5:317013/1‑147 (MQ=255)
     gATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACg                                                                                                                                     >  4:42958/1‑149 (MQ=255)
         ccGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  4:199251/148‑1 (MQ=255)
         ccGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAg                                                                                                                                  <  6:193957/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                                                                                                                         >  2:284754/1‑149 (MQ=255)
                 ttCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         >  8:178472/1‑149 (MQ=255)
                 ttCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACAcag                                                                                                                         >  6:360824/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                                                                                                                                       >  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                                                                                                                            >  5:193019/1‑125 (MQ=255)
                                      aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAaca                                                                                                                            <  6:193019/125‑1 (MQ=255)
                                      aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  5:245295/1‑149 (MQ=255)
                                      aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   <  8:218462/149‑1 (MQ=255)
                                      aaaaGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGc                                                                                                   >  3:89335/1‑149 (MQ=255)
                                               gAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTgg                                                                                          >  3:26883/1‑149 (MQ=255)
                                                aaaCGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGc                                                                                         <  6:176002/149‑1 (MQ=255)
                                                aaaCGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGc                                                                                         <  1:65205/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                                                             <  1:284754/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                                                                         <  1:200895/127‑1 (MQ=255)
                                                                                      gAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCCTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAg                                                                         >  2:200895/1‑127 (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                                             <  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)
                                                                                                                                            gcggggcggTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGtctc                                      <  4:331193/108‑1 (MQ=255)
                                                                                                                                            gcggggcggTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGtctc                                      >  3:331193/1‑108 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
AGCCGGATACCGAGAACTTCAAGCGCTACGCCGACTTCAAAAGCCTGGAAACGTTCCGCAAGTTCGTTGCCGGCTTTTTCATGCGCGAGCGTTGAAATTGGCTCGCGCAGTTCCCTGCAGGAGATCGCCCTGCCCTTCGGGCGGGGCGGTCACGCAGCGAACACAGGCGGCAGGCGGCTAGACAGGCGTGCATTGGCCGCACAGGCGCGTCAGCCCAGGAACAGCTTGTACGCCGGGTTGTCAGTCTCGTCCCAGTACCGGTACCCCATCTCGTCCAGCGTCGA  >  NC_002947/3902941‑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: