Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,278,469 G→C 7.5% intergenic (+209/‑21) PP_4652 → / → purP membrane protein/adenine:H+ symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,278,4690GC 7.5% 233.1 / 6.3 79intergenic (+209/‑21)PP_4652/purPmembrane protein/adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (40/33);  new base C (2/4);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGCGAGAC  >  NC_002947/5278330‑5278617
                                                                                                                                           |                                                                                                                                                    
gAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCtt                                                                                                                                             >  2:222019/1‑149 (MQ=255)
  aaaaCCCATTCCCCACCGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTgg                                                                                                                                           <  4:32110/149‑1 (MQ=255)
       ccATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTggg                                                                                                                                          <  5:139/145‑1 (MQ=255)
       ccATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTggg                                                                                                                                          >  6:139/1‑145 (MQ=255)
       ccATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGc                                                                                                                                      >  5:79355/1‑149 (MQ=255)
          ttCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGTGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCggg                                                                                                                                   <  7:122785/149‑1 (MQ=255)
          ttCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCggg                                                                                                                                   >  5:119472/1‑149 (MQ=255)
             cccAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCAt                                                                                                                                >  1:187875/1‑149 (MQ=255)
             cccAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCAt                                                                                                                                >  2:7496/1‑149 (MQ=255)
              ccAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCTGGGGCGAGTTGGGGAGCGGGCagg                                                                                                                               >  6:215460/1‑147 (MQ=255)
              ccAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATg                                                                                                                               >  1:65122/1‑149 (MQ=255)
               cAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCg                                                                                                                                     >  7:203037/1‑142 (MQ=255)
               cAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCg                                                                                                                                     <  8:203037/142‑1 (MQ=255)
                 aCGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                          <  5:10482/135‑3 (MQ=255)
                 aCGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                          >  6:10482/1‑133 (MQ=255)
                  cGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTgg                                                                                                                           <  1:222019/149‑1 (MQ=255)
                       aTAGGCCCCCTCCATTCCTGCTGGTATTATGCGCGGCTTTGTTCCTCCCCACACAAAACCCCGTCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      <  8:256862/149‑1 (MQ=255)
                       aTAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      <  5:102806/149‑1 (MQ=255)
                       aTAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      <  6:79355/149‑1 (MQ=255)
                       aTAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      >  6:102806/1‑149 (MQ=255)
                        tAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCAAGGCGCCCCCGGCGCCATgaa                                                                                                                                          <  2:251100/128‑3 (MQ=255)
                         aGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCt                                                                                                                    >  5:91339/1‑149 (MQ=255)
                           gCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGc                                                                                                                              >  1:45782/1‑137 (MQ=255)
                           gCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGc                                                                                                                              <  2:45782/137‑1 (MQ=255)
                               ccGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTc                                                                                                               >  3:154830/1‑148 (MQ=255)
                               ccGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTc                                                                                                               <  4:154830/148‑1 (MQ=255)
                                   catgcatgCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCAAGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACt                                                                                                           >  2:59571/1‑148 (MQ=255)
                                       catgcTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAagag                                                                                                      <  8:85583/149‑1 (MQ=255)
                                       catgcTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAagag                                                                                                      >  8:3642/1‑149 (MQ=255)
                                         tgcTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGca                                                                                                    <  7:3642/149‑1 (MQ=255)
                                           cTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGcaca                                                                                                  <  1:7496/149‑1 (MQ=255)
                                           cTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGcaca                                                                                                  <  2:187875/149‑1 (MQ=255)
                                              gTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACacaa                                                                                               >  3:132320/1‑149 (MQ=255)
                                              gTATGATGCGCGGCTTTTTTCCTCCCCACAAAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGcacacat                                                                                               >  5:80075/1‑148 (MQ=255)
                                               tATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACacaac                                                                                              <  5:268314/149‑1 (MQ=255)
                                                  gatgCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACacaa                                                                                               <  5:36905/145‑1 (MQ=255)
                                                  gatgCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACacaa                                                                                               >  6:36905/1‑145 (MQ=255)
                                                       gcgGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGt                                                                                      <  6:119472/149‑1 (MQ=255)
                                                            ttttttCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCa                                                                                 >  2:159898/1‑149 (MQ=255)
                                                                  cctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCgag                                                                            <  7:75079/148‑1 (MQ=255)
                                                                   ctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                          <  1:294030/85‑3 (MQ=255)
                                                                   ctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                          >  2:294030/1‑83 (MQ=255)
                                                                      cccACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATtct                                                                       >  1:166600/1‑149 (MQ=255)
                                                                       ccACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATtctc                                                                      >  7:286633/1‑149 (MQ=255)
                                                                        cacacaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGcacac                                                                                                 >  8:105537/1‑121 (MQ=255)
                                                                        cacacaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGc                                                                                                     <  8:187850/117‑1 (MQ=255)
                                                                        cacacaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGc                                                                                                     >  7:187850/1‑117 (MQ=255)
                                                                        cacacaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTcg                                                                     >  1:100932/1‑149 (MQ=255)
                                                                        cacacaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTcg                                                                     >  2:228604/1‑149 (MQ=255)
                                                                           acaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGc                                                                                                     <  8:69499/114‑1 (MQ=255)
                                                                           acaAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGc                                                                                                     >  7:69499/1‑114 (MQ=255)
                                                                               aaCCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGt                                                              <  2:65122/149‑1 (MQ=255)
                                                                                 ccACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTca                                                            <  5:215460/149‑1 (MQ=255)
                                                                                  cACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTcac                                                           >  6:196247/1‑149 (MQ=255)
                                                                                  cACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTcac                                                           >  3:60835/1‑148 (MQ=255)
                                                                                                  gTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATgg                                          <  4:60835/149‑1 (MQ=255)
                                                                                                       tgCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTa                                      <  2:166600/149‑1 (MQ=255)
                                                                                                       tgCTTTGCTGTTGGGGTCGATACATTCATGGCACCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCt                                       >  5:31559/1‑148 (MQ=255)
                                                                                                         cTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTAc                                     <  2:100932/148‑1 (MQ=255)
                                                                                                          tttGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACAt                                   >  6:200494/1‑149 (MQ=255)
                                                                                                            tGCTGTTGGTGTCCATACATTCATCGCGCCGCGGGCGCCTTGGGGAGCGGGCATGCTGCAAAGGCTGTATCAACTAAGAGCACACAACACCATCCTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATc                                  <  7:252194/148‑1 (MQ=255)
                                                                                                               tGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGt                              >  4:84566/1‑149 (MQ=255)
                                                                                                                gTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      >  1:68772/1‑60 (MQ=255)
                                                                                                                gTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                      <  2:68772/60‑1 (MQ=255)
                                                                                                                   ggggTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATc                                  <  7:91412/141‑1 (MQ=255)
                                                                                                                   ggggTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATc                                  >  8:91412/1‑141 (MQ=255)
                                                                                                                   ggggTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCGTCCTGGCCATGGCCTACATCCTGTTCGt                          >  5:180539/1‑149 (MQ=255)
                                                                                                                       tCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAAc                      >  5:124838/1‑149 (MQ=255)
                                                                                                                         gATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAAcc                     <  5:196247/148‑1 (MQ=255)
                                                                                                                           tACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGa                  <  1:228604/149‑1 (MQ=255)
                                                                                                                              aTTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACGCGa                  <  6:84547/146‑1 (MQ=255)
                                                                                                                              aTTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGa                  >  5:84547/1‑146 (MQ=255)
                                                                                                                                 cATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATcc            >  8:104164/1‑149 (MQ=255)
                                                                                                                                  aTGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTc                            <  6:268420/132‑1 (MQ=255)
                                                                                                                                  aTGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTc                            >  5:268420/1‑132 (MQ=255)
                                                                                                                                   tGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCgg                                                                  >  6:57430/1‑93 (MQ=255)
                                                                                                                                   tGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCgg                                                                  <  5:57430/93‑1 (MQ=255)
                                                                                                                                      cgcCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGc       <  5:200494/149‑1 (MQ=255)
                                                                                                                                       gcCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCCACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGc       >  7:251663/1‑148 (MQ=255)
                                                                                                                                           ggggCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGACCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGCGAGAc  <  7:104164/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGCGAGAC  >  NC_002947/5278330‑5278617

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: