Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,344,688:1 +C 5.5% intergenic (‑39/+36) PP_2061 ← / ← PP_2062 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,344,6881.C 5.5% 293.2 / 8.9 72intergenic (‑39/+36)PP_2061/PP_2062hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (38/30);  new base C (2/2);  total (40/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.43e-01

AGTCGTCGACAATAAATACCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAGGTGTTTG  >  NC_002947/2344541‑2344835
                                                                                                                                                     |                                                                                                                                                   
aGTCGTCGACAATAAATACCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGtt                                                                                                                                                    >  4:420679/1‑149 (MQ=255)
  tcgtcgACAATAAATACCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTtc                                                                                                                                                  <  5:334452/149‑1 (MQ=255)
       gACAATAAATACCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTtcgtcg                                                                                                                                              <  5:438561/148‑1 (MQ=255)
          aataaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCtt                                                                                                                                          >  4:367660/1‑149 (MQ=255)
          aataaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCt                                                                                                                                           >  4:29853/1‑148 (MQ=255)
          aataaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCt                                                                                                                                           >  3:172998/1‑148 (MQ=255)
           ataaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCtt                                                                                                                                          >  5:389557/1‑148 (MQ=255)
           ataaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTc                                                                                                                                         <  6:214476/149‑1 (MQ=255)
            taaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCa                                                                                                                                        >  8:110300/1‑149 (MQ=255)
             aaataCCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCACCTTCAt                                                                                                                                       >  8:375632/1‑149 (MQ=255)
                  ccccGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATg                                                                                                                                  <  3:367660/149‑1 (MQ=255)
                   cccGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTAt                                                                                                                                   >  3:165483/1‑147 (MQ=255)
                   cccGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGa                                                                                                                                 <  7:172032/149‑1 (MQ=255)
                   cccGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGa                                                                                                                                 <  7:79633/149‑1 (MQ=255)
                    ccGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAg                                                                                                                                >  5:209059/1‑149 (MQ=255)
                          tcATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTa                                                                                                                          >  6:196005/1‑149 (MQ=255)
                             tAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTAttt                                                                                                                       >  4:128147/1‑149 (MQ=255)
                             tAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTAttt                                                                                                                       <  6:237955/149‑1 (MQ=255)
                              aGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTg                                                                                                                      >  2:34422/1‑149 (MQ=255)
                                gCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGt                                                                                                                     >  4:357599/1‑148 (MQ=255)
                                gCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTa                                                                                                                    <  1:34422/149‑1 (MQ=255)
                                      taACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACCACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCa                                                                                                               <  7:110300/148‑1 (MQ=255)
                                       aaCCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGt                                                                                                             >  7:321495/1‑149 (MQ=255)
                                        aCCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCGGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGt                                                                                                             >  1:108025/1‑148 (MQ=255)
                                        aCCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGt                                                                                                             >  2:70615/1‑148 (MQ=255)
                                        aCCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGt                                                                                                             <  4:172998/148‑1 (MQ=255)
                                        aCCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTg                                                                                                            <  3:420679/149‑1 (MQ=255)
                                             cAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATggg                                                                                                       <  2:363988/149‑1 (MQ=255)
                                                 gcTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGa                                                                                                   >  5:2950/1‑149 (MQ=255)
                                                 gcTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGa                                                                                                   >  5:9019/1‑149 (MQ=255)
                                                 gcTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGa                                                                                                   >  4:184825/1‑149 (MQ=255)
                                                        ttATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGAttttttt                                                                                            <  6:209059/149‑1 (MQ=255)
                                                           ttGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTgcg                                                                                         <  3:357599/149‑1 (MQ=255)
                                                               tGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTAt                                                                                     >  1:418376/1‑149 (MQ=255)
                                                                 gCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTAtt                                                                                    <  1:70615/148‑1 (MQ=255)
                                                                   aTTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTAttgtt                                                                                 >  1:335696/1‑149 (MQ=255)
                                                                           tCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGtt                                                                           >  1:102428/1‑148 (MQ=255)
                                                                            ccAGATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTCGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGa                                                                         <  4:253637/148‑1 (MQ=255)
                                                                              acaTTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGataa                                                                      >  8:236654/1‑149 (MQ=255)
                                                                                    cAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTg                                                                 <  2:418376/148‑1 (MQ=255)
                                                                                     aGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGAt                                                               <  5:378467/149‑1 (MQ=255)
                                                                                     aGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGAt                                                               >  5:404973/1‑149 (MQ=255)
                                                                                           gcTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTc                                                                              <  6:40858/128‑1 (MQ=255)
                                                                                           gcTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTc                                                                              >  5:40858/1‑128 (MQ=255)
                                                                                               gCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATg                                                     >  6:184730/1‑149 (MQ=255)
                                                                                                cGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGc                                                    <  3:184825/149‑1 (MQ=255)
                                                                                                 gTGAGTATTCAAGTGGGCGTCAGTAGAGGCGCCGAAGACGACCTGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGcc                                                   <  5:424601/149‑1 (MQ=255)
                                                                                                 gTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGc                                                    >  7:310164/1‑148 (MQ=255)
                                                                                                  tGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTAATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATATTATTGATGTCTACGATTccc                                                  >  5:353985/1‑149 (MQ=255)
                                                                                                      tATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGc                                              <  5:184730/149‑1 (MQ=255)
                                                                                                         tCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTta                                           >  2:115443/1‑149 (MQ=255)
                                                                                                           aaGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTtat                                          >  3:227238/1‑148 (MQ=255)
                                                                                                           aaGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCACCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTtat                                          <  7:375632/148‑1 (MQ=255)
                                                                                                            agagGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATaa                                        <  3:29853/149‑1 (MQ=255)
                                                                                                             gagGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTtata                                         >  1:193080/1‑147 (MQ=255)
                                                                                                                 gCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCtt                                   <  4:165483/149‑1 (MQ=255)
                                                                                                                   gTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCtttt                                 >  1:106979/1‑149 (MQ=255)
                                                                                                                   gTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCtttt                                 <  2:106979/149‑1 (MQ=255)
                                                                                                                       gTAGAGGCGCCGAAGACGACCAGGGTCGTTTGGTGGCCTTCATTTATGAGTGTGTATTTTTAAATAAGTGATGGGCAGATTTTTTTGGGCTTTTGTTGTGGTTGATAATATTGATGTCTACGATTCCCTTGATTATAAAGCttattgc                              >  6:293912/1‑143 (MQ=255)
                                                                                                                        tgaagGCGACGAAACGACCCTGGTCGTCTTCGGCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCAt                            >  6:420743/4‑149 (MQ=38)
                                                                                                                        tAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCAt                            >  7:180528/1‑149 (MQ=255)
                                                                                                                                cgacgaaacgacCCTGGTCGTCTTCGGCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCtt                        >  2:344698/8‑147 (MQ=255)
                                                                                                                                    aaGACGACCAGGGTCGTTTCGTCGACAACATTTATGAGTGTGTATTTGTAAATCAGTGATGCGATGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTAATAAAGCTTTTGGCATCATTGAGTTCg                 >  8:187130/1‑149 (MQ=255)
                                                                                                                                       acgacCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGtt              <  2:108025/149‑1 (MQ=255)
                                                                                                                                       acgacCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGtt              >  2:148111/1‑149 (MQ=255)
                                                                                                                                        cgacCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGtt              <  6:389557/148‑1 (MQ=255)
                                                                                                                                         gacCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGttt             >  6:253050/1‑148 (MQ=255)
                                                                                                                                           ccctggTCGTCTTCGGCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTa           <  1:344698/145‑1 (MQ=255)
                                                                                                                                            cctggTCGTCTTCGGCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTaa          <  5:420743/146‑1 (MQ=255)
                                                                                                                                            cAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGATTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAg         <  6:9019/149‑1 (MQ=255)
                                                                                                                                            cAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGATTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAg         <  6:2950/149‑1 (MQ=255)
                                                                                                                                                 tCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAGGTGtt    >  4:180601/1‑149 (MQ=255)
                                                                                                                                                   gTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAGGTGTTTg  >  2:102538/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGTCGTCGACAATAAATACCCCGTTCTCATAGGCGATATAACCCGCAGCGCTAACGTTATTGCTGGCATTGAAATTCCACATTTCAGCGCCGCTGGCGTGAGTATTCAAGAGGGCGTCAGTAGAGGCGCCGAAGACGACCAGGGTCGTTTCGTCGCCTTCATTTATGAGTGTGTATTTGTAAATCAGTGATGGGCTGATTTTTTTGCGCTATTGTTGTCGTTGATAATATTGATGTCTACGATGCCCTTGCTTATAAAGCTTTTGGCATCCTTGAGTTCGGTTTTAAGGTGTTTG  >  NC_002947/2344541‑2344835

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: