Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,387,274 T→C 6.0% I43T (ATC→ACC)  PP_1207 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,387,2740TC 6.0% 192.1 / 2.4 66I43T (ATC→ACC) PP_1207hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/34);  new base C (2/2);  total (30/36)
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.09e-01

GCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTG  >  NC_002947/1387129‑1387413
                                                                                                                                                 |                                                                                                                                           
gCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCa                                                                                                                                           <  2:24054/148‑1 (MQ=255)
 ccTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGGTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCTATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGt                                                                                                                                         <  5:51422/149‑1 (MQ=255)
 ccTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTCTGCAAGTTGCTGATGAGCAAGGATGCCAATTACCCATGGTTCATTATGGTGCCCAAGCGCGCTGGTATCAGt                                                                                                                                         >  7:62162/1‑149 (MQ=255)
  cTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAgcg                                                                                                                                        >  3:15839/1‑147 (MQ=255)
       ggAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTg                                                                                                                                   >  4:348433/1‑149 (MQ=255)
         agaTATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgg                                                                                                                                 <  4:15839/149‑4 (MQ=255)
               gCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCt                                                                                                                           >  1:197829/1‑149 (MQ=255)
                ccGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCt                                                                                                                           <  8:205970/148‑1 (MQ=255)
                        gTCCTCGATTCGCGGTTGCTGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   <  3:23952/148‑1 (MQ=255)
                         tCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGACCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCa                                                                                                                 >  7:338169/1‑149 (MQ=255)
                           cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGa                                                                                                                              <  5:143751/134‑1 (MQ=255)
                           cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGa                                                                                                                              >  6:143751/1‑134 (MQ=255)
                           cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   <  3:196384/145‑1 (MQ=255)
                           cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCgg                                                                                                                   >  4:196384/1‑145 (MQ=255)
                           cTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCaga                                                                                                               <  1:227918/149‑1 (MQ=255)
                              gATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGc                                                                                                             >  2:114724/1‑148 (MQ=255)
                                      tttGCAGCAGGATTCGCGGGTGTTGGGAGAGTTTCCGGTGTGCCAGTTGCTGGTGAGCAAGGATGCCAATTACCCATGCTTCATTCTGGTGCCGAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGt                                                                                                    <  5:109604/149‑1 (MQ=255)
                                         gcagcaGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTtg                                                                                                  <  3:348433/148‑1 (MQ=255)
                                          cagcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgggc                                                                                                                               <  6:359261/118‑6 (MQ=255)
                                          cagcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCTCACTGATACCAGCGCGCttgggc                                                                                                                               >  5:359261/1‑113 (MQ=255)
                                            gcagGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTgg                                                                                               >  7:134017/1‑148 (MQ=255)
                                              agGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAg                                                                                            >  8:54694/1‑149 (MQ=255)
                                                gATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCag                                                                                                                <  8:275952/127‑1 (MQ=255)
                                                gATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCag                                                                                                                >  7:275952/1‑127 (MQ=255)
                                                    cGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAac                                                                                       <  4:288226/148‑1 (MQ=255)
                                                     gCTGGTGTTGGTAGATTTTCCGCTGTGGCAGTTGCTGCTGAGCAAGGATGCCTATTTCCCTTGGTTCTTTCTTGTGGCCAAGCGCTGTTGTATCAGTGAGCTGTTCTATCTGGTTTCGGGAGAGCAACTTCAGTTGTGGCTGGATacca                                                                                     <  8:62162/149‑1 (MQ=255)
                                                     gCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAacca                                                                                     <  8:338169/149‑1 (MQ=255)
                                                                         cGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGtcg                                                                 <  8:188455/149‑3 (MQ=255)
                                                                         cGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTg                                                                 <  4:79017/149‑1 (MQ=255)
                                                                                ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGaaa                                                              >  2:249234/1‑145 (MQ=255)
                                                                                ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGaaa                                                              <  1:249234/145‑1 (MQ=255)
                                                                                ccAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           <  8:134017/148‑1 (MQ=255)
                                                                                       ctgctgAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATgg                                                   <  6:103121/149‑1 (MQ=255)
                                                                                            gAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGGGCTGGTATCAGTGAGCTGTTCGAGCTGGAt                                                                                                                       <  1:31685/76‑1 (MQ=255)
                                                                                            gAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGAt                                                                                                                       >  2:31685/1‑76 (MQ=255)
                                                                                              gCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGAc                                            <  2:197829/149‑1 (MQ=255)
                                                                                                  ggATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGAc                                            >  4:60891/1‑145 (MQ=255)
                                                                                                       ccAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAAc                                   <  3:118757/149‑1 (MQ=255)
                                                                                                         aaTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGt                                 >  4:39233/1‑149 (MQ=255)
                                                                                                         aaTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGt                                 >  8:12369/1‑149 (MQ=255)
                                                                                                          aTTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTGGAAGGCCAGCTATGGGGGCGACAAGATGAACGTg                                >  8:322083/1‑149 (MQ=255)
                                                                                                           ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTAt                                                     >  5:300203/1‑127 (MQ=255)
                                                                                                           ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTAt                                                     <  6:300203/127‑1 (MQ=255)
                                                                                                           ttACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTgg                               <  7:12369/149‑1 (MQ=255)
                                                                                                               ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGTAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           >  4:319169/1‑117 (MQ=255)
                                                                                                               ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGTAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGcc                                                           <  3:319169/117‑1 (MQ=255)
                                                                                                               ccATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCACCTATGGGGCCGACAAGATGAACGTGGCTa                            >  7:239580/1‑148 (MQ=255)
                                                                                                                 aTGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGACCGTGGCGACCt                         >  8:235326/1‑149 (MQ=255)
                                                                                                                   ggTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGa                                        >  7:66358/1‑132 (MQ=255)
                                                                                                                   ggTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGa                                        <  8:66358/132‑1 (MQ=255)
                                                                                                                      tCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTggg                     <  7:658/148‑1 (MQ=255)
                                                                                                                      tCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTggg                     >  8:658/1‑148 (MQ=255)
                                                                                                                         ttCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACaa                                          <  6:39952/124‑1 (MQ=255)
                                                                                                                         ttCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACaa                                          >  5:39952/1‑124 (MQ=255)
                                                                                                                            tGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTg              >  3:149278/1‑149 (MQ=255)
                                                                                                                              gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCATAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTg              <  4:324502/147‑1 (MQ=255)
                                                                                                                              gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCTACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            >  8:70360/1‑149 (MQ=255)
                                                                                                                              gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGACCGTGGCTACCTTGGGCAACGTg              <  7:54694/147‑1 (MQ=255)
                                                                                                                              gTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            <  4:149278/149‑1 (MQ=255)
                                                                                                                               tGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTc           <  3:60891/149‑1 (MQ=255)
                                                                                                                                gCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGt            >  3:155888/1‑147 (MQ=255)
                                                                                                                                  ccAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGTCAGGAAACGTCCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGc        <  1:114724/149‑1 (MQ=255)
                                                                                                                                  ccAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGATGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGc        >  6:105353/1‑149 (MQ=255)
                                                                                                                                        gcgcTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTg  <  4:155888/149‑1 (MQ=255)
                                                                                                                                              gTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTtgt                                                                                                 <  4:55497/48‑1 (MQ=255)
                                                                                                                                              gTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTtgt                                                                                                 >  3:55497/1‑48 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GCCTGAAGGAGATATGCCGTGTTCGTCCTCGATTCGCGTTTGCAGCAGGATTCGCTGGTGTTGGGAGAGTTTCCGCTGTGCCAGTTGCTGCTGAGCAAGGATGCCAATTACCCATGGTTCATTCTGGTGCCCAAGCGCGCTGGTATCAGTGAGCTGTTCGAGCTGGATGCGGCAGAGCAACAGCAGTTGTGGCAGGAAACCACCTTGTTGGCAGAGGCGTTGAAAGCCAGCTATGGGGCCGACAAGATGAACGTGGCTACCTTGGGCAACGTGGTCAGCCAGTTG  >  NC_002947/1387129‑1387413

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