Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,634,566 Δ1 bp 6.2% coding (985/1872 nt) PP_2304 → peptidyl‑prolyl cis‑trans isomerase D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,634,5660G. 6.2% 208.2 / 6.9 65coding (985/1872 nt)PP_2304peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base G (28/33);  new base . (2/2);  total (30/35)
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.92e-01

ACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACAAA  >  NC_002947/2634420‑2634709
                                                                                                                                                   |                                                                                                                                               
aCGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGAcc                                                                                                                                              <  2:138120/150‑1 (MQ=255)
     cAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGcc                                                                                                                                         >  2:876659/1‑150 (MQ=255)
      aGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCt                                                                                                                                        >  1:540825/1‑150 (MQ=255)
      aGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCt                                                                                                                                        >  2:33613/1‑150 (MQ=255)
           aaGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGaa                                                                                                                                   <  2:1344870/150‑1 (MQ=255)
            aGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAag                                                                                                                                  <  2:896344/150‑1 (MQ=255)
             ggCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTc                                                                                                                                      >  2:10967/1‑145 (MQ=255)
             ggCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTc                                                                                                                                      <  1:10967/145‑1 (MQ=255)
             ggCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGaa                                                                                                                                   <  1:588018/148‑1 (MQ=255)
             ggCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGaa                                                                                                                                   >  2:588018/1‑148 (MQ=255)
              gCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAagag                                                                                                                                <  2:874560/150‑1 (MQ=255)
               cccGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCgg                                                                                                                                           <  2:1613997/138‑1 (MQ=255)
               cccGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCgg                                                                                                                                           >  1:1613997/1‑138 (MQ=255)
                ccGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGgc                                                                                                                              >  2:1225866/1‑150 (MQ=255)
                      cGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTa                                                                                                                        >  2:646385/1‑150 (MQ=255)
                          gagaTCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAg                                                                                                                    <  1:1592795/150‑1 (MQ=255)
                          gagaTCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGACTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCGGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAg                                                                                                                    <  1:1567283/150‑1 (MQ=255)
                              tCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGc                                                                                                                >  1:17924/1‑150 (MQ=255)
                                     gcGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATgg                                                                                                         >  2:50641/1‑150 (MQ=255)
                                      cGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGAGCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGTCCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATggg                                                                                                        <  1:1565181/150‑1 (MQ=255)
                                          tGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAgg                                                                                                    <  2:705960/150‑1 (MQ=255)
                                                         aCTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTgg                                                                                     <  1:876659/150‑2 (MQ=255)
                                                           tttGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATg                                                                                                          >  2:1095838/1‑127 (MQ=255)
                                                           tttGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATg                                                                                                          <  1:1095838/127‑1 (MQ=255)
                                                               ccGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCg                                                                               >  2:1155872/1‑150 (MQ=255)
                                                                     tGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTg                                                                         <  2:540825/150‑1 (MQ=255)
                                                                      ggCCAAGGAGTTATCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTgg                                                                        <  1:1155872/150‑1 (MQ=255)
                                                                             gAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCAt                                                                 >  2:31157/1‑150 (MQ=255)
                                                                               gTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAagag                                                                                                                                <  2:1401531/85‑1 (MQ=255)
                                                                               gTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAagag                                                                                                                                >  1:1401531/1‑85 (MQ=255)
                                                                               gTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCt                                                               <  1:1225866/150‑1 (MQ=255)
                                                                                 tctctcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGa                                                             <  2:17924/150‑1 (MQ=255)
                                                                                 tctctcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGa                                                             <  1:33613/150‑1 (MQ=255)
                                                                                     tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                           <  2:1535585/148‑1 (MQ=255)
                                                                                     tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATc                                                           >  1:1535585/1‑148 (MQ=255)
                                                                                     tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCaa                                                         >  1:243819/1‑150 (MQ=255)
                                                                                     tcAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCaa                                                         >  2:563970/1‑150 (MQ=255)
                                                                                        ggACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGAt                                                            >  2:583660/1‑144 (MQ=255)
                                                                                        ggACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGAt                                                            <  1:583660/144‑1 (MQ=255)
                                                                                        ggACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAgct                                                      >  2:1596641/1‑150 (MQ=255)
                                                                                            ccAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGctgctg                                                  >  2:1128878/1‑150 (MQ=255)
                                                                                             cAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTgg                                                 <  2:243819/150‑1 (MQ=255)
                                                                                             cAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTgg                                                 <  1:1375719/150‑1 (MQ=255)
                                                                                               ggCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTa                                                                                                  >  2:1069212/1‑99 (MQ=255)
                                                                                               ggCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTa                                                                                                  <  1:1069212/99‑1 (MQ=255)
                                                                                                  tCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTc                                            <  1:646385/150‑1 (MQ=255)
                                                                                                   ccGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTcc                                           >  1:1468379/1‑150 (MQ=255)
                                                                                                     gCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAg                                         >  2:5110/1‑150 (MQ=255)
                                                                                                              ggtggtGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGaa                                <  1:1493654/150‑1 (MQ=255)
                                                                                                              ggtggtGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGaa                                <  1:50641/150‑1 (MQ=255)
                                                                                                                tggtgACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGttt                                                                          <  2:59984/106‑1 (MQ=255)
                                                                                                                tggtgACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGttt                                                                          >  1:59984/1‑106 (MQ=255)
                                                                                                                  gtgACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTAc                            <  1:5110/150‑1 (MQ=255)
                                                                                                                       cTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGc                       >  2:647633/1‑150 (MQ=255)
                                                                                                                        tCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGa                                                             >  2:1134492/1‑111 (MQ=255)
                                                                                                                        tCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCGGCACCGGTGCGCACCGAGTTTGGCTACCATCTGa                                                             <  1:1134492/111‑1 (MQ=255)
                                                                                                                          ggCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTc                    <  1:1128878/150‑1 (MQ=255)
                                                                                                                             tttGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGATGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGcc                 >  2:20256/1‑150 (MQ=255)
                                                                                                                                gCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGc              <  1:483259/150‑1 (MQ=255)
                                                                                                                                  cGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCt            >  1:1148960/1‑150 (MQ=255)
                                                                                                                                  aaggccGGGTCGTATACACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTgg                                                 >  1:1601724/5‑113 (MQ=255)
                                                                                                                                  aaggccGGGTCGTATACACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTgg                                                 <  2:1601724/109‑1 (MQ=255)
                                                                                                                                   aggccGGGTCGTATACACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAgg                                        <  1:1267263/118‑1 (MQ=255)
                                                                                                                                   aggccGGGTCGTATACACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAgg                                        >  2:1267263/4‑121 (MQ=255)
                                                                                                                                             gtATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACaaa  <  2:1148960/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
ACGCCCAGGCCAAGGCCCGCGCCGAAGAGATCGAGCAGCGTCTGGCCAAGGGTGAGGACTTTGCCGCTTTGGCCAAGGAGTTCTCTCAGGACCCAGGCTCCGCCAACACCGGTGGTGACCTCGGCTTTGCCGGCCCGGGTGTATACGACCCGGCCTTCGAAGAGGCGCTGTACAAGCTGCAGGATGGGCAGGTATCTGCACCGGTGCGCACCGAGTTTGGCTACCATCTGATCAAGCTGCTGGGCGTCCAGGCGCCAGAAGTACCGAGCTTCGCCAGCCTGAAGGACAAA  >  NC_002947/2634420‑2634709

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