Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,185,785 C→G 5.6% R608R (CGG→CGC PP_1936 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,185,7850CG 5.6% 213.0 / 3.1 71R608R (CGG→CGCPP_1936membrane protein
Reads supporting (aligned to +/- strand):  ref base C (37/30);  new base G (2/2);  total (39/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 = 4.11e-01

TGTGTCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGTTTTAGGTAA  >  NC_002947/2185643‑2185932
                                                                                                                                              |                                                                                                                                                    
tgtgTCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGAGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCa                                                                                                                                                >  7:123795/1‑149 (MQ=255)
  tgtCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATg                                                                                                                                              <  1:111831/149‑1 (MQ=255)
    tCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCt                                                                                                                                            >  5:73525/1‑149 (MQ=255)
    tCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAACGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCGGGCTCATgat                                                                                                                                            <  5:19356/149‑3 (MQ=255)
      gCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTg                                                                                                                                          <  2:62014/149‑1 (MQ=255)
       cAGGCTGCCATGGAATCGTCCGATAGTTGCATATCGGTGGAGGTATTCAATGAAGGTGTTGATAGCATCTTCAGGGTCAATATTGTTGATGCTGGAAGTGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCttta                                                                                                                                         <  8:123795/149‑3 (MQ=255)
       cAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGATGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGa                                                                                                                                         <  7:286324/149‑1 (MQ=255)
       cAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGa                                                                                                                                         >  5:272106/1‑149 (MQ=255)
        aGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGAt                                                                                                                                        <  3:61362/149‑1 (MQ=255)
        aGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGAt                                                                                                                                        <  1:28679/149‑1 (MQ=255)
        aGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGAt                                                                                                                                        <  7:72386/149‑1 (MQ=255)
                 tCGATTCGTCAGATAGATGCATATCGGGGGATTTATTGAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGaa                                                                                                                               <  6:272106/149‑1 (MQ=255)
                 tCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTTATGCTTGATAAGATTGaa                                                                                                                               <  5:175163/149‑1 (MQ=255)
                 tCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTCGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGaa                                                                                                                               <  4:305316/149‑1 (MQ=255)
                           cGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGtttatctc                                                                                                                     >  1:325190/1‑143 (MQ=255)
                               aGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGtt                                                                                                                 >  1:247038/1‑148 (MQ=255)
                                   gCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAAt                                                                                                            >  5:140367/1‑149 (MQ=255)
                                           ggggAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGa                                                                                                    >  2:237607/1‑149 (MQ=255)
                                            gggATGTTTTAAATCAATGTGTTTATATCATATTTAGCGGCAATATTGTTGAAGCTGGAATATCTTTTGACGATCATCCTCCAAAATTGATCATGATCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGc                                                                                                      <  8:232783/146‑1 (MQ=255)
                                                ggTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTGTTTTGATGTTAAAATGTTAGCGAttttt                                                                                               >  8:106899/1‑149 (MQ=255)
                                                  tATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGAtttttt                                                                                              >  7:143697/1‑148 (MQ=255)
                                                    ttCAATCAAGGTGTTGATAGCATCTTAAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGATTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCAc                                                                                           >  5:363276/1‑149 (MQ=255)
                                                       aatcaaGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCAt                                                                                        >  1:305345/1‑149 (MQ=255)
                                                       aatcaaGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGATGGAAGATCTGTTGCCGAGCATCCTCCAAAATTGATAATGAGCCCGCTCATGATTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACAAt                                                                                        >  5:80509/1‑149 (MQ=255)
                                                         tcaaGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCAtt                                                                                       >  2:61161/1‑148 (MQ=255)
                                                          caaGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATtca                                                                                     >  3:150955/1‑149 (MQ=255)
                                                             ggTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATtcata                                                                                   >  5:330020/1‑147 (MQ=255)
                                                              gtgtTGATAGCATCATCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGAAGTTAAAATGGTAGCGATTTTTTCAACATTCAacac                                                                                 >  8:167054/1‑149 (MQ=255)
                                                                gatGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACt                                                                                <  1:81008/146‑1 (MQ=255)
                                                                 ttGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGa                                                                              <  2:219500/149‑1 (MQ=255)
                                                                 ttGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGa                                                                              <  6:140367/149‑1 (MQ=255)
                                                                     tAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAACAGCTGTAGAAGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGtt                                                                           >  5:82818/1‑148 (MQ=255)
                                                                                gCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTg                                                               <  7:167054/149‑1 (MQ=255)
                                                                                          ggTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAgcgc                                                      >  4:42760/1‑148 (MQ=255)
                                                                                              gAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGc                                                                     >  6:184189/1‑129 (MQ=255)
                                                                                              gAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGc                                                                     <  5:184189/129‑1 (MQ=255)
                                                                                              gAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGCAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAg                                                  <  6:73525/148‑1 (MQ=255)
                                                                                                aGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTaa                                                                         <  2:356424/123‑1 (MQ=255)
                                                                                                  cTGGAAGAGCTGTTGACGACCAGCCTACAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGTTATc                                             <  1:271463/149‑1 (MQ=255)
                                                                                                   tGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATcc                                            >  2:195793/1‑149 (MQ=255)
                                                                                                    ggAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCt                                           >  4:353692/1‑149 (MQ=255)
                                                                                                     gAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCt                                           <  1:237607/148‑1 (MQ=255)
                                                                                                      aaGATCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGa                                         <  5:366089/149‑1 (MQ=255)
                                                                                                       agagCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAg                                        <  2:298214/149‑1 (MQ=255)
                                                                                                             gTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGa                                         >  6:75656/1‑142 (MQ=255)
                                                                                                             gTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGa                                         <  5:75656/142‑1 (MQ=255)
                                                                                                             gTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCg                                  >  8:294515/1‑149 (MQ=255)
                                                                                                             gTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCg                                  >  8:143270/1‑149 (MQ=255)
                                                                                                                gACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGATTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCgggg                               >  5:101795/1‑149 (MQ=255)
                                                                                                                 aCGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGt                              <  8:143697/149‑1 (MQ=255)
                                                                                                                   gAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTTGCGATTTTTTCACCATTGATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGATGCTCCGGGTTgg                            >  8:197477/1‑149 (MQ=255)
                                                                                                                   gAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTgg                            >  5:230841/1‑149 (MQ=255)
                                                                                                                   gAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTgg                            >  1:290355/1‑149 (MQ=255)
                                                                                                                   gAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTgg                            >  3:130906/1‑149 (MQ=255)
                                                                                                                      cagcCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCTCGAAGGTATCCTGAGGCTCCGGGGTGGCa                          >  4:349059/1‑148 (MQ=255)
                                                                                                                        gcCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAAt                       <  2:247038/149‑1 (MQ=255)
                                                                                                                           tCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTATTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCt                     >  7:352277/1‑148 (MQ=255)
                                                                                                                             cAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTaa                   >  1:5129/1‑148 (MQ=255)
                                                                                                                             cAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTaa                   <  6:363276/148‑1 (MQ=255)
                                                                                                                              aaaaTTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAAc                  >  1:364389/1‑148 (MQ=255)
                                                                                                                                 aTTGATCATGATCCCGCTCATGCTTGATAAGATTTAAGTTTTTTGATGTTAAAATGGTAGAGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATc              <  5:58589/149‑1 (MQ=255)
                                                                                                                                  ttGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCa             <  2:305345/149‑1 (MQ=255)
                                                                                                                                    aagcATGAGCGGGCTCATGCTTGATAAGATTGAAGTTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTaa                   >  1:66917/4‑142 (MQ=255)
                                                                                                                                    aagcATGAGCGGGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTaa                   <  2:66917/138‑1 (MQ=255)
                                                                                                                                    aagcATGAGCGGGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGt           >  4:142639/4‑149 (MQ=255)
                                                                                                                                     aTCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTACACCATTCATCACTGAGATAAAAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGATACGGGGTGGCAAATCTAACAATCAGtt          >  7:104029/1‑149 (MQ=255)
                                                                                                                                      tCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGttt         <  2:5129/149‑1 (MQ=255)
                                                                                                                                      tCATGAGCCCGCTCATGCTTGATAAGATTCAAGTTTTTTCATGTTAAAATGGTATCGATTTTTTCACCATTCATCACTGAGTTAACATCTTCTAGGAACACCGCGAAGGTATCCTGAGTCTACGGGGTGGCAAATCTAACAATAAGatt         >  2:351491/1‑149 (MQ=255)
                                                                                                                                        aTGAGCCCGCTCATGCATGATAAGATTGAACTTTTTTGATGTTAAAATGGTAGCGAATTTTTCACAATTCATCACTCAGTTAACAGCGTCATGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACCATcctgtt         >  4:69536/1‑142 (MQ=255)
                                                                                                                                            gCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGTTTTAGGt    <  2:364389/148‑1 (MQ=255)
                                                                                                                                            gCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTAATCACTGAGTTAACAGCTTCTTGGTAATGCGCGAAGGTTTCATGTTGATCCGTGTTGGAAAATCTAACATTAAGTTTTAGGt    >  3:229709/1‑148 (MQ=255)
                                                                                                                                             cccGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGTTTTAGGTaa  <  1:61161/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGTGTCGCAGGCTGCCATCGAATCGTCCGATAGATGCATATCGGGGGAGGTATTCAATCAAGGTGTTGATAGCATCTTCAGCGGCAATATGGTTGAAGCTGGAAGAGCTGTTGACGAGCAGCCTCCAAAATTGATCATGAGCCCGCTCATGCTTGATAAGATTGAAGTTTTTTGATGTTAAAATGGTAGCGATTTTTTCACCATTCATCACTGAGTTAACAGCTTCTTGGTACAGCGCGAAGGTATCCTGAGGCTCCGGGGTGGCAAATCTAACAATCAGTTTTAGGTAA  >  NC_002947/2185643‑2185932

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: