Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,015,893 Δ1 bp 100% intergenic (‑156/+78) PP_5679 ← / ← PP_4421 hypothetical protein/aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,015,8930A.88.2% 236.1 / 27.0 68intergenic (‑156/+78)PP_5679/PP_4421hypothetical protein/aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base . (32/28);  total (36/32)
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.39e-01
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

CTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCT  >  NC_002947/5015747‑5016041
                                                                                                                                                  |                                                                                                                                                    
ctcGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaa                                                                                                                                                     >  6:168643/1‑148 (MQ=255)
    tttGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaa                                                                                                                                                 >  1:274649/1‑148 (MQ=255)
      tGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaa                                                                                                                                               <  5:122469/148‑1 (MQ=255)
       gAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACaaaaaaaaag                                                                                                                                             >  1:16473/1‑148 (MQ=255)
                     atgaATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTCCGGCCTTGTGACCTGATGACAGCGTCTTGGCGTTGAGTGTGATAGCACTACAAAAAAAAAAGCCCGCCGAAAGGc                                                                                                                              >  3:243385/1‑149 (MQ=255)
                        taaaaGCCCGCGCCGGAGACCGGTGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGc                                                                                                                          >  4:95901/4‑149 (MQ=255)
                        aaTAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCggg                                                                                                                           <  5:168643/148‑1 (MQ=255)
                        aaTAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCggg                                                                                                                           <  5:58093/148‑1 (MQ=255)
                           aaGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCtt                                                                                                                        >  8:197442/1‑148 (MQ=255)
                             gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAATTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTTTAGGCGGGCttt                                                                                                                       <  4:270637/149‑2 (MQ=255)
                             gCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTaa                                                                                                                      >  2:93316/1‑148 (MQ=255)
                              cccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa                                                                                                                    <  6:339958/149‑1 (MQ=255)
                               ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTTAAAc                                                                                                                   >  7:292502/1‑149 (MQ=255)
                               ccGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGGCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTaaaa                                                                                                                    <  4:220592/149‑1 (MQ=255)
                                  cACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGaa                                                                                                                <  2:274649/149‑1 (MQ=255)
                                   aCCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAACGCGGGCTTAAAACg                                                                                                                  >  3:229607/1‑146 (MQ=255)
                                             ggCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCCGCAAAAAGCGGGCATAAAACGAAGAGGGGAAt                                                                                                       >  2:317100/1‑148 (MQ=255)
                                                gcgGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAg                                                                                                   <  1:72629/149‑1 (MQ=255)
                                                  gggCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTg                                                                                                >  4:281814/1‑149 (MQ=255)
                                                    gCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTgg                                                                                               >  7:109594/1‑148 (MQ=255)
                                                      ttttACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcat                                                                                            >  1:292111/1‑149 (MQ=255)
                                                           cTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAACAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                        <  1:93316/149‑1 (MQ=255)
                                                           cTTGCCGGGCCTTGGACCGTATCGGCCCCGAAAGTACGTCCGTGTGAACTGAGGACAGCGGCTTGGCTTTGATTGTGATCGCACTACAAAAAAAAAAGCCCGCCTACAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                        <  4:243385/149‑1 (MQ=255)
                                                              tCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGc                                                                                      <  8:109594/147‑1 (MQ=255)
                                                               cggggCCTTGCACTGTATCTCCCCCTAATTTACGTCTTTGTTAACTCTGCACTGCGCCTTGCGTTTGATTGTTATCGCACTCAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGcc                                                                                     <  1:317100/146‑1 (MQ=255)
                                                                 gggCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAAc                                                                                  >  3:306166/1‑149 (MQ=255)
                                                                    ccTTGGACAGTATCTGACCCGAAAGTACGTCCTTTTTAACTGATGACAGCGTCTTGGCTTTGATTTTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGtt                                                                               <  3:160503/149‑1 (MQ=255)
                                                                      ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                        <  2:206423/137‑1 (MQ=255)
                                                                      ttGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGCCTTAAAACGAAGAGGGTAATTCAGCTGGCATCAtt                                                                                        >  1:206423/1‑137 (MQ=255)
                                                                         gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCa                                                                                                    <  7:340065/122‑1 (MQ=255)
                                                                         gACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCa                                                                                                    >  8:340065/1‑122 (MQ=255)
                                                                                   gACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCt                                                                 >  2:363387/1‑148 (MQ=255)
                                                                                    aCCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCtt                                                                <  3:95901/147‑1 (MQ=255)
                                                                                        gAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTTATGGCTCAGACCTTCAGc                                                            >  7:328104/1‑148 (MQ=255)
                                                                                         aaaGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtc                                                          >  2:178109/1‑148 (MQ=255)
                                                                                             tACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCtcatc                                                       >  8:278593/1‑147 (MQ=255)
                                                                                                gTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAgg                                                   <  4:306166/149‑1 (MQ=255)
                                                                                                gTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGa                                                  >  5:25686/1‑149 (MQ=255)
                                                                                                 tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAaga                                                                                                              <  6:5742/88‑1 (MQ=255)
                                                                                                 tCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAaga                                                                                                              >  5:5742/1‑88 (MQ=255)
                                                                                                    ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGATGGCATCATAACCAACGTTTATGGCTCATACCTTAAGCTAATCAAGGACg                                                >  1:119031/1‑149 (MQ=255)
                                                                                                    ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                            <  1:202072/119‑1 (MQ=255)
                                                                                                    ttGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGAt                                                                            >  2:202072/1‑119 (MQ=255)
                                                                                                      gtgAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTc                                             <  1:178109/148‑1 (MQ=255)
                                                                                                          aCTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTACAGCTCATCAAGGACGCAAAGca                                         >  6:194602/1‑148 (MQ=255)
                                                                                                              atgaCAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCt                                    <  2:16473/149‑1 (MQ=255)
                                                                                                                gaCAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCCGCACAGCt                                    >  3:50119/1‑149 (MQ=255)
                                                                                                                  cAGCGTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                <  2:292111/149‑1 (MQ=255)
                                                                                                                   aGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCgt                                <  1:363387/149‑1 (MQ=255)
                                                                                                                      gTCTTGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGc                            <  6:25686/149‑1 (MQ=255)
                                                                                                                          tGGCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATg                        >  1:348187/1‑149 (MQ=255)
                                                                                                                           ggCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATg                        <  2:40263/148‑1 (MQ=255)
                                                                                                                           ggCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTAAGACCTTCAGCTCATCAAGGACGCACAGCACAGCTTCGTGTGCGATGc                       >  6:221228/1‑149 (MQ=255)
                                                                                                                            gCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATg                        <  4:50119/149‑1 (MQ=255)
                                                                                                                            gCTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGc                       <  7:278593/148‑1 (MQ=255)
                                                                                                                             cTTTGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCt                      >  6:243895/1‑148 (MQ=255)
                                                                                                                              tttGATTGTGATAGCACTACAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCt                      >  1:101523/1‑148 (MQ=255)
                                                                                                                              tttGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcatcat                                                                                         >  3:216057/1‑80 (MQ=255)
                                                                                                                              tttGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGcatcat                                                                                         <  4:216057/80‑1 (MQ=255)
                                                                                                                              tttGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtga                    >  1:99089/1‑149 (MQ=255)
                                                                                                                              tttGATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGAGTGCGATGCtga                    >  1:195445/1‑149 (MQ=255)
                                                                                                                                 gATTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgatg                  <  3:281814/148‑1 (MQ=255)
                                                                                                                                  aTTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCtgatga                 <  2:99089/148‑1 (MQ=255)
                                                                                                                                  aTTGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTgatgat                >  6:256603/1‑149 (MQ=255)
                                                                                                                                   ttGTGATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTgatgat                <  5:124651/148‑1 (MQ=255)
                                                                                                                                      tgATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGa            <  5:243895/149‑1 (MQ=255)
                                                                                                                                       gATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCGTAAAACGAAGAGGGTAATTAAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCCCAGCTTAGTGTGCGATGCTGATGATGc              >  2:281469/1‑148 (MQ=255)
                                                                                                                                       gATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGAt           >  4:42638/1‑149 (MQ=255)
                                                                                                                                       gATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGAt           >  6:249889/1‑149 (MQ=255)
                                                                                                                                       gATAGCACTAC‑‑AAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGAt           >  2:94380/1‑149 (MQ=255)
                                                                                                                                               aaaaaaaaaaaaaaGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACtt     <  1:318317/146‑1 (MQ=255)
                                                                                                                                                 cAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg                                                    <  3:185669/100‑1 (MQ=255)
                                                                                                                                                 cAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAg                                                    >  4:185669/1‑100 (MQ=255)
                                                                                                                                                  aaaaaaaaaaaGCCCGCCTAAAGGCGGGGTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCATCGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCt  >  7:274132/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTCGTTTGAGTCCACGATGAGATGAATAAGCCCGCACCGGTCTCCGGCGCGGGCTTTTACTTTCCGGGCCTTGGACAGTATCTGACCCGAAAGTACGTCCTTGTGAACTGATGACAGCGTCTTGGCTTTGATTGTGATAGCACTACAAAAAAAAAAAGCCCGCCTAAAGGCGGGCTTAAAACGAAGAGGGTAATTCAGCTGGCATCATTGCCAACGTTGATGGCTCAGACCTTCAGCTCATCAAGGACGCTCAGCACAGCTTCGTGTGCGATGCTGATGATGCGATCCACTTCCT  >  NC_002947/5015747‑5016041

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: