Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,508,967 G→T 6.3% T184T (ACG→ACT)  PP_3106 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,508,9670GT 6.3% 200.1 / 2.1 63T184T (ACG→ACT) PP_3106hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/32);  new base T (2/2);  total (29/34)
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.48e-01

GCCGGATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGCTGCAATCCGG  >  NC_002947/3508818‑3509116
                                                                                                                                                     |                                                                                                                                                     
gCCGGATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACg                                                                                                                                                       <  1:2121936/150‑1 (MQ=255)
 ccGGATTTGGCACGCTGGCAAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCACACAGCTACCTTACCCCGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACgg                                                                                                                                                      >  1:527672/1‑150 (MQ=255)
  cGGATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGc                                                                                                                                                     <  2:1491640/150‑1 (MQ=255)
  cGGATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTGCGAGCACACGGc                                                                                                                                                     <  1:228673/150‑1 (MQ=255)
   ggATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGcc                                                                                                                                                    >  1:1274459/1‑150 (MQ=255)
       ttGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACgg                                                                                                                                                      <  1:1262730/144‑1 (MQ=255)
       ttGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACgg                                                                                                                                                      >  2:1262730/1‑144 (MQ=255)
                tGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGa                                                                                                                                       <  2:1164926/150‑1 (MQ=255)
                tGGCCAAGTACGAATTTCGTGTAAGCAACCCGAACAAATCGAACAGCTACATTACCCAGTACCGCGACGTCGATTTCAACTGTCTGCAGCGCTTGCTTGAGAGCGAAGCGATCTTCTACGACTTCGACCACACGGCCGACGCCCACCTgt                                                                                                                                       <  2:527672/150‑2 (MQ=255)
                    cAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGt                                                                                                                                   >  2:875436/1‑150 (MQ=255)
                     aaGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGtc                                                                                                                                  <  2:708098/150‑1 (MQ=255)
                      aGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGtca                                                                                                                                 <  1:875436/150‑1 (MQ=255)
                              tttCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         <  2:2146979/150‑1 (MQ=255)
                              tttCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTCTCTTCAGTGCTTTCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         <  2:573481/150‑1 (MQ=255)
                                         caagcCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCAc                                                                                                              >  1:1262559/1‑150 (MQ=255)
                                                cTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCttgttg                                                                                                       <  1:1690748/150‑1 (MQ=255)
                                                 tCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGc                                                                                                      >  1:504740/1‑150 (MQ=255)
                                                    aaTCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACTGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGc                                                                                                      >  1:218925/1‑147 (MQ=255)
                                                    aaTCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACTGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGc                                                                                                      >  1:218620/1‑147 (MQ=255)
                                                    aaTCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACTGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGc                                                                                                      <  2:218620/147‑1 (MQ=255)
                                                    aaTCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACTGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGc                                                                                                      <  2:218925/147‑1 (MQ=255)
                                                      tCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTg                                                                                                 >  1:1459773/1‑150 (MQ=255)
                                                         cacaGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGccc                                                                                              <  1:880728/150‑1 (MQ=255)
                                                         cacaGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGccc                                                                                              <  1:1284003/150‑1 (MQ=255)
                                                                  aTTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCAAGTCAGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCg                                                                                     >  1:126682/1‑150 (MQ=255)
                                                                       ccAGTACCGCGAGAGCGATTTCAACTTNGTGCAGCGNTTGNTTNAGANCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGAt                                                                                >  1:2190186/1‑150 (MQ=255)
                                                                         aGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGAttg                                                                              >  1:954047/1‑149 (MQ=255)
                                                                           tACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGc                                                                            <  2:2237714/150‑1 (MQ=255)
                                                                                        aTTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGt                                                               <  2:1274459/150‑1 (MQ=255)
                                                                                        aTTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGt                                                               >  1:623764/1‑150 (MQ=255)
                                                                                         tttCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTc                                                              >  2:1509191/1‑150 (MQ=255)
                                                                                                tttGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTAcc                                                       >  1:606664/1‑150 (MQ=255)
                                                                                                tttGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTAcc                                                       >  1:2296974/1‑150 (MQ=255)
                                                                                                      cAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTcc                                                                                                                                              <  1:874501/57‑1 (MQ=255)
                                                                                                      cAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTcc                                                                                                                                              >  2:874501/1‑57 (MQ=255)
                                                                                                       aGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCg                                                >  2:542600/1‑150 (MQ=255)
                                                                                                       aGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCg                                                >  2:138970/1‑150 (MQ=255)
                                                                                                        gcgcTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGc                                               <  1:2194786/150‑1 (MQ=255)
                                                                                                          gcttgcttGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTg                                             <  2:1262559/150‑1 (MQ=255)
                                                                                                           cttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGaa                                                                                           <  2:201603/103‑1 (MQ=255)
                                                                                                           cttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGaa                                                                                           >  1:201603/1‑103 (MQ=255)
                                                                                                           cttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGa                                            <  2:2296974/150‑1 (MQ=255)
                                                                                                            ttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         >  1:819233/1‑72 (MQ=255)
                                                                                                            ttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         <  2:819233/72‑1 (MQ=255)
                                                                                                            ttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCAcc                                                                                                             >  1:870532/1‑84 (MQ=255)
                                                                                                            ttgcttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCAcc                                                                                                             <  2:870532/84‑1 (MQ=255)
                                                                                                            ttgcttgAGAGCGAAGGGCTTTTATACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATAACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGAt                                           >  1:1154595/1‑150 (MQ=255)
                                                                                                            ttgctt‑aGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGAtt                                          >  2:1783988/1‑150 (MQ=255)
                                                                                                                ttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         <  1:2236573/68‑1 (MQ=255)
                                                                                                                ttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGAt                                                                                                                         >  2:2236573/1‑68 (MQ=255)
                                                                                                                ttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCa                                                                                                                  <  2:1601628/75‑1 (MQ=255)
                                                                                                                ttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCa                                                                                                                  >  1:1601628/1‑75 (MQ=255)
                                                                                                                ttgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCa                                       <  2:251930/150‑1 (MQ=255)
                                                                                                                 tgAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCAt                                      <  1:361346/150‑1 (MQ=255)
                                                                                                                          aGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTg                             <  2:504740/150‑1 (MQ=255)
                                                                                                                             gcgtttCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCa                          <  2:232167/147‑1 (MQ=255)
                                                                                                                             gCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCa                          <  1:138970/150‑1 (MQ=255)
                                                                                                                             gCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCCGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCCAAACCGCTGATCCCATTCCTCAGTGGCa                          >  1:601149/1‑150 (MQ=255)
                                                                                                                                    tactacTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAAcg                   >  1:1829088/1‑150 (MQ=255)
                                                                                                                                           tCGAGCACACGGCCGACTCCCACCTGATGGTCATCCCCGACGTCTCCCTCACCTTTTTGCCGCTGCCCGAAAACATTCTGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGc            <  1:542600/150‑1 (MQ=255)
                                                                                                                                                 acacGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGa                                                                                 >  2:1943662/1‑75 (MQ=255)
                                                                                                                                                 acacGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGa                                                                                 <  1:1943662/75‑1 (MQ=255)
                                                                                                                                                 acacGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGCTGCAAt      >  2:2092314/1‑150 (MQ=255)
                                                                                                                                                  cacGGCCGACTCCCCCCGGATGTGCATCGNCGTTGANTCCNGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGCTGCAATc     <  2:2190186/150‑1 (MQ=255)
                                                                                                                                                   acGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGc            <  1:749824/142‑1 (MQ=255)
                                                                                                                                                   acGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGc            >  2:749824/1‑142 (MQ=255)
                                                                                                                                                    cGGCCGACCCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGCTGCAATCCg   <  2:1459773/150‑1 (MQ=255)
                                                                                                                                                     ggCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGACGGTTACAGAAACCGCTGATTCGATTACTCAGTGGCAGTCAACGCGGCAGCTGCAATCCgg  >  1:1969090/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGGATTTGGCACGCTGGCCAAGTACGAATTTCGTGTAAGCAAGCCGCTCAAATCGCACAGCTACATTACCCAGTACCGCGAGAGCGATTTCAACTTTGTGCAGCGCTTGCTTGAGAGCGAAGGGCTTTTCTACTACTTCGAGCACACGGCCGACTCCCACCTGATGGTCATCACCGATGACTCCAGCACCTTGTTGCCGCTGCCCGAACAGCCGCAGATTCGCTACCACTCCGCGTCGGTTACCGAAACCGCTGATTCCATTACTCAGTGGCAGTCAACGCGGCAGCTGCAATCCGG  >  NC_002947/3508818‑3509116

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑