Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,155,974 C→T 6.5% A273T (GCT→ACT)  PP_3659 ← transcriptional regulator, AraC family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,155,9740CT 6.5% 181.1 / 2.3 61A273T (GCT→ACT) PP_3659transcriptional regulator, AraC family
Reads supporting (aligned to +/- strand):  ref base C (27/30);  new base T (2/2);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GAATGCCCTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAG  >  NC_002947/4155831‑4156118
                                                                                                                                               |                                                                                                                                                
gAATGCCCTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATa                                                                                                                                               <  5:167857/147‑1 (MQ=255)
 aaTGCCCTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATaga                                                                                                                                             <  7:20598/148‑1 (MQ=255)
 aaTGCCCTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATa                                                                                                                                               >  8:20598/1‑146 (MQ=255)
      ccTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTc                                                                                                                                       >  1:79046/1‑149 (MQ=255)
          tgGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCt                                                                                                                                    <  4:19664/148‑1 (MQ=255)
           ggAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTg                                                                                                                                   >  5:230817/1‑148 (MQ=255)
            gtaCGCACTGTTCTCAATGAATCCAACCTCTGACGAAACGCCATCCACTGTCATACCTGCCAGGCTCATCAGACTGATGGAAATATCCCGCCTTAGTTTATCTTTGTGCCGTTGGTAGGTGAAACCCTCAGCATAGATCCTTCGCTGTa                                                                                                                                 <  1:247415/147‑1 (MQ=255)
            gAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTa                                                                                                                                 <  3:94396/149‑1 (MQ=255)
            gAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTa                                                                                                                                 <  3:237109/149‑1 (MQ=255)
             aaCGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGATGTAg                                                                                                                                >  3:264080/1‑149 (MQ=255)
                gCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGt                                                                                                                             <  1:17963/149‑1 (MQ=255)
                 cACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTg                                                                                                                            <  8:98777/149‑1 (MQ=255)
                  aCTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTgg                                                                                                                           <  2:103145/149‑1 (MQ=255)
                         tcAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGa                                                                                                                    <  4:264080/149‑1 (MQ=255)
                           aaTGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTgg                                                                                                                     <  5:120556/146‑1 (MQ=255)
                           aaTGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGAc                                                                                                                   <  6:230817/148‑1 (MQ=255)
                                aGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCt                                                                                                              <  8:159684/148‑1 (MQ=255)
                                       ctctGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTcgc                                                                                                      >  5:134887/1‑149 (MQ=255)
                                         ctGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCgg                                                                                                    >  4:231279/1‑149 (MQ=255)
                                          tGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGTATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCggg                                                                                                   >  6:214617/1‑149 (MQ=255)
                                            aCGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCgg                                                                                                    <  3:231279/146‑1 (MQ=255)
                                            aCGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGc                                                                                                  <  2:79046/148‑1 (MQ=255)
                                              gCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGTATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAg                                                                                               <  5:214617/149‑1 (MQ=255)
                                                aaCGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTgg                                                                                                                     <  1:83939/125‑1 (MQ=255)
                                                aaCGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTgg                                                                                                                     >  2:83939/1‑125 (MQ=255)
                                                       tCCACTGTCATAGCTGCCCGTCAAAGCAGACTGATGGAAATATCCCGACTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGGTGTAGGGTGGAAGTGGACATGCTCAAGTCGCGGGCGATCTCTTCGAg                                                                                      <  1:165790/149‑1 (MQ=255)
                                                        ccACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAg                                                                                                          >  1:10206/1‑128 (MQ=255)
                                                        ccACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAg                                                                                                          <  2:10206/128‑1 (MQ=255)
                                                         cACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGtt                                                                                    <  1:236522/149‑1 (MQ=255)
                                                          aCTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGtt                                                                                    >  1:116733/1‑148 (MQ=255)
                                                            tGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGTATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAg                                                                                 >  8:188075/1‑149 (MQ=255)
                                                                       cccGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACgtgtg                                                                       <  1:293458/148‑1 (MQ=255)
                                                                        ccGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATc                                                                                         >  1:277261/1‑129 (MQ=255)
                                                                        ccGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATc                                                                                         <  2:277261/129‑1 (MQ=255)
                                                                            cacaGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGgaga                                                                 <  2:116733/149‑1 (MQ=255)
                                                                             agagcagACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGgagag                                                                <  1:122842/147‑1 (MQ=255)
                                                                               agcagACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGagagag                                                              >  2:14468/1‑149 (MQ=255)
                                                                                      tGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATg                                                                                                                >  6:128123/1‑92 (MQ=255)
                                                                                        aTGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTaa                                                     >  8:20728/1‑149 (MQ=255)
                                                                                          ggAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTGGCGGGCGAGCTCATCGAGTTCAg                                                                                 >  1:280812/1‑119 (MQ=255)
                                                                                          ggAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAg                                                                                 <  2:280812/119‑1 (MQ=255)
                                                                                              atatCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATcgccg                                               >  6:159854/1‑149 (MQ=255)
                                                                                                   ccGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATcgccgc                                              <  5:232372/145‑1 (MQ=255)
                                                                                                   ccGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATcgccgc                                              >  6:232372/1‑145 (MQ=255)
                                                                                                    cGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGAt                                         <  7:20728/149‑1 (MQ=255)
                                                                                                              tATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTc                               >  8:213627/1‑149 (MQ=255)
                                                                                                               aTCTTTGAGCCGTTGGTAGGTGAAACCCTCAGTATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTc                               <  7:188075/148‑1 (MQ=255)
                                                                                                               aTCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCg                              >  1:45741/1‑149 (MQ=255)
                                                                                                                  tttGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGAt                                         >  2:78669/1‑135 (MQ=255)
                                                                                                                  tttGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGGTTTAATCGCCGCCGGAt                                         <  1:78669/135‑1 (MQ=255)
                                                                                                                        ccGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTgtcgt                     >  3:26955/1‑149 (MQ=255)
                                                                                                                           ttGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTcgtcgt                  >  6:221943/1‑149 (MQ=255)
                                                                                                                             ggTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGtt                 <  5:221943/148‑1 (MQ=255)
                                                                                                                              gTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTAc               >  1:132219/1‑149 (MQ=255)
                                                                                                                               tAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGtt                 >  7:296716/1‑146 (MQ=255)
                                                                                                                               tAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGtt                 <  8:296716/146‑1 (MQ=255)
                                                                                                                               tAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACg              >  8:190959/1‑149 (MQ=255)
                                                                                                                                     aaaCCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGATAGCAGTTTTAATCGCTGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTAtt         <  6:134887/148‑1 (MQ=255)
                                                                                                                                         ccTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTcac     >  2:326977/1‑148 (MQ=255)
                                                                                                                                           tCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAg  >  5:340578/1‑149 (MQ=255)
                                                                                                                                            cAGCATAGACCCTTCGCTGTAGGGTGGAACTGGACATACTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGAGGAGAGACCAGATTATAACACAGCCGGATCATAATGATCGTACTGTCGTCGTTACGGTATTTCACAAg  >  1:284400/1‑148 (MQ=255)
                                                                                                                                               cATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTc       >  6:129631/1‑140 (MQ=255)
                                                                                                                                               cATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTc       <  5:129631/140‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
GAATGCCCTGTGGAACGCACTGGTCTCAATGAAGCCAACCTCTGACGCAACGCCATCCACTGTCATACCTGCCCGGCACAGCAGACTGATGGAAATATCCCGCCTTAGATTATCTTTGAGCCGTTGGTAGGTGAAACCCTCAGCATAGAGCCTTCGCTGTAGGGTGGAACTGGACATGCTCAAGTCGCGGGCGAGCTCATCGAGTTCAGGCCACGTGTGTGGAGAGAGCAGTTTTAATCGCCGCCGGATCATAATGCTCGTACTGTCGTCGTTACGGTATTTCACAAG  >  NC_002947/4155831‑4156118

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: ‑