Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 853,763 T→C 6.5% R62R (CGT→CGC prmC → release factor‑(glutamine‑N5) methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947853,7630TC 6.5% 178.1 / 2.3 64R62R (CGT→CGCprmCrelease factor‑(glutamine‑N5) methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/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 = 5.66e-01

GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGCCCGCCAAG  >  NC_002947/853616‑853910
                                                                                                                                                   |                                                                                                                                                    
gCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGt                                                                                                                                                      <  4:47016/148‑1 (MQ=255)
  cTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTggc                                                                                                                                                  >  1:221799/1‑149 (MQ=255)
  cTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTggc                                                                                                                                                  >  4:148561/1‑149 (MQ=255)
       tCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGa                                                                                                                                             <  7:261185/149‑1 (MQ=255)
          cccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGc                                                                                                                                           >  1:319544/1‑148 (MQ=255)
          cccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttc                                                                                                                                          <  1:82997/149‑5 (MQ=255)
           ccACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCACGgcggcgcttca                                                                                                                                         <  2:50842/149‑6 (MQ=255)
            cACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAg                                                                                                                                        >  7:152609/1‑149 (MQ=255)
            cACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAg                                                                                                                                        >  7:148216/1‑149 (MQ=255)
             aCTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGt                                                                                                                                       <  4:154645/149‑1 (MQ=255)
               tGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGGCTCGCCGCCACGgcggcgcttcagg                                                                                                                                       >  8:240091/1‑140 (MQ=255)
                 aGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGc                                                                                                                                    <  8:152609/148‑1 (MQ=255)
                 aGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGc                                                                                                                                    <  8:148216/148‑1 (MQ=255)
                     ccTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCtata                                                                                                                               >  4:44403/1‑149 (MQ=255)
                              cGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGc                                                                                                                      >  6:338575/1‑149 (MQ=255)
                                     ctgctgGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCGTTGGGCTGCAgg                                                                                                                <  3:44403/148‑1 (MQ=255)
                                     ctgctgGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAg                                                                                                                 <  7:275355/147‑1 (MQ=255)
                                        ctggctgTGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGCCCCGAGCGTATTGTCAGCAGCGAACATGCCGAGACCTATGCCGGCAACCTGAAGCCCCGCCCTGCCGGCGAGCCACTGGCCTATATCCTTCGGCTGCAGGCCtg                                                                                                            >  2:153967/1‑148 (MQ=255)
                                            ctgCGTCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                        <  2:319544/149‑1 (MQ=255)
                                            ctgCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                        <  2:306592/149‑1 (MQ=255)
                                             tgCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                        >  8:128818/1‑148 (MQ=255)
                                                  ccATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTgg                                                                                                        >  1:246828/1‑143 (MQ=255)
                                                     tCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                               >  3:102050/1‑149 (MQ=255)
                                                     tCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                               <  2:221799/149‑1 (MQ=255)
                                                       ggTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCt                                                                                             <  4:102050/149‑1 (MQ=255)
                                                        gTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                               >  7:169056/1‑146 (MQ=255)
                                                        gTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGAc                                                                                               <  8:169056/146‑1 (MQ=255)
                                                        gTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTg                                                                                            >  6:110173/1‑149 (MQ=255)
                                                           aaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCTGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAg                                                                                         <  7:128818/149‑1 (MQ=255)
                                                           aaTCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCCTGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGa                                                                                          <  1:95216/149‑1 (MQ=255)
                                                             tCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGt                                                                                       >  5:78364/1‑149 (MQ=255)
                                                               gcgcAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                                      >  6:92888/1‑132 (MQ=255)
                                                               gcgcAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                                      <  5:92888/132‑1 (MQ=255)
                                                                 gcAGCTACCTGCATACCTGGCCCGAGCGTTTTGTCTGCAGCGAAGATGCCGAGACGTATGCCGGCTACCTGAAGCGCCGCCGGGGCGGCGAGCCAGTGGCCTATATGCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTcgcg                                                                                   <  1:153967/149‑1 (MQ=255)
                                                                       aCCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCa                                                                              >  3:154492/1‑148 (MQ=255)
                                                                        ccTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCTGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCa                                                                              <  5:338575/147‑1 (MQ=255)
                                                                         cTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGa                                                                                                    >  6:159815/1‑124 (MQ=255)
                                                                         cTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGa                                                                                                    <  5:159815/124‑1 (MQ=255)
                                                                          tGCATACCTGGCCCGAGCGTATTGTCGGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTAGATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATAcc                                                                          >  4:9429/1‑149 (MQ=255)
                                                                                    gCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                                      <  3:259077/111‑1 (MQ=255)
                                                                                    gCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGaa                                                                                                      >  4:259077/1‑111 (MQ=255)
                                                                                       cGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGc                                                             <  4:154492/149‑1 (MQ=255)
                                                                                        gAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGcc                                                            >  6:342295/1‑149 (MQ=255)
                                                                                               ttGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATAc                                                      <  2:246828/148‑1 (MQ=255)
                                                                                                 gTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCg                                                    <  6:78364/148‑1 (MQ=255)
                                                                                                    agcagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATccc                                                                  <  4:44960/131‑1 (MQ=255)
                                                                                                    agcagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATccc                                                                  >  3:44960/1‑131 (MQ=255)
                                                                                                    agcagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAgct                                                >  3:178372/1‑149 (MQ=255)
                                                                                                    agcagcGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAgct                                                <  5:110173/149‑1 (MQ=255)
                                                                                                           aaGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATACCGCGGCCGGATACCGAGCTGCTGGTc                                         >  3:247144/1‑149 (MQ=255)
                                                                                                            cgctGCCGAGCCCTCGGCCGGCTACCTGACGCCCCGCCGGGTCGGCGAGCCAGTGGCCTATCTCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCg                                        <  4:247144/146‑1 (MQ=255)
                                                                                                                        cTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCtt                             >  8:283009/1‑148 (MQ=255)
                                                                                                                           tGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGaa                          <  5:342295/148‑1 (MQ=255)
                                                                                                                             ccGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACt                        <  3:215607/148‑1 (MQ=255)
                                                                                                                              cGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGc                      <  4:178372/149‑1 (MQ=255)
                                                                                                                                gCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCac                    >  2:12088/1‑148 (MQ=255)
                                                                                                                                     cTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTc                                         <  1:63748/123‑1 (MQ=255)
                                                                                                                                     cTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTc                                         >  2:63748/1‑123 (MQ=255)
                                                                                                                                      tGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTgg                                                                                           >  2:107600/1‑72 (MQ=255)
                                                                                                                                      tGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGctgctg                                            >  7:306793/1‑119 (MQ=255)
                                                                                                                                      tGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGctgctg                                            <  8:306793/119‑1 (MQ=255)
                                                                                                                                        gctcgccgcCACGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCg                                                           >  6:229468/4‑102 (MQ=255)
                                                                                                                                        aaGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCt              <  7:283009/147‑1 (MQ=255)
                                                                                                                                            gccgccGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTcgccc        >  6:261672/1‑149 (MQ=255)
                                                                                                                                                  gTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGCCCGCCAAg  >  2:163571/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCTGAATCGCCCACTGAGCGCCTGGATGCCGAACTGCTGCTGGCTGCGGCCATCGGTAAATCGCGCAGCTACCTGCATACCTGGCCCGAGCGTATTGTCAGCAGCGAAGATGCCGAGACCTATGCCGGCTACCTGAAGCGCCGCCGTGGCGGCGAGCCAGTGGCCTATATCCTTGGGCTGCAGGGCTTCTGGAAGATCGACCTGGAGGTCGCGCCGCATACCCTGATCCCGCGGCCGGATACCGAGCTGCTGGTCGAGACGGCGCTTGAACTGCAACCCGCTTCGCCCGCCAAG  >  NC_002947/853616‑853910

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: