Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,424,882 G→C 5.1% H21D (CAC→GAC)  PP_5647 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,424,8820GC 5.1% 253.5 / 4.3 79H21D (CAC→GAC) PP_5647hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (40/35);  new base C (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

TAGTTTCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAA  >  NC_002947/4424735‑4425025
                                                                                                                                                     |                                                                                                                                               
tAGTTTCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGtga                                                                                                                                                <  7:163913/149‑1 (MQ=255)
   tttCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGtgatga                                                                                                                                             >  8:266598/1‑149 (MQ=255)
   tttCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGtgatga                                                                                                                                             >  8:138225/1‑149 (MQ=255)
           ggTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGt                                                                                                                                      <  7:215537/148‑1 (MQ=255)
              cTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCATCACGACGaaatccttg                                                                                                                                  <  3:127868/149‑8 (MQ=255)
                      gTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGTGCTTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTgg                                                                                                                            >  1:323092/1‑149 (MQ=255)
                      gTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGt                                                                                                                          >  2:215749/1‑149 (MQ=255)
                       tAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTc                                                                                                                         >  2:94042/1‑149 (MQ=255)
                           ccTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGcc                                                                                                                     >  4:160312/1‑149 (MQ=255)
                                 ctACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGAc                                                                                                                                            >  8:97819/1‑120 (MQ=255)
                                 ctACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGAc                                                                                                                                            <  7:97819/120‑1 (MQ=255)
                                   aCAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTc                                                                                                               >  3:392250/1‑147 (MQ=255)
                                   aCAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGACTGTACGCGAAAGACGGGTATGCACGTCTATACTACTGCCTTCCTCAAGGATTTCGTCGTGATGACGAAGGAGAAGCGATGGGTCTGCCCAGGTcc                                                                                                              >  5:181538/1‑148 (MQ=255)
                                    cAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTc                                                                                                               >  5:254881/1‑146 (MQ=255)
                                       ttCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAAt                                                                                                         <  1:94042/149‑1 (MQ=255)
                                            tCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTc                                                                                                    <  3:61290/149‑1 (MQ=255)
                                             ccccGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcg                                                                                                   >  4:197759/1‑149 (MQ=255)
                                              cccGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcgc                                                                                                  >  7:360890/1‑149 (MQ=255)
                                                 gCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTcgccgc                                                                                               <  4:392250/149‑1 (MQ=255)
                                                     ttGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                             <  7:138225/147‑1 (MQ=255)
                                                          gTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                      <  1:58656/149‑1 (MQ=255)
                                                          gTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                      <  1:215749/149‑1 (MQ=255)
                                                           tACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCa                                                                                      >  5:165073/1‑148 (MQ=255)
                                                             cacCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCatat                                                                                   <  7:266598/149‑1 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAgg                                                                                                                 <  7:235581/118‑1 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAgg                                                                                                                 >  8:235581/1‑118 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAg                                                                                                                  >  4:24884/1‑117 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAg                                                                                                                  <  3:24884/117‑1 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACAtctc                                                                                       >  1:325192/1‑144 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACAtctc                                                                                       <  2:325192/144‑1 (MQ=255)
                                                              acCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCatat                                                                                   >  7:91660/1‑148 (MQ=255)
                                                                cccGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGATGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCaa                                                                                >  5:152830/1‑149 (MQ=255)
                                                                 ccGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCaa                                                                                >  3:329737/1‑148 (MQ=255)
                                                                     tGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGCCATCTCATATCAATctcc                                                                           >  3:113124/1‑149 (MQ=255)
                                                                       aaGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcct                                                                          <  3:160312/148‑1 (MQ=255)
                                                                            cTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTcttact                                                                    >  3:393393/1‑149 (MQ=255)
                                                                                 t‑‑ttCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATcc                                                               >  3:270179/1‑149 (MQ=255)
                                                                                       gCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                            <  8:91660/148‑1 (MQ=255)
                                                                                       gCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                            <  4:270179/148‑1 (MQ=255)
                                                                                        ccTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAg                                                            >  1:6104/1‑147 (MQ=255)
                                                                                         cTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAg                                                         >  5:178872/1‑149 (MQ=255)
                                                                                         cTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAg                                                         <  6:165073/149‑1 (MQ=255)
                                                                                            tCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                             <  7:289978/110‑1 (MQ=255)
                                                                                            tCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGa                                                                                             >  8:289978/1‑110 (MQ=255)
                                                                                              cgcgAAAGACGGGTTTGCACTTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATTAAGAAGGTGAAGCTCTGGGTCTGCACATGTCCAGAATTTTTCTCCTCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGtcgc                                                     <  6:181538/148‑1 (MQ=255)
                                                                                              cgcgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGcgcg                                                    >  4:392211/1‑149 (MQ=255)
                                                                                              cgcgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGcgcg                                                    <  4:329737/149‑1 (MQ=255)
                                                                                                cgAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGt                                                  <  6:178872/149‑1 (MQ=255)
                                                                                                   aaGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTaaa                                               <  3:197759/149‑1 (MQ=255)
                                                                                                   aaGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTaaa                                               <  2:323092/149‑1 (MQ=255)
                                                                                                      aCGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAg                                                                                                                                        >  8:241590/1‑57 (MQ=255)
                                                                                                      aCGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAg                                                                                                                                        <  7:241590/57‑1 (MQ=255)
                                                                                                        gggTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCACAGTAGGCGCGGTAAagagag                                          >  8:110809/1‑149 (MQ=255)
                                                                                                         ggTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAgg                                         <  1:342793/149‑1 (MQ=255)
                                                                                                           tttGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGATGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGga                                       <  6:152830/149‑1 (MQ=255)
                                                                                                                aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgatagat                                  >  1:216737/1‑149 (MQ=255)
                                                                                                                aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgatagat                                  <  3:392211/149‑1 (MQ=255)
                                                                                                                aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgataga                                   >  7:311941/1‑148 (MQ=255)
                                                                                                                aCGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGgataga                                   >  6:270149/1‑148 (MQ=255)
                                                                                                                  gTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGatat                                >  7:135429/1‑149 (MQ=255)
                                                                                                                      aTACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCAt                            <  4:248339/149‑1 (MQ=255)
                                                                                                                          tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcctc                                                                         <  3:218416/100‑1 (MQ=255)
                                                                                                                          tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATctcctc                                                                         >  4:218416/1‑100 (MQ=255)
                                                                                                                          tCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                        >  5:381354/1‑149 (MQ=255)
                                                                                                                                cGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGGCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                  <  5:82974/149‑1 (MQ=255)
                                                                                                                                 gTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCTCTTAAAt                 <  5:270149/149‑1 (MQ=255)
                                                                                                                                     tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAagagag                                          <  7:361536/120‑1 (MQ=255)
                                                                                                                                     tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAagagag                                          >  8:361536/1‑120 (MQ=255)
                                                                                                                                     tCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGa             >  5:257460/1‑149 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATc                                                                <  8:183106/96‑1 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATc                                                                >  7:183106/1‑96 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                  >  5:323862/1‑142 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTaaa                  <  6:323862/142‑1 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGAcc           >  5:186363/1‑149 (MQ=255)
                                                                                                                                       aaGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGAcc           >  2:107026/1‑149 (MQ=255)
                                                                                                                                          gATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGt        <  2:6104/149‑1 (MQ=255)
                                                                                                                                          accttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTctta                                                                      <  4:334261/84‑1 (MQ=255)
                                                                                                                                          accttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTctta                                                                      >  3:334261/4‑87 (MQ=255)
                                                                                                                                           ccttCGTCATCACGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTg       >  3:376498/3‑149 (MQ=255)
                                                                                                                                           aTTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTg       <  6:254881/149‑1 (MQ=255)
                                                                                                                                               cgtcgtGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGaa    >  5:243821/1‑148 (MQ=255)
                                                                                                                                                gtcgtGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGaaaa  >  3:326193/1‑149 (MQ=255)
                                                                                                                                                    tgatgaCGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                        >  2:264577/1‑123 (MQ=255)
                                                                                                                                                    tgatgaCGAAGGTGAAGAGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCg                        <  1:264577/123‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
TAGTTTCGCTTGGTCTGGTGCTGTAAGCCTCTTCTACAATTCGGTCCCCGCTTTTGGAGTACACCCCGATGAAGAGCTGCGT‑‑TTCGCCTGTCCGCGAAAGACGGGTTTGCACGTCGATACTCCTGCCGTCCTCAAGGATTTCGTCGTGATGACGAAGGTGAAGCGCTGGGTCTGCCCAGGTCCAGAATTTTTCGCCGCGACATCTCATATCAATCTCCTCTTACTTATCCCAGTAGGCGCGGTAAAGAGAGGGATAGATATCCATTTCGCTTAAATGCGACCGGTGGAAAA  >  NC_002947/4424735‑4425025

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: