Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,223 A→C 7.0% T8P (ACC→CCC)  PP_3559 → glycine betaine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,038,2230AC 7.0% 204.3 / 5.0 68T8P (ACC→CCC) PP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (28/35);  new base C (3/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCCGTCATGAGCGAAAAACGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGG  >  NC_002947/4038080‑4038369
                                                                                                                                               |                                                                                                                                                  
gCCGTCATGAGCGAAAAACGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTg                                                                                                                                               <  2:2569/149‑1 (MQ=255)
  cGTCATGAGCGAAAAACGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCa                                                                                                                                             <  2:149371/149‑1 (MQ=255)
       tGAGCGAAAAACGCACGCCCCCGCTCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGtttt                                                                                                                                        <  1:149633/149‑1 (MQ=255)
         aGCGAAAAACGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTc                                                                                                                                       <  4:174377/148‑1 (MQ=255)
              aaaaCGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGc                                                                                                                                 <  1:44872/149‑1 (MQ=255)
                     acgcCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGGGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTCTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAg                                                                                                                           >  2:181238/1‑148 (MQ=255)
                      cgcCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCa                                                                                                                         <  5:324822/149‑1 (MQ=255)
                                 gAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGCCCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTg                                                                                                              >  8:99989/1‑149 (MQ=255)
                                             aCCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggttg                                                                                                  <  4:346509/149‑1 (MQ=255)
                                             aCCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggttg                                                                                                  >  6:75008/1‑149 (MQ=255)
                                             aCCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggtt                                                                                                   <  6:279056/148‑1 (MQ=255)
                                              ccTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggttg                                                                                                  <  5:75008/148‑1 (MQ=255)
                                                    tGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTg                                                                                            >  8:151485/1‑148 (MQ=255)
                                                      cgcgATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAAc                                                                                         <  8:79168/149‑1 (MQ=255)
                                                         gATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGCCCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTAc                                                                                      >  8:84181/1‑149 (MQ=255)
                                                          aTAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACg                                                                                     <  1:130694/149‑1 (MQ=255)
                                                             aTCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGCCCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGc                                                                                   <  7:84181/148‑1 (MQ=255)
                                                              tCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGGATTGAACTACGGCGa                                                                                 >  2:233821/1‑149 (MQ=255)
                                                                gcgAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGTTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGTTTGGTATTGAACTACGGCGacc                                                                               <  5:53142/149‑1 (MQ=255)
                                                                    ggTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggtt                                                                                                   >  5:316346/1‑125 (MQ=255)
                                                                    ggTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggtt                                                                                                   <  6:316346/125‑1 (MQ=255)
                                                                     gtgtGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTg                                                                          >  7:172745/1‑149 (MQ=255)
                                                                     gtgtGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTg                                                                          >  5:76003/1‑149 (MQ=255)
                                                                        tgGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTgcg                                                                        >  3:158436/1‑148 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGCCCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGaccac                                                                             <  6:240318/141‑1 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGCCCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGaccac                                                                             >  5:240318/1‑141 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggttggt                                                                                                >  5:202476/1‑122 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACtggttggt                                                                                                <  6:202476/122‑1 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGc                                                                                   >  5:202867/1‑135 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGc                                                                                   <  6:202867/135‑1 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAg                                                                     >  5:248250/1‑149 (MQ=255)
                                                                          gAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAg                                                                     >  5:51409/1‑149 (MQ=255)
                                                                               cgtGGCGGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGGGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAAGTACGGCGACCACCTGCGGTGTTTGt                                                                <  6:200106/147‑1 (MQ=255)
                                                                               ccTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGt                                                                <  4:158436/149‑1 (MQ=255)
                                                                               ccTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGt                                                                <  6:76003/149‑1 (MQ=255)
                                                                               ccTGGCGGCCGGCGGACCTCGCCGCCAAGGTCGACGCCAGCCTGTCCGCAGCGGTTTCCCGCAGACCCCCCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTCAACTACGGCGACCACCAGCCCAGTTTGt                                                                <  8:172745/149‑1 (MQ=255)
                                                                                cTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTc                                                               >  6:106203/1‑149 (MQ=255)
                                                                                   gCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCt                                                              <  6:51409/147‑1 (MQ=255)
                                                                                        cggcggACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGac                                                                                >  8:34675/1‑124 (MQ=255)
                                                                                        cggcggACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGac                                                                                <  7:34675/124‑1 (MQ=255)
                                                                                          gcggACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGAc                                                           <  3:159418/143‑1 (MQ=255)
                                                                                          gcggACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGAc                                                           >  4:159418/1‑143 (MQ=255)
                                                                                               cATCGCC‑CCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCaa                                               >  3:110551/1‑149 (MQ=255)
                                                                                                 tCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACg                                                                                     >  5:258319/1‑110 (MQ=255)
                                                                                                 tCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACg                                                                                     <  6:258319/110‑1 (MQ=255)
                                                                                                    ccgccAAGGTCGAGGCCAGCCTGTGAGCAGCGGGTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAActgc                                           <  5:90523/149‑1 (MQ=255)
                                                                                                     cgccAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAActgct                                          <  7:289721/149‑1 (MQ=255)
                                                                                                        cAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGt                                       >  1:51413/1‑149 (MQ=255)
                                                                                                        cAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGt                                       >  3:284622/1‑149 (MQ=255)
                                                                                                            gTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGCAGTTGCAGCAACTGGTcaga                                   <  7:151485/149‑4 (MQ=255)
                                                                                                               gAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCGACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTtgctg                                                  <  6:353255/131‑1 (MQ=255)
                                                                                                               gAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCGACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTtgctg                                                  >  5:353255/1‑131 (MQ=255)
                                                                                                                aGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTgg                               <  2:51413/149‑1 (MQ=255)
                                                                                                                  gCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAActgc                                           <  2:249677/135‑1 (MQ=255)
                                                                                                                  gCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAActgc                                           >  1:249677/1‑135 (MQ=255)
                                                                                                                    cAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGaa                           >  2:194573/1‑149 (MQ=255)
                                                                                                                    cAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACACCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGaa                           >  6:161446/1‑149 (MQ=255)
                                                                                                                     aGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAAc                          <  4:110551/149‑1 (MQ=255)
                                                                                                                     aGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAAc                          <  6:248250/149‑1 (MQ=255)
                                                                                                                         tgtgAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGt                      >  2:116358/1‑149 (MQ=255)
                                                                                                                             agcagcGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTgcg                  >  7:60609/1‑149 (MQ=255)
                                                                                                                                agcGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTTGGCCTg                                >  8:324662/1‑132 (MQ=255)
                                                                                                                                agcGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTTGGCCTg                                <  7:324662/132‑1 (MQ=255)
                                                                                                                                   ggTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGtt            <  7:99989/149‑1 (MQ=255)
                                                                                                                                      ttCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGcc         >  4:45878/1‑149 (MQ=255)
                                                                                                                                           gCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTtggt    <  1:194573/149‑1 (MQ=255)
                                                                                                                                            cAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACTGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTtggtg   >  1:241887/1‑149 (MQ=255)
                                                                                                                                             aGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTtggtgg  <  8:60609/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GCCGTCATGAGCGAAAAACGCACGCCCCCGCGCGAAGCCGCCCTGACCTTCCTGCGCGATAATCGCGAGGTGTGGAAAGCCTGGCTGCCGGCGGACATCGCCGCCAAGGTCGAGGCCAGCCTGTGAGCAGCGGTTTCCCGCAGACCCTGCAGTTTTCCTTTGCCGACAGCATCAACCGCCTGGTCGACTGGTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGG  >  NC_002947/4038080‑4038369

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: ‑