Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,633,135 C→T 8.3% L203L (CTG→TTG)  PP_3201 → BNR domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,633,1350CT 8.3% 172.7 / 4.2 60L203L (CTG→TTG) PP_3201BNR domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (24/31);  new base T (3/2);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

GCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTC  >  NC_002947/3632987‑3633279
                                                                                                                                                    |                                                                                                                                                
gCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGccc                                                                                                                                                  <  1:115430/149‑1 (MQ=255)
gCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGccc                                                                                                                                                  >  5:76188/1‑149 (MQ=255)
   gttgttGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTgg                                                                                                                                               <  6:164710/149‑1 (MQ=255)
   gttgttCGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGCTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGAACGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTgg                                                                                                                                               >  8:123129/1‑149 (MQ=255)
    ttgttgGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTATATGCCCTGgc                                                                                                                                              >  3:210464/1‑149 (MQ=255)
       ttgGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGgcgcg                                                                                                                                           >  1:95017/1‑149 (MQ=255)
         ggCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGgcgcgc                                                                                                                                          <  8:81702/148‑1 (MQ=255)
          gCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGcc                                                                                                                                         <  2:228192/148‑1 (MQ=255)
               gAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAgg                                                                                                                                     >  1:194786/1‑147 (MQ=255)
                aaTTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                                                                                                                                    <  3:371196/147‑1 (MQ=255)
                  tttCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAAc                                                                                                                                 >  6:24868/1‑148 (MQ=255)
                      aCCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGgt                                                                                                                            >  1:113320/1‑149 (MQ=255)
                        cgctGCCCGAACGGGCCGGGTCGTCGGTGCCTAGGGGGGGGCAATGCCCAGCTAGGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAg                                                                                                                                      <  1:381347/134‑1 (MQ=255)
                         gATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtgca                                                                                                                         <  4:104361/149‑3 (MQ=38)
                          ctGCCCGCCACGGCCTGGTCGTCGGTGCCTACGGGGGGGCCATGCGCACCGACGACGGCGTGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTAc                                                                                                                        <  6:77211/148‑1 (MQ=255)
                              ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCg                                                                                                                    >  3:297685/1‑149 (MQ=255)
                              ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtg                                                                                                                           >  1:8503/1‑142 (MQ=39)
                              ccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGGTTGCCTTGGCGCGCCAGGGCATACAGgtg                                                                                                                           <  2:8503/142‑1 (MQ=39)
                               cGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGc                                                                                                                   >  6:40846/1‑149 (MQ=255)
                                   aCGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGc                                                                                                               <  5:24868/149‑1 (MQ=255)
                                        cTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAg                                                                                                            >  3:353941/1‑147 (MQ=255)
                                        cTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAg                                                                                                            >  8:153774/1‑147 (MQ=255)
                                           gtcgtcGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAgg                                                                                                        <  5:322281/148‑1 (MQ=255)
                                            tcgtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                      <  1:368240/149‑1 (MQ=255)
                                            tcgtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGc                                                                                                      <  6:191098/149‑1 (MQ=255)
                                              gtcgGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCct                                                                                                    >  6:62993/1‑149 (MQ=255)
                                                      cTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTATATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtc                                                                                            <  4:210464/149‑1 (MQ=255)
                                                       tACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTtcg                                                                                           >  5:184319/1‑149 (MQ=255)
                                                              tGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                    <  7:123129/149‑1 (MQ=255)
                                                               ggCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGt                                                                                    >  2:328932/1‑148 (MQ=255)
                                                                  cATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                <  5:40846/149‑1 (MQ=255)
                                                                  cATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATg                                                                                <  5:311934/149‑1 (MQ=255)
                                                                   aTGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATgg                                                                               >  3:138630/1‑149 (MQ=255)
                                                                      cgcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCg                                                                             <  4:353941/148‑1 (MQ=255)
                                                                       gcACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGt                                                                           >  7:201289/1‑149 (MQ=255)
                                                                                        cgcAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCt                                                          <  2:95017/149‑1 (MQ=255)
                                                                                             cTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACg                                                     >  4:371487/1‑149 (MQ=255)
                                                                                              tGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACgg                                                    >  2:250264/1‑149 (MQ=255)
                                                                                               ggCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACggg                                                   >  2:169718/1‑149 (MQ=255)
                                                                                                gCAGTCGGAGCTCGGCCGACGGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGc                                                  <  5:265918/149‑1 (MQ=255)
                                                                                                gCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACggg                                                   <  8:201289/148‑1 (MQ=255)
                                                                                                   gTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCt                                               <  4:297685/149‑1 (MQ=255)
                                                                                                      ggAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTAc                                             <  5:212438/148‑1 (MQ=255)
                                                                                                       gagaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGa                                           <  6:184319/149‑1 (MQ=255)
                                                                                                        agaTGGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCGTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                          <  6:82689/149‑1 (MQ=255)
                                                                                                        agaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                          <  3:203129/149‑1 (MQ=255)
                                                                                                        agaTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGAt                                          <  3:371487/149‑1 (MQ=255)
                                                                                                          aTCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATgg                                        >  7:270368/1‑149 (MQ=255)
                                                                                                             ggCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCa                                      <  3:22108/148‑1 (MQ=255)
                                                                                                               ccGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGc                                    >  4:94825/1‑148 (MQ=255)
                                                                                                                         aaCGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTcgccgc                         >  3:256030/1‑149 (MQ=255)
                                                                                                                                   gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGtt                <  3:94825/148‑1 (MQ=255)
                                                                                                                                   gACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGtt                <  1:250264/148‑1 (MQ=255)
                                                                                                                                    aCTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGc              >  6:57189/1‑149 (MQ=255)
                                                                                                                                     cTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGcc             <  2:113320/149‑1 (MQ=255)
                                                                                                                                      tGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGccc            >  2:218380/1‑149 (MQ=255)
                                                                                                                                          ccTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATg        <  7:153774/149‑1 (MQ=255)
                                                                                                                                           ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATgg       >  1:229005/6‑149 (MQ=255)
                                                                                                                                           ctggttGCCTTGGCGCGCCAGGGCATACAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCAGATgg       >  6:211853/6‑149 (MQ=255)
                                                                                                                                              tATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCg    <  4:138630/149‑1 (MQ=255)
                                                                                                                                                tGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTc  >  2:266212/1‑149 (MQ=255)
                                                                                                                                                tGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTc  <  1:218380/149‑1 (MQ=255)
                                                                                                                                                   ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCAt                                                                     >  2:54251/1‑79 (MQ=255)
                                                                                                                                                   ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCAt                                                                     <  1:54251/79‑1 (MQ=255)
                                                                                                                                                   ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCACTGCACGGGCCCTACGATGGCAGCTACTTcgccg                          <  7:226955/122‑1 (MQ=255)
                                                                                                                                                   ccTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCACTGCACGGGCCCTACGATGGCAGCTACTTcgccg                          >  8:226955/1‑122 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGCCGGCGAGCAGGGCCTGCTGCTTCGCTCGCGTGATGGCGGTGAGCATTACGAGGCGCTGCACGGGCCCTACGATGGCAGCTACTTCGCCGCCACCGTGTTGCCCGATGGCCGTC  >  NC_002947/3632987‑3633279

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: