Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,360,770 T→C 11.6% I118V (ATC→GTC) ‡ PP_3835 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,360,7700TC11.6% 171.4 / 12.9 69I118V (ATC→GTC) ‡PP_3835two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (31/30);  new base C (5/3);  total (36/33)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

GTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCGCA  >  NC_002947/4360624‑4360912
                                                                                                                                                  |                                                                                                                                               
gttgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCt                                                                                                                                               >  5:88918/1‑149 (MQ=255)
gttgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCt                                                                                                                                               >  1:25816/1‑149 (MQ=255)
gttgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCt                                                                                                                                               >  3:329028/1‑149 (MQ=255)
 ttgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCGGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTg                                                                                                                                              <  6:84549/149‑1 (MQ=255)
  tgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTgg                                                                                                                                             >  5:401519/1‑149 (MQ=255)
  tgAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTgg                                                                                                                                             >  4:244017/1‑149 (MQ=255)
      gTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCAc                                                                                                                                         <  3:351712/149‑1 (MQ=255)
       tCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaag                                                                                                                                        <  1:100526/149‑3 (MQ=255)
         gTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaagg                                                                                                                                       >  5:369567/1‑145 (MQ=255)
         gTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGATCATTGGCCTGGCCaagg                                                                                                                                       <  6:369567/148‑4 (MQ=255)
          tGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtc                                                                                                                                     <  3:321561/149‑1 (MQ=255)
          tGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtc                                                                                                                                     <  4:329028/149‑1 (MQ=255)
            cgcCAGACCTCCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtca                                                                                                                                    <  3:26105/148‑1 (MQ=255)
            cgcCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCGGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtcat                                                                                                                                   <  2:166043/149‑1 (MQ=255)
              ccAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTtcatca                                                                                                                                 <  1:268268/149‑1 (MQ=255)
                  accaccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCAACACGGa                                                                                                                             >  5:195477/1‑149 (MQ=255)
                   ccaccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGAt                                                                                                                            <  7:337141/149‑1 (MQ=255)
                    caccGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGAt                                                                                                                            >  4:90385/1‑148 (MQ=255)
                            aaCTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAg                                                                                                                   >  6:321961/1‑149 (MQ=255)
                                 gCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTgcc                                                                                                               <  6:418028/148‑1 (MQ=255)
                                          ttCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGccc                                                                                                     <  6:355313/149‑1 (MQ=255)
                                            cAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCa                                                                                                    <  5:7820/148‑1 (MQ=255)
                                               gCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAAc                                                                                                >  1:284870/1‑149 (MQ=255)
                                                 ttCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTg                                                                                              >  1:190546/1‑149 (MQ=255)
                                                     gcgcCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTg                                                                                         >  6:240431/1‑149 (MQ=255)
                                                        ccTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCg                                                                                      >  5:437107/1‑149 (MQ=255)
                                                         cTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGt                                                                                     >  8:219281/1‑149 (MQ=255)
                                                          tGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTc                                                                                    >  6:9461/1‑149 (MQ=255)
                                                              aCAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTg                                                                                <  6:195477/149‑1 (MQ=255)
                                                                 gCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGAGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTgcg                                                                              >  2:376115/1‑148 (MQ=255)
                                                                       ctctcACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGc                                                                       <  6:404288/145‑1 (MQ=255)
                                                                       cgcgcACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGc                                                                       <  1:42066/149‑1 (MQ=255)
                                                                       cgcgcACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGc                                                                       <  5:321961/149‑1 (MQ=255)
                                                                                cGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCg                                                              >  2:349078/1‑149 (MQ=255)
                                                                                  tCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTc                                                            <  5:240431/149‑1 (MQ=255)
                                                                                    gAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGAGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTcc                                                           <  1:376115/148‑1 (MQ=255)
                                                                                     aGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTg                                                         <  3:377670/149‑1 (MQ=255)
                                                                                           tgcGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCcac                                                   <  4:390308/149‑1 (MQ=255)
                                                                                             cgcCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCCCTTTCATCACGGATCACATCCAGTGCCGCCACGGCCCAGCCCTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCccac                                                 >  5:350835/1‑149 (MQ=255)
                                                                                             cgcCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGTGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCcac                                                   <  1:303651/147‑1 (MQ=255)
                                                                                             cgcCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgcgc                                                      >  1:419709/1‑144 (MQ=255)
                                                                                             cgcCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTgcgc                                                      <  2:419709/144‑1 (MQ=255)
                                                                                                     cgcgACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTc                                          >  7:381915/1‑148 (MQ=255)
                                                                                                       cgACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATTATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAg                                       <  6:401519/149‑1 (MQ=255)
                                                                                                        gACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAgc                                      >  2:118571/1‑149 (MQ=255)
                                                                                                           gctCAGTGATATCACGTGTCACCTAGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCAGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGcc                                   >  3:221553/3‑149 (MQ=255)
                                                                                                               aGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTgg                                <  5:9461/148‑1 (MQ=255)
                                                                                                               aGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGc                               >  8:358144/1‑149 (MQ=255)
                                                                                                               aGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGc                               >  4:284328/1‑149 (MQ=255)
                                                                                                                    aatCACGTGTCACCTTGGCCAGGCCAATGATCTTGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCAcag                           <  1:349078/148‑1 (MQ=255)
                                                                                                                      tCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCAcag                           >  7:385502/1‑146 (MQ=255)
                                                                                                                        aCGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCa                                        <  8:152829/131‑1 (MQ=255)
                                                                                                                        aCGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCa                                        >  7:152829/1‑131 (MQ=255)
                                                                                                                         cGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACg                     <  1:118571/149‑1 (MQ=255)
                                                                                                                          gtgtCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGc                    >  8:302384/1‑149 (MQ=255)
                                                                                                                            atgaaagtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCt                   >  4:266378/8‑148 (MQ=255)
                                                                                                                               aCCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCTTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAgg               >  4:331212/1‑149 (MQ=255)
                                                                                                                                ccTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAggg              >  5:360747/1‑149 (MQ=255)
                                                                                                                                aagtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAggg              >  2:296166/4‑149 (MQ=255)
                                                                                                                                 agtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGggc             >  5:16170/3‑149 (MQ=255)
                                                                                                                                 agtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGggc             >  3:95915/3‑149 (MQ=255)
                                                                                                                                  gtGGCCAGATCATTGGCCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGggcg            <  3:266378/148‑1 (MQ=255)
                                                                                                                                    tctcaggccaATGATCTGGCCACTTTCATCACGGATCAAATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCAGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACATCAGACGCTCAGggcgg           <  2:181832/145‑1 (MQ=255)
                                                                                                                                    ggccaggccaATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGggcgg           <  3:244017/148‑1 (MQ=255)
                                                                                                                                        aggccaATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATc      <  2:284870/149‑1 (MQ=255)
                                                                                                                                        aggccaATGATCTGGCCACTTTCATCAAGGATCACATCCAGTGACGCCAAGGACCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGggcgg           >  1:181832/1‑144 (MQ=255)
                                                                                                                                         ggccaATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATcc     >  1:174180/1‑149 (MQ=255)
                                                                                                                                           ccaATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCg    >  8:10555/1‑148 (MQ=255)
                                                                                                                                            caATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGGCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCGca  <  4:221553/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GTTGAAGTCGTGCGCCAGACCACCGGTCAACTGGCCGAGGGCTTCAAGCTTCTGCGCCTGGAACAGCTGGGCGCGCACTGCGTCGAGCTGCTGCGCCGTCTCGCGACGGTCAGTGATATCACGTGTCACCTTGGCCAGGCCAATGATCTGGCCACTTTCATCACGGATCACATCCAGTGCCGCCAAGGCCCAGAACTGCGTGCCGTCCTTGCGCACCCGCCAGCCTTCGTCCTGCGCCACACCTTGCTCCAGCGCCTGGCACAGCAGACGCTCAGGGCGGCCATCCGCA  >  NC_002947/4360624‑4360912

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: