Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,989,556 T→C 20.2% C89C (TGT→TGC PP_5247 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,989,5560TC20.2% 136.1 / 28.7 70C89C (TGT→TGCPP_5247hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/28);  new base C (7/7);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TCCCCGATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTGTATCGGCC  >  NC_002947/5989408‑5989704
                                                                                                                                                     |                                                                                                                                                    
tCCCCGATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACtgt                                                                                                                                                      <  7:141798/149‑1 (MQ=255)
  cccGATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTTGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTgg                                                                                                                                                    >  3:265201/1‑149 (MQ=255)
  cccGATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTgg                                                                                                                                                    >  8:143803/1‑149 (MQ=255)
     gATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAAGCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTa                                                                                                                                                 <  8:88539/149‑1 (MQ=255)
          aCCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAg                                                                                                                                            >  3:248697/1‑149 (MQ=255)
            cccTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGcc                                                                                                                                          >  6:71058/1‑149 (MQ=255)
             ccTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAAAACCTGCCCGGTGCCTACGACTGCGGCTACAAAGCCt                                                                                                                                         >  6:318303/1‑149 (MQ=255)
                 cccAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGa                                                                                                                                     <  2:354143/149‑1 (MQ=255)
                  ccAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATTTCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGa                                                                                                                                     <  1:103585/149‑1 (MQ=255)
                  ccAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTg                                                                                                                                       <  1:45890/146‑1 (MQ=255)
                  ccAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTg                                                                                                                                       >  2:45890/1‑146 (MQ=255)
                   cAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGa                                                                                                                                     >  5:218558/1‑147 (MQ=255)
                   cAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGa                                                                                                                                     <  6:218558/147‑1 (MQ=255)
                     acacCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCAt                                                                                                                                 <  7:143803/149‑1 (MQ=255)
                           cacaACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGAt                                                                                                                           <  4:248697/149‑1 (MQ=255)
                            acaACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGAt                                                                                                                           >  2:329327/1‑148 (MQ=255)
                             caACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTTGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGa                                                                                                                         <  4:265201/149‑1 (MQ=255)
                             caACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGa                                                                                                                         <  5:284426/149‑1 (MQ=255)
                                 gCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTg                                                                                                                                       >  1:311124/1‑131 (MQ=255)
                                 gCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTg                                                                                                                                       <  4:280461/131‑1 (MQ=255)
                                 gCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTg                                                                                                                                       <  2:311124/131‑1 (MQ=255)
                                 gCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTg                                                                                                                                       >  3:280461/1‑131 (MQ=255)
                                 gCTCTTTCAGCCGCTGCGATGCCGGCAATCCCCCATTATTCACCTTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCTCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGAcaca                                                                                                                     <  4:260860/149‑1 (MQ=255)
                                         aGCCGCGGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCCGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGGGGCTACAAAGCCTTGCAGCATCTGGATGACACAGGctctt                                                                                                              <  7:52711/148‑5 (MQ=255)
                                          gCCGCTTCGATGCCCGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGc                                                                                                            <  4:121025/149‑1 (MQ=255)
                                                         gCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAAc                                                                                             <  5:72960/149‑1 (MQ=255)
                                                         gCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAAc                                                                                             >  6:241317/1‑149 (MQ=255)
                                                                    ttattCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACt                                                                                  >  6:17826/1‑149 (MQ=255)
                                                                      attCACCGTGAACCCCGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTc                                                                                 <  6:25982/148‑1 (MQ=255)
                                                                       ttCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCg                                                                                >  7:262642/1‑148 (MQ=255)
                                                                       ttCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCg                                                                                >  4:279999/1‑148 (MQ=255)
                                                                          aCCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGaaa                                                                             <  5:241317/148‑1 (MQ=255)
                                                                                  cccAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTAGGAGCATCTGGATGACACAGGCAAGCGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAg                                                                    >  7:234470/1‑149 (MQ=255)
                                                                                   ccAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAAc                                                                                             <  1:49419/123‑1 (MQ=255)
                                                                                   ccAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAAc                                                                                             >  2:49419/1‑123 (MQ=255)
                                                                                     aGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCTGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggt                                                                 <  1:329327/149‑1 (MQ=255)
                                                                                          tctcGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggttggtt                                                            >  5:117160/1‑149 (MQ=255)
                                                                                           ctcGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAgg                                                                   >  4:279058/1‑141 (MQ=255)
                                                                                           ctcGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAgg                                                                   <  3:279058/141‑1 (MQ=255)
                                                                                              gACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGa                                                          <  5:17826/147‑1 (MQ=255)
                                                                                              gACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAgg                                                        >  3:63313/1‑149 (MQ=255)
                                                                                                cccACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAgg                                                        <  1:159343/147‑1 (MQ=255)
                                                                                                 ccACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCAt                                                     >  6:51314/1‑149 (MQ=255)
                                                                                                  cACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCAtt                                                    >  2:64174/1‑149 (MQ=255)
                                                                                                     gATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCt                                                 <  5:51314/149‑1 (MQ=255)
                                                                                                      aTGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTg                                                >  7:24838/1‑149 (MQ=255)
                                                                                                           ccTTGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAATCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGAt                                             <  5:71058/148‑1 (MQ=255)
                                                                                                           ccTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATg                                            >  8:341320/1‑149 (MQ=255)
                                                                                                               gTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGa                                         <  6:346935/148‑1 (MQ=255)
                                                                                                                tccATGTCGCCCAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGAtt                                       <  5:318303/147‑1 (MQ=255)
                                                                                                                tGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGAtt                                       <  4:63313/149‑1 (MQ=255)
                                                                                                                  cATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGa                                     <  8:24838/149‑1 (MQ=255)
                                                                                                                    tGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGaca                                                                                                                       <  4:266350/65‑1 (MQ=255)
                                                                                                                    tGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGaca                                                                                                                       >  3:266350/1‑65 (MQ=255)
                                                                                                                    tGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGaaa                                                                             <  5:281345/107‑1 (MQ=255)
                                                                                                                    tGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGaaa                                                                             >  6:281345/1‑107 (MQ=255)
                                                                                                                         cccAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGa                                                                                                  >  2:334549/1‑81 (MQ=255)
                                                                                                                         cccAATACCTGCGCGGTGCCTACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGa                                                                                                  <  1:334549/81‑1 (MQ=255)
                                                                                                                             aTACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggtt                                                                >  7:240483/1‑111 (MQ=255)
                                                                                                                             aTACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggtt                                                                <  8:240483/111‑1 (MQ=255)
                                                                                                                                ccTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTTGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATg                        >  6:292927/1‑148 (MQ=255)
                                                                                                                                ccTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTTGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATg                        >  2:311995/1‑148 (MQ=255)
                                                                                                                                 cTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGaa                      <  3:279999/149‑1 (MQ=255)
                                                                                                                                    cgcgGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACt                    <  1:64174/148‑1 (MQ=255)
                                                                                                                                           ccTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGcc            <  6:117160/149‑1 (MQ=255)
                                                                                                                                            cTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATAGGATTGAGTCACAGCCCATGAACTGAGTCGCCt           >  6:167330/1‑149 (MQ=255)
                                                                                                                                             tACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggtt                                                                >  7:69969/1‑95 (MQ=255)
                                                                                                                                             tACGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGggtt                                                                <  8:69969/95‑1 (MQ=255)
                                                                                                                                              aCGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTg          >  7:103004/1‑148 (MQ=255)
                                                                                                                                               cGACTGCGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTGTa        >  1:45118/1‑149 (MQ=255)
                                                                                                                                                  cTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTGTATCg     >  1:71386/1‑149 (MQ=255)
                                                                                                                                                     tgGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTGTATCGGcc  >  8:27282/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCCCGATCCACCCCTCCCCAACACCACACAACGCTCTTTCAGCCGCTGCGATGCCGGCCATCCCCCATTATTCACCGTGAACCCAGGCGTCTCGACCCACGATGCCCTGGTGCATGTCGCCCAATACCTGCGCGGTGCCTACGACTGTGGCTACAAAGCCTTGGAGCATCTGGATGACACAGGCAAGAGCCTGTTCTGGAGCAACCTGAACGCACTCGAAATGGCAGAAGGGTTGGTTGAGGCATTGCTGGATGGGATTGAGTCACAGCCCATGAACTGAGTCGCCTGTATCGGCC  >  NC_002947/5989408‑5989704

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: