Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,900,065 A→G 7.7% intergenic (+26/‑472) PP_4308 → / → PP_4309 AsnC family transcriptional regulator/NCS1 family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,900,0650AG 7.7% 185.1 / 4.2 65intergenic (+26/‑472)PP_4308/PP_4309AsnC family transcriptional regulator/NCS1 family transporter
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base G (3/2);  total (33/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.97e-01

CATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGA  >  NC_002947/4899920‑4900211
                                                                                                                                                 |                                                                                                                                                   
cATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGc                                                                                                                                                  >  8:103103/1‑149 (MQ=255)
cATGTGCGCTCGATGGAGCAGCTGCACACCCTGCTCGACCTCATGAATGTGCTGGCCGACACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAg                                                                                                                                                    <  2:315023/147‑1 (MQ=255)
    tgCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGc                                                                                                                                              >  5:325757/1‑149 (MQ=255)
           gATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGACTCCTGTGCTGGCCTCTTCGCAGGCGTGCCcgcg                                                                                                                                         >  4:65448/1‑147 (MQ=255)
             tGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                     <  5:298150/149‑1 (MQ=255)
              ggAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAagag                                                                                                                                    <  6:309252/149‑1 (MQ=255)
                agcagcTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGc                                                                                                                                  >  6:197249/1‑149 (MQ=255)
                             ccggcgcGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                       <  2:82383/125‑1 (MQ=255)
                             ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                       >  6:165054/1‑147 (MQ=255)
                             ccTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGaa                                                                                                                       >  1:189560/1‑147 (MQ=255)
                                    gACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGACTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGa                                                                                                                <  3:65448/147‑1 (MQ=255)
                                     aCCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTAGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTa                                                                                                            >  5:75598/1‑149 (MQ=255)
                                       cTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACGCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCAAGTAAAGGAAAAGATGAGTagt                                                                                                          >  7:264302/1‑147 (MQ=255)
                                        tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                       <  4:94360/120‑1 (MQ=255)
                                        tCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGaa                                                                                                                                       >  3:94360/1‑120 (MQ=255)
                                          aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAg                                                                                                                                                    <  3:2545/105‑1 (MQ=255)
                                          aTGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAg                                                                                                                                                    >  4:2545/1‑105 (MQ=255)
                                                    tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                                   >  7:218368/1‑143 (MQ=255)
                                                    tGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTa                                                                                                   <  8:218368/143‑1 (MQ=255)
                                                      gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAgg                                                                                                                                   <  2:46663/110‑1 (MQ=255)
                                                      gCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAgg                                                                                                                                   >  1:46663/1‑110 (MQ=255)
                                                        cGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCa                                                                                                                                <  7:244640/111‑1 (MQ=255)
                                                        cGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCa                                                                                                                                >  8:244640/1‑111 (MQ=255)
                                                         gagaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACg                                                                                         >  2:170767/1‑148 (MQ=255)
                                                          agaCCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGt                                                                                        >  4:107380/1‑148 (MQ=255)
                                                                     gcgATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTAAACAAAAACGTAACCATAgtggc                                                                            >  8:288990/1‑145 (MQ=255)
                                                                       gATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                     >  6:314115/1‑91 (MQ=255)
                                                                       gATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAaga                                                                                                                                     <  5:314115/91‑1 (MQ=255)
                                                                           atcaAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAtttt                                                                      >  2:83602/1‑149 (MQ=255)
                                                                              aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                   <  6:325757/149‑1 (MQ=255)
                                                                              aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                   <  2:189560/149‑1 (MQ=255)
                                                                              aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                   <  5:197249/149‑1 (MQ=255)
                                                                              aagaagACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCAGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgt                                                                   <  6:75598/149‑1 (MQ=255)
                                                                                gaagaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTtgtgt                                                                 >  3:141087/1‑149 (MQ=255)
                                                                                   gaCCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAACAGGCCAGTACAGGAAAAACTGAGGTATTTTTTTTACACAAAAACGTAACCATAGGTGCCATTTTtgtgtg                                                                >  1:102252/1‑148 (MQ=39)
                                                                                      cccGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTCAGAATTTTTTTTAAACAAAAAAGTAACCCTAGGTTCCATTTTTTTTTTCACtt                                                           >  1:260775/1‑149 (MQ=38)
                                                                                        cGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTttc                                                         <  7:103103/149‑1 (MQ=255)
                                                                                             aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAg                                                    <  1:170767/149‑1 (MQ=255)
                                                                                             aaGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAg                                                    >  2:96361/1‑149 (MQ=255)
                                                                                                cggcggTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTgg                                                 <  1:83602/149‑1 (MQ=255)
                                                                                                  gcggTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGgcg                                               <  4:141087/149‑1 (MQ=255)
                                                                                                          cgccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAAtt                                       <  5:122042/149‑1 (MQ=255)
                                                                                                          caccAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAAtt                                       >  8:20510/3‑149 (MQ=255)
                                                                                                              aaTGGATTGATTGGCTCCTGTGCTTGCCTCTTCGCAGGCGTGCCCGCGAAGATGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCttt                                   <  7:7860/149‑1 (MQ=255)
                                                                                                               aTGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCTCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                  <  6:201802/149‑1 (MQ=255)
                                                                                                                   attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAtttt                                                                      <  1:183140/109‑1 (MQ=255)
                                                                                                                   attgattgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCAtttt                                                                      >  2:183140/1‑109 (MQ=255)
                                                                                                                        ttgGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTc                         >  4:246396/1‑149 (MQ=255)
                                                                                                                          ggCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTcc                        <  5:165054/148‑1 (MQ=255)
                                                                                                                           gCTCCTGTGCTGGCCTGTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACATGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCa                       <  7:288990/148‑1 (MQ=255)
                                                                                                                                 gtgCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCTTAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTa                >  4:300602/1‑149 (MQ=255)
                                                                                                                                   gCTGGCCTCTTCGCAGGCGTGCCCTCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                                   <  1:232221/128‑1 (MQ=255)
                                                                                                                                   gCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGtt                                   >  2:232221/1‑128 (MQ=255)
                                                                                                                                       gCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTAtttt            <  7:20510/147‑1 (MQ=255)
                                                                                                                                        ccTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCTTAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc         <  3:300602/149‑1 (MQ=255)
                                                                                                                                         ctctTCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                  <  6:247527/123‑1 (MQ=38)
                                                                                                                                         ctctTCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTa                                  >  5:247527/1‑123 (MQ=38)
                                                                                                                                         ctctTCGCAGGCGTGTCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                                    >  3:293884/1‑121 (MQ=255)
                                                                                                                                         ctctTCGCAGGCGTGTCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGt                                    <  4:293884/121‑1 (MQ=255)
                                                                                                                                         ctctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAc         >  7:225893/1‑148 (MQ=255)
                                                                                                                                          tctTCGCGGGCACGCCTGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTATTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc        >  5:210018/1‑148 (MQ=255)
                                                                                                                                          tctTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGAcc        >  7:233890/1‑148 (MQ=255)
                                                                                                                                            ttCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCaa                                         >  6:22055/1‑113 (MQ=38)
                                                                                                                                            ttCGCGGGCACGCCTGCGAAGAGGCCAGCACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCaa                                         <  5:22055/113‑1 (MQ=38)
                                                                                                                                               gCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGa  <  3:246396/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CATGTGCGCTCGATGGAGCAGCTGGACACCCTGCTCGACCTCATGAATGTGCTGGCCGAGACCAATACCGCGATCATCAAGAAGACCCCGGTCAAGCGGCGGTTGCCGCCAATGGATTGATTGGCTCCTGTGCTGGCCTCTTCGCAGGCGTGCCCGCGAAGAGGCCAGTACAGGAAAAGCTGAGTATTTTTTGTACACAAAAACGTAACCATAGGTGCCATTTTTGTGTGCACTTTCTTCAGTGGCGCCCCAATTCGTTACACACTTTCCACCCTCTATTTTGACCAAATGA  >  NC_002947/4899920‑4900211

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