Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,041,202 A→G 9.9% intergenic (+180/‑482) PP_5680 → / → PP_4447 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,041,2020AG 9.9% 215.9 / 10.0 80intergenic (+180/‑482)PP_5680/PP_4447Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/32);  new base G (6/2);  total (46/34)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTG  >  NC_002947/5041057‑5041347
                                                                                                                                                  |                                                                                                                                                  
cTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGcc                                                                                                                                                 >  5:51526/1‑149 (MQ=255)
  ccACACCCAGTCTCACATCTGACCCTGGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGccc                                                                                                                                                >  1:17665/1‑149 (MQ=255)
  ccACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccg                                                                                                                                               >  3:191619/1‑149 (MQ=255)
    acacCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                             >  6:8605/1‑149 (MQ=255)
      acCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                             >  2:231315/1‑147 (MQ=255)
       cccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCTCGGCGGGCTTAAGGCGCGCCgatg                                                                                                                                         <  2:332621/149‑3 (MQ=255)
        ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                          <  8:171640/148‑2 (MQ=255)
        ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                          >  7:171640/1‑147 (MQ=255)
        ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCgag                                                                                                                                         >  7:68216/1‑148 (MQ=255)
        ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGAGGCCAGCACGTCTGCGATGAGGACATGACCTCTTCAGCGGCGCGCCGTAGGGCCGCCgag                                                                                                                                         >  6:338546/1‑148 (MQ=255)
           gTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                     <  5:8605/149‑1 (MQ=255)
           gTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                     <  1:346161/149‑1 (MQ=255)
              tcACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                                   >  8:76869/1‑148 (MQ=255)
               cacaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                                   >  7:159637/1‑147 (MQ=255)
                acttCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                                   <  7:66040/143‑1 (MQ=255)
                acaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                                   >  1:9823/1‑146 (MQ=255)
                       aCCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGAcccc                                                                                                                         <  6:338770/149‑1 (MQ=255)
                       aCCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGAc                                                                                                                            <  8:159637/146‑1 (MQ=255)
                                   gACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGAccc                                                                                                              <  1:43672/149‑1 (MQ=255)
                                    aCGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGAcccc                                                                                                             >  6:185502/1‑149 (MQ=255)
                                      gCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCt                                                                                                            >  2:295320/1‑148 (MQ=255)
                                        tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                     <  8:101091/121‑1 (MQ=255)
                                        tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                     >  7:101091/1‑121 (MQ=255)
                                                   ttGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                               <  2:17665/148‑1 (MQ=255)
                                                     gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                             <  1:231315/148‑1 (MQ=255)
                                                         cacaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                        >  3:148854/1‑149 (MQ=255)
                                                          acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                        >  4:381026/1‑148 (MQ=255)
                                                          acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                        >  8:279582/1‑148 (MQ=255)
                                                          acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttagg                                                                                       <  6:51526/149‑4 (MQ=255)
                                                            agagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                                     <  5:14412/149‑1 (MQ=255)
                                                                ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                      <  4:57661/144‑1 (MQ=255)
                                                                ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                      >  3:57661/1‑144 (MQ=255)
                                                                ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCt                                                                                 <  8:8261/149‑1 (MQ=255)
                                                                     gggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                             <  6:15971/148‑1 (MQ=255)
                                                                      ggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGt                                                                           >  4:104252/1‑149 (MQ=255)
                                                                        aTAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTct                                                                         <  7:76869/149‑1 (MQ=255)
                                                                         tAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctg                                                                        <  8:295038/149‑1 (MQ=255)
                                                                            ccTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                             >  2:174416/1‑125 (MQ=255)
                                                                              ttAAGGCGTAATTGTCTGATGCCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCtgctgc                                                                    >  3:375534/1‑148 (MQ=255)
                                                                                 aGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                                <  2:9823/149‑1 (MQ=255)
                                                                                   gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                             <  8:307894/134‑1 (MQ=255)
                                                                                   gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                             >  7:307894/1‑134 (MQ=255)
                                                                                   gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                              >  2:388435/1‑149 (MQ=255)
                                                                                        aTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATa                                                         <  7:369898/149‑1 (MQ=255)
                                                                                         ttGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAg                                                        >  1:392942/1‑149 (MQ=255)
                                                                                              tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGCAGGCCGATAGCCGAg                                                   >  6:271489/1‑149 (MQ=255)
                                                                                                aTACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGCCCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTc                                                 >  4:68197/1‑149 (MQ=255)
                                                                                                    cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                      <  4:218338/140‑1 (MQ=255)
                                                                                                    cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                      >  3:218338/1‑140 (MQ=255)
                                                                                                    cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                             <  7:201803/149‑1 (MQ=255)
                                                                                                      acacCTCTGCGATGAGGACATGACCTCTTCACCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaacc                                            >  3:19284/1‑148 (MQ=255)
                                                                                                          ctctGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGCTAGCCGAGGCCCACCAATTc                                       >  1:182225/1‑149 (MQ=255)
                                                                                                               cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCTGATAGCCGAGTCCAACCAATTCTAt                                    <  5:338546/147‑1 (MQ=255)
                                                                                                               cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAtt                                   <  3:104252/148‑1 (MQ=255)
                                                                                                                 aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                                >  1:217594/1‑149 (MQ=255)
                                                                                                                 aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                                >  3:72780/1‑149 (MQ=255)
                                                                                                                  tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                               >  3:243418/1‑149 (MQ=255)
                                                                                                                  tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                               <  2:217594/149‑1 (MQ=255)
                                                                                                                    aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                              <  1:295320/148‑1 (MQ=255)
                                                                                                                    aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                              <  5:185502/148‑1 (MQ=255)
                                                                                                                    aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGa                             <  4:375534/149‑1 (MQ=255)
                                                                                                                          tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGa                                                    >  4:335551/1‑120 (MQ=255)
                                                                                                                          tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGc                       <  2:392942/149‑1 (MQ=255)
                                                                                                                           tggctctTCTCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc                      >  8:249695/4‑149 (MQ=38)
                                                                                                                            aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                         >  6:236529/1‑145 (MQ=255)
                                                                                                                            aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                         <  5:236529/145‑1 (MQ=255)
                                                                                                                             ccTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCt                     >  4:103081/1‑148 (MQ=255)
                                                                                                                              ctctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCACCCAATAATATTTTTCGACCTGGCCTGa                   >  6:30807/1‑149 (MQ=255)
                                                                                                                               tctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc                  >  3:361985/1‑149 (MQ=255)
                                                                                                                                  ttctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg                >  7:360367/4‑149 (MQ=255)
                                                                                                                                   tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg                >  3:112802/3‑148 (MQ=255)
                                                                                                                                   tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg                >  3:118277/3‑148 (MQ=255)
                                                                                                                                     tCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt            <  8:68216/149‑1 (MQ=255)
                                                                                                                                        gcgcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga         <  3:381026/149‑1 (MQ=255)
                                                                                                                                          gcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc       >  3:152855/1‑149 (MQ=255)
                                                                                                                                              cTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgtgt   >  2:391848/1‑149 (MQ=255)
                                                                                                                                              cTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgtgt   >  5:397350/1‑149 (MQ=255)
                                                                                                                                              cTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTATCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCgagt   >  4:171660/1‑149 (MQ=255)
                                                                                                                                               ttAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga         >  5:393347/1‑142 (MQ=255)
                                                                                                                                               ttAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga         <  6:393347/142‑1 (MQ=255)
                                                                                                                                                 aaGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGtgtg  <  4:72780/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGCGTGTG  >  NC_002947/5041057‑5041347

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: