Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,066 A→T 10.8% intergenic (+31/+260) PP_2003 → / ← PP_2004 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,0660AT10.8% 163.2 / 11.0 65intergenic (+31/+260)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (27/31);  new base T (2/5);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.49e-01

GTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAG  >  NC_002947/2272925‑2273214
                                                                                                                                               |                                                                                                                                                    
gtggtgGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcgcg                                                                                                                                               <  6:255148/149‑1 (MQ=255)
 tggtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCcgcg                                                                                                                                               >  6:138397/1‑148 (MQ=255)
   gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                                            <  2:382986/149‑1 (MQ=255)
   gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGaaa                                                                                                                                            <  1:394918/149‑1 (MQ=255)
    tggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAg                                                                                                                                           <  2:119182/149‑1 (MQ=255)
       aTGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGc                                                                                                                                        <  3:370332/149‑1 (MQ=255)
        tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAgggtc                                                                                                                                       <  3:581/149‑3 (MQ=255)
        tGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGcc                                                                                                                                       >  1:166985/1‑149 (MQ=255)
         gACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                                      >  6:401334/1‑149 (MQ=255)
          aCAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCg                                                                                                                                      <  5:138397/148‑1 (MQ=255)
            aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCa                                                                                                                                                      <  4:262791/130‑1 (MQ=255)
            aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCa                                                                                                                                                      >  3:262791/1‑130 (MQ=255)
            aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                                                   >  3:57539/1‑149 (MQ=255)
                aGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGttaacg                                                                                                                               >  2:63213/1‑145 (MQ=255)
                 gggTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGTTTACAGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAg                                                                                                                                <  3:249758/149‑1 (MQ=25)
                  ggctCGCAGCCTGCGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCt                                                                                                                             <  4:231644/146‑1 (MQ=255)
                               tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATccc                                                                                                                  <  3:105192/147‑1 (MQ=255)
                                   cccGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTa                                                                                                            <  5:69871/149‑1 (MQ=255)
                                          gccgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCAGAACCCACCCATTAGCCCGCACCCCCCGAAAAGGCCCGTTCAGGGGGTTAACAAACCCAGTTAAAGAccc                                                                                                     >  2:143825/1‑147 (MQ=255)
                                            cgAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGccc                                                                                                   >  5:117788/1‑149 (MQ=255)
                                                    cACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCCACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATa                                                                                            >  8:383165/1‑148 (MQ=255)
                                                       ccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCATTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATGCCCGGTTACCGACGCCCCGCGATAACtt                                                                                        >  3:347249/1‑149 (MQ=255)
                                                       ccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAAACCACCCATTAGCCGGCAGCCCGCGAAAGGGCCGGTTGAGGCTGTAAACATCCCCGGTTACAGACGCCCCGCGATAAAtt                                                                                        >  1:153100/1‑149 (MQ=255)
                                                         aTTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTc                                                                                      >  1:108363/1‑149 (MQ=255)
                                                          ttCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCt                                                                                     <  3:301992/149‑1 (MQ=255)
                                                            cAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTg                                                                                    <  7:18536/148‑1 (MQ=255)
                                                              ggTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAac                                                                                 <  5:401334/149‑1 (MQ=255)
                                                               gTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCACCTGTTACAGACGCCCCGCGATAACTTTCTGAaca                                                                                <  2:166985/149‑1 (MQ=255)
                                                                          tGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACTCCTGCTACGt                                                                      <  6:117788/148‑1 (MQ=255)
                                                                            cccAGCCGGCCGAGGCCGCCGCCCCCTCACTTTGATGTGTCGCCTGAACCGACCCTTTCGCCGTCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGttt                                                                    <  3:280881/148‑1 (MQ=255)
                                                                             ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                                            <  3:292433/75‑1 (MQ=25)
                                                                             ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                                            >  4:292433/1‑75 (MQ=25)
                                                                             ccAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACttt                                                                                       <  7:75347/128‑1 (MQ=255)
                                                                                  cggccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTc                                                             <  7:383165/149‑1 (MQ=255)
                                                                                    gccgAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCAt                                                           <  1:143825/149‑1 (MQ=255)
                                                                                          gccgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCtt                                                     <  4:57539/149‑1 (MQ=255)
                                                                                            cgccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTa                                                                                                            <  5:16924/92‑1 (MQ=255)
                                                                                             gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTTCTACGTTTACCACTCATAAGCTTGtc                                                  <  6:412030/149‑1 (MQ=255)
                                                                                             gccgccCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtc                                                  <  2:153100/149‑1 (MQ=255)
                                                                                                 cccccTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatca                                              >  8:321343/1‑149 (MQ=255)
                                                                                                 cccccTCCCTTTGATGGGTCGCCTGAAACGACCCTTTCGCCGGCAGCCCGCGAAATGGCCGGTTCAGCCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatca                                              <  7:104129/149‑1 (MQ=255)
                                                                                                       gccgttGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCACCGCGATAACTTTCTGAACACCTGCTAAGTTTACCACTCATAAGCTTGTCATCATCCCt                                         <  1:349847/144‑1 (MQ=255)
                                                                                                        cctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAACGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCtt                                        >  2:881/1‑148 (MQ=255)
                                                                                                         ctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCTCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCa                                      >  3:161743/1‑149 (MQ=255)
                                                                                                         ctTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGt                                                   <  7:294174/136‑1 (MQ=255)
                                                                                                             gATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATccc                                  >  2:121134/1‑149 (MQ=255)
                                                                                                                 ggTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGTGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCt                              <  1:63213/149‑1 (MQ=255)
                                                                                                                   tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTa                                                                   <  8:95269/110‑1 (MQ=255)
                                                                                                                   tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTa                                                                   >  7:95269/1‑110 (MQ=255)
                                                                                                                   tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATACCTTCATCCCACCTGc                            >  8:363060/1‑149 (MQ=255)
                                                                                                                   tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACATGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGc                            >  6:210681/1‑149 (MQ=255)
                                                                                                                   tCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACAGCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGATTACCACTCATAAGCTTGTCATCATCCCTTCATCCCCCCTGc                            >  7:182961/1‑149 (MQ=255)
                                                                                                                        ccTGAACCGACCCTTTCTCCGTCAGCCCGCGAAAGGGCCGGTTAAGGCTGTAAAAATCACCTGTTACAGACGCAACGCGATAACTTTCTGAACACCTGCTACTTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCtt                          <  3:43878/148‑1 (MQ=255)
                                                                                                                        ccTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGTCCGGTTCAGGCTGTAAAAATCCCCTTTTACAGACGCCCAGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTAATCATCACTTCATACCACATGCtt                          >  4:43878/1‑148 (MQ=255)
                                                                                                                         cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATc                                    >  3:56355/1‑137 (MQ=255)
                                                                                                                         cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATc                                    <  4:56355/137‑1 (MQ=255)
                                                                                                                         cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg                        >  8:245744/1‑149 (MQ=255)
                                                                                                                               cGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTg                    >  5:199236/1‑147 (MQ=255)
                                                                                                                                ggcccTTTCGCGGGCTGCCGGCGAAAGGGTCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACAAATCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTAAGGCTGaat                 >  2:102433/3‑148 (MQ=38)
                                                                                                                                 aCCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAgg                >  2:366647/1‑149 (MQ=255)
                                                                                                                                 aCCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAgg                >  2:15496/1‑149 (MQ=255)
                                                                                                                                    cTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGa             <  1:366647/149‑1 (MQ=255)
                                                                                                                                     tttCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatc                                               >  7:241996/1‑114 (MQ=255)
                                                                                                                                     tttCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGtcatc                                               <  8:241996/114‑1 (MQ=255)
                                                                                                                                         gCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATACCTTCATCCCACCTGCTTAGGCTGAAGGCGATAcc         >  3:210048/1‑148 (MQ=255)
                                                                                                                                          ccGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTa       <  3:86268/149‑1 (MQ=255)
                                                                                                                                              cAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACaaa   <  4:347249/149‑1 (MQ=255)
                                                                                                                                               aGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg  <  7:338607/149‑1 (MQ=255)
                                                                                                                                               aGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg  >  8:384648/1‑149 (MQ=255)
                                                                                                                                               aGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAg  <  7:321343/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GTGGTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAGGCTGAAGGCGATACCTACAAAG  >  NC_002947/2272925‑2273214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: