Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 354,548 T→A 5.6% P208P (CCA→CCT cbcW ← choline/betaine/carnitine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947354,5480TA 5.6% 217.1 / 2.9 71P208P (CCA→CCTcbcWcholine/betaine/carnitine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (32/35);  new base A (2/2);  total (34/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 = 9.73e-01

CCAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAACACCACGGTG  >  NC_002947/354400‑354695
                                                                                                                                                    |                                                                                                                                                   
ccAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCtt                                                                                                                                                     <  3:159618/149‑1 (MQ=255)
  aGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCATga                                                                                                                                                   <  2:257230/149‑2 (MQ=255)
    gCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCa                                                                                                                                                 >  6:158103/1‑149 (MQ=255)
      ccGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATc                                                                                                                                               <  3:78569/149‑1 (MQ=255)
       cGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATcg                                                                                                                                              >  3:148728/1‑149 (MQ=255)
       cGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATcg                                                                                                                                              <  4:330675/149‑1 (MQ=255)
          ttCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGCCAGCGACAGCATGATGCATTGAGTGACCACGGCGGCGATGCTTGGCATCGCGt                                                                                                                                           >  4:207730/1‑149 (MQ=255)
          ttCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATcgc                                                                                                                                             >  4:162480/1‑147 (MQ=255)
          ttCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATcgc                                                                                                                                             <  3:162480/147‑1 (MQ=255)
          ttCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGt                                                                                                                                           >  6:45219/1‑149 (MQ=255)
            cGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGt                                                                                                                                           <  8:230482/147‑1 (MQ=255)
                  ccAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGt                                                                                                                                   >  2:77854/1‑149 (MQ=255)
                    aGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGTCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTc                                                                                                                                 >  3:141267/1‑149 (MQ=255)
                           tatCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACgg                                                                                                                          <  1:77854/149‑1 (MQ=255)
                           tatCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACgg                                                                                                                          <  4:148728/149‑1 (MQ=255)
                            atCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACggg                                                                                                                         >  3:246826/1‑149 (MQ=255)
                              cggcggTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACAAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGt                                                                                                                        >  1:226074/1‑148 (MQ=255)
                               ggcggTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTa                                                                                                                       <  3:322206/148‑1 (MQ=255)
                                     gtTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCt                                                                                                                <  7:137322/149‑1 (MQ=255)
                                       tCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCAAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTga                                                                                                              >  5:161681/1‑149 (MQ=255)
                                        cAGTGCGTTGTCCACAGGTTTGCCCAGGCCTTCAGGGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATGGCGTGCGGCAGTTCGATACGGGTAAGCAGCTgt                                                                                                              <  6:178255/148‑2 (MQ=255)
                                        cAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCTACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTga                                                                                                              <  7:54365/148‑1 (MQ=255)
                                          gTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCTGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTgacga                                                                                                           <  6:161681/149‑1 (MQ=255)
                                           tGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTgacgac                                                                                                          >  1:197967/1‑149 (MQ=255)
                                                  aCCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGTTAAGCAGCTGACGACGCGatca                                                                                                   <  6:25857/149‑1 (MQ=255)
                                                  aCCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGag                                                                                                     >  5:366892/1‑147 (MQ=255)
                                                     acaGGTTTGCCCAGGTCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGcc                                                                                                <  4:141267/149‑1 (MQ=255)
                                                           ttGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCaccctc                                                                                          >  8:101305/1‑144 (MQ=255)
                                                              cccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCa                                                                                                                                     <  1:44080/103‑1 (MQ=255)
                                                              cccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCa                                                                                                                                     >  2:44080/1‑103 (MQ=255)
                                                               ccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGagcagc                                                                                                 >  6:244110/1‑138 (MQ=255)
                                                               ccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGagcagc                                                                                                 <  5:244110/138‑1 (MQ=255)
                                                               ccAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTg                                                                                      >  3:82795/1‑149 (MQ=255)
                                                                 aGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGcc                                                                                    <  2:197967/149‑1 (MQ=255)
                                                                    ccATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGc                                                                                 <  8:168972/149‑1 (MQ=255)
                                                                      aTCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGt                                                                               <  5:166696/149‑1 (MQ=255)
                                                                      aTCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGt                                                                               <  5:158103/149‑1 (MQ=255)
                                                                      aTCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCATGGCATCGCGTGCGGCAGTTCGATACGGGTAAGAAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGt                                                                               >  1:276550/1‑149 (MQ=255)
                                                                       tCAGCGCCAACCAGGGCGGCAATGACCCCCATCGACAGCGACAGCTTGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTc                                                                              >  1:252871/1‑149 (MQ=255)
                                                                              cAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAAc                                                                       <  6:366892/149‑1 (MQ=255)
                                                                              cAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAAc                                                                       <  5:45219/149‑1 (MQ=255)
                                                                                 ccAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCaaa                                                                                             <  2:48810/124‑1 (MQ=255)
                                                                                 ccAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCaaa                                                                                             >  1:48810/1‑124 (MQ=255)
                                                                                  cAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCtt                                                                   >  1:195453/1‑149 (MQ=255)
                                                                                   aGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTg                                                                  >  4:321976/1‑149 (MQ=255)
                                                                                    gggCGGCAATGACCAGCATCGACAGCGACAGCATGCTGCTTTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGc                                                                 <  3:207730/149‑1 (MQ=255)
                                                                                               accaccATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCg                                                        >  6:332197/1‑147 (MQ=255)
                                                                                                 caccaTCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCg                                                        >  3:324401/1‑145 (MQ=255)
                                                                                                  accaTCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGAt                                                   <  3:321976/149‑1 (MQ=255)
                                                                                                            gcgacagcATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCaggtag                                         <  2:195453/149‑1 (MQ=255)
                                                                                                              gacagcATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCaggtaggt                                       >  4:321030/1‑149 (MQ=255)
                                                                                                                   cATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                                  >  7:205797/1‑149 (MQ=255)
                                                                                                                   cATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                                  <  1:97133/149‑1 (MQ=255)
                                                                                                                      gatgCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGc                                                                                                                                         <  1:230951/43‑1 (MQ=255)
                                                                                                                      gatgCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGc                                                                                                                                         >  2:230951/1‑43 (MQ=255)
                                                                                                                            ttGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGgcg                         >  5:155460/1‑149 (MQ=255)
                                                                                                                               aGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCgcag                      >  6:143546/1‑149 (MQ=255)
                                                                                                                               aGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCgca                       <  2:226074/148‑1 (MQ=255)
                                                                                                                                   acacCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAAt                  <  3:323926/149‑1 (MQ=255)
                                                                                                                                       cggcggCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGc               >  5:180388/1‑148 (MQ=255)
                                                                                                                                        ggcggcGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCgcag                      <  1:248561/140‑1 (MQ=255)
                                                                                                                                        ggcggcGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCgaag                      >  2:248561/1‑140 (MQ=255)
                                                                                                                                        gcacgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGc                                                       >  3:117341/5‑107 (MQ=255)
                                                                                                                                        gcacgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGc                                                       <  4:117341/103‑1 (MQ=255)
                                                                                                                                         gcggcgATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGaa            <  6:180388/149‑1 (MQ=255)
                                                                                                                                         gcggcgATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGaa            >  7:215396/1‑149 (MQ=255)
                                                                                                                                          cggcgATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAac           <  4:324401/149‑1 (MQ=255)
                                                                                                                                          acgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGc                                                                                 >  8:64243/3‑79 (MQ=255)
                                                                                                                                          acgcgATGCCAAGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGc                                                                                 <  7:64243/77‑1 (MQ=255)
                                                                                                                                           ggcgATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATCGGCGCAGCAATGGCGAaca          <  8:215396/149‑1 (MQ=255)
                                                                                                                                                 gCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                                  >  4:125660/1‑119 (MQ=255)
                                                                                                                                                 gCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGc                                  <  3:125660/119‑1 (MQ=255)
                                                                                                                                                   ttGGCATCGCGTGCGGCAGTTCGGTACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAACACCACGGTg  >  5:300303/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCAGGCCCGCTTCGAAGCCCAGGGAGATATCGGCGGTGTTCAGTGCGTTGACCACAGGTTTGCCCAGGCCATCAGCGCCAACCAGGGCGGCAATGACCACCATCGACAGCGACAGCATGATGCATTGAGTGACACCGGCGGCGATGCTTGGCATCGCGTGCGGCAGTTCGATACGGGTAAGCAGCTGACGACGCGAGCAGCCAAAGGCCTTGCCGGCGTCCATCAACTCTTGCGGCACGTCGCAGATGCCCAGGTAGGTCAGGCGAATGGGCGCAGCAATGGCGAACACCACGGTG  >  NC_002947/354400‑354695

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: