Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,278,468 G→C 6.5% intergenic (+208/‑22) 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,4680GC 6.5% 244.3 / 5.7 79intergenic (+208/‑22)PP_4652/purPmembrane protein/adenine:H+ symporter
Reads supporting (aligned to +/- strand):  ref base G (38/33);  new base C (2/3);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGC  >  NC_002947/5278330‑5278612
                                                                                                                                          |                                                                                                                                                
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                                                                                                                                     >  6:10482/1‑133 (MQ=255)
                 aCGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                     <  5:10482/135‑3 (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:102806/1‑149 (MQ=255)
                       aTAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                 <  6:79355/149‑1 (MQ=255)
                        tAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCAAGGCGCCCCCGGCGCCATgaa                                                                                                                                     <  2:251100/128‑3 (MQ=255)
                         aGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCt                                                                                                               >  5:91339/1‑149 (MQ=255)
                           gCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGc                                                                                                                         <  2:45782/137‑1 (MQ=255)
                           gCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGc                                                                                                                         >  1:45782/1‑137 (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                                                                                             <  2:187875/149‑1 (MQ=255)
                                           cTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGcaca                                                                                             <  1:7496/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                                                                                          >  6:36905/1‑145 (MQ=255)
                                                  gatgCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACacaa                                                                                          <  5:36905/145‑1 (MQ=255)
                                                       gcgGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGt                                                                                 <  6:119472/149‑1 (MQ=255)
                                                            ttttttCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCa                                                                            >  2:159898/1‑149 (MQ=255)
                                                                  cctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCgag                                                                       <  7:75079/148‑1 (MQ=255)
                                                                   ctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                     >  2:294030/1‑83 (MQ=255)
                                                                   ctccCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCCCCGGCGCCATgaa                                                                                                                                     <  1:294030/85‑3 (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                                                                                                                 <  2:68772/60‑1 (MQ=255)
                                                                                                                gTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAgg                                                                                                                 >  1:68772/1‑60 (MQ=255)
                                                                                                                   ggggTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATc                             >  8:91412/1‑141 (MQ=255)
                                                                                                                   ggggTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATc                             <  7:91412/141‑1 (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)
                                                                                                                                          |                                                                                                                                                
GAAAAACCCATTCCCCAACGGCCATAGGCCCCCGCCATGCATGCTGGTATGATGCGCGGCTTTTTTCCTCCCCACACAAAACCACGGCATCCGGTACGGTCTGTGCTTTGCTGTTGGGGTCGATACATTCATGGCGCCGGGGGCGCCTTGGGGAGCGGGCATGCTGGAAAGGCTGTTTCAACTAAGAGCACACAACACCAACGTGCGCACCGAGATTCTCGCGGGCGTCACCACCTTCCTGGCCATGGCCTACATCCTGTTCGTCAACCCGAGCATCCTCGGC  >  NC_002947/5278330‑5278612

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: