Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,765,248 C→A 6.1% T138K (ACG→AAG) ‡ PP_1574 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,765,2480CA 6.1% 212.2 / 3.1 65T138K (ACG→AAG) ‡PP_1574hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/29);  new base A (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

ACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGCCTG  >  NC_002947/1765101‑1765395
                                                                                                                                                   |                                                                                                                                                   
aCCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCAc                                                                                                                                                     <  5:215040/148‑1 (MQ=255)
  ccGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAg                                                                                                                                                  >  8:167019/1‑149 (MQ=255)
              gagaAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTcc                                                                                                                                      >  7:251570/1‑149 (MQ=255)
                       ccGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATcggc                                                                                                                             >  6:92005/1‑149 (MQ=255)
                       ccGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATcggc                                                                                                                             >  4:217239/1‑149 (MQ=255)
                          gagaGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTc                                                                                                                                            <  8:87853/131‑1 (MQ=255)
                                aCCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgac                                                                                                                     >  4:117657/1‑148 (MQ=255)
                                  cAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTcgacga                                                                                                                   >  8:5607/1‑148 (MQ=255)
                                    aGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcg                                                                                                                >  6:54079/1‑149 (MQ=255)
                                    aGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATcg                                                                                                                >  3:104080/1‑149 (MQ=255)
                                               tctcTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTc                                                                                                                                            <  4:199003/110‑1 (MQ=255)
                                               tctcTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTc                                                                                                                                            >  3:199003/1‑110 (MQ=255)
                                                 tctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAg                                                                                                    >  4:49201/1‑148 (MQ=255)
                                                 tctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAg                                                                                                    >  7:91008/1‑148 (MQ=255)
                                                 tctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAg                                                                                                    >  3:26683/1‑148 (MQ=255)
                                                  ctTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGa                                                                                                   >  6:277395/1‑148 (MQ=255)
                                                     gATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAc                                                                                               >  8:135599/1‑149 (MQ=255)
                                                      aTGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAAcc                                                                                              <  4:104080/149‑1 (MQ=255)
                                                       tGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCt                                                                                             <  4:26683/149‑1 (MQ=255)
                                                         gCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGCGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  7:167019/149‑1 (MQ=255)
                                                         gCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  6:93832/149‑1 (MQ=255)
                                                         gCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  3:117657/149‑1 (MQ=255)
                                                         gCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCCGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  3:257799/149‑1 (MQ=255)
                                                         gCGACCTCTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  4:288498/149‑1 (MQ=255)
                                                              cTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTg                                                                                      >  2:180136/1‑149 (MQ=255)
                                                                 tCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGt                                                                                   <  3:217239/149‑1 (MQ=255)
                                                                    gTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGt                                                                                   <  3:218252/146‑1 (MQ=255)
                                                                    gTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGt                                                                                   >  4:218252/1‑146 (MQ=255)
                                                                       gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGtt                                                                              >  1:88550/1‑148 (MQ=255)
                                                                       gcgcTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTc                                                                             <  7:135599/149‑1 (MQ=255)
                                                                           tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                           >  3:58443/1‑147 (MQ=255)
                                                                           tgctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAg                                                                           <  4:58443/147‑1 (MQ=255)
                                                                             ctgGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCt                                                                       >  7:6318/1‑149 (MQ=255)
                                                                                gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGt                                                                               <  2:234964/138‑1 (MQ=255)
                                                                                gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGt                                                                               >  1:234964/1‑138 (MQ=255)
                                                                                gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTgg                                                                    <  3:262023/149‑1 (MQ=255)
                                                                                gTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAACAAAACCTGAGCCTGCGTGAGTTCAGGTCTTgg                                                                    >  2:284175/1‑149 (MQ=255)
                                                                                     tATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgag                                                                 >  1:67633/1‑147 (MQ=255)
                                                                                      aTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                    <  4:211740/127‑1 (MQ=255)
                                                                                      aTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCg                                                                                    >  3:211740/1‑127 (MQ=255)
                                                                                      aTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGgag                                                                 >  8:72529/1‑146 (MQ=255)
                                                                                      aTTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGaa                                                               <  7:72529/148‑1 (MQ=255)
                                                                                       ttCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAg                                                              <  8:91008/148‑1 (MQ=255)
                                                                                          gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTAc                                                           <  5:277395/148‑1 (MQ=255)
                                                                                          gggAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACa                                                          <  5:92005/149‑1 (MQ=255)
                                                                                          gggAGGTCAACGACCGGTGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACa                                                          <  8:288525/149‑1 (MQ=255)
                                                                                               gTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                            >  2:82196/1‑126 (MQ=255)
                                                                                               gTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCa                                                                            <  1:82196/126‑1 (MQ=255)
                                                                                                   aCGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           >  4:174407/1‑107 (MQ=255)
                                                                                                     gACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAACAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGt                                               <  1:284175/149‑1 (MQ=255)
                                                                                                            gAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGCGAGAAGTACAGGCAGCTCCGTGGCTCGc                                        <  1:180136/149‑1 (MQ=255)
                                                                                                             aaaGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGt                                               <  6:296704/141‑1 (MQ=255)
                                                                                                             aaaGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGt                                               >  5:296704/1‑141 (MQ=255)
                                                                                                                            aCCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATgcg                        >  3:207114/1‑149 (MQ=255)
                                                                                                                                 gCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTg                    >  7:34622/1‑148 (MQ=255)
                                                                                                                                     gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGa                            >  8:156903/1‑136 (MQ=255)
                                                                                                                                     gacgaGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGa                            <  7:156903/136‑1 (MQ=255)
                                                                                                                                     gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           >  8:67789/1‑73 (MQ=37)
                                                                                                                                     gacgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGa                                                                                           <  7:67789/73‑1 (MQ=37)
                                                                                                                                      acgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATg                                  <  2:302383/129‑1 (MQ=255)
                                                                                                                                      acgaGCTCGTGCCAGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATg                                  >  1:302383/1‑129 (MQ=255)
                                                                                                                                           ctctGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTg                                              <  2:206924/112‑1 (MQ=255)
                                                                                                                                           ctctGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTg                                              >  1:206924/1‑112 (MQ=255)
                                                                                                                                            tctGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACg        >  4:75998/1‑149 (MQ=255)
                                                                                                                                                gCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGc     <  8:6318/148‑1 (MQ=255)
                                                                                                                                                  aCGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGCctg  <  5:54079/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACCCGGTTGAACTGGAGAAGGACCCGGAGAGCACCAAGCGCCTGGGATCTCTTGATGGCGACCTGTCCGTCGCGCTGCTGGTAGCTATTCGGGAGGTCAACGACCTGGGAAAGACGAAGAGCTCACCGAGCTCGACGAGCTCTGGCACGAGCTCGTCCTCTCCGGCATCGGCGGCTCGACGATCGCGCAGGCCAAAGAAAACCTGAGCCTGCGTGAGTTCAGGTCTTGGGAGAAGTACAGGCAGCTCCGTGGCTCGCTGCATGTCGGGATGCGCGTGGAGCGCGCCACGGGCCTG  >  NC_002947/1765101‑1765395

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