Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913836,7560TG33.3% 69.0 / 10.9 42F136V (TTT→GTT) hcxBhydroxycarboxylate dehydrogenase B
Reads supporting (aligned to +/- strand):  ref base T (5/23);  new base G (14/0);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 2.20e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTG  >  NC_000913/836617‑836888
                                                                                                                                           |                                                                                                                                    
tGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGt                                                                                                                                      >  2:663449/1‑140 (MQ=255)
 ggTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGtt                                                                                                                                     <  1:946259/140‑1 (MQ=255)
 ggTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGtt                                                                                                                                     <  1:437323/140‑1 (MQ=255)
       gTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGtctc                                                                                                                               >  1:1028590/1‑140 (MQ=255)
             gCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCa                                                                                                                         <  1:1130536/140‑1 (MQ=255)
                     aGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTTTCCACttttttt                                                                                                                 >  2:1477094/1‑139 (MQ=255)
                     aGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACttttttt                                                                                                                 >  2:797339/1‑139 (MQ=255)
                       cGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGc                                                                                                               <  2:888298/140‑1 (MQ=255)
                            gcgcTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                                                                                          <  2:1418939/140‑1 (MQ=255)
                            gcgcTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGt                                                                                                          <  2:34391/140‑1 (MQ=255)
                            gcgcTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTTTTAGCGTggg                                                                                                          >  2:1441421/1‑139 (MQ=255)
                               cTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCgg                                                                                                       <  2:1035344/140‑1 (MQ=255)
                                 ggggATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTa                                                                                                     >  2:95329/1‑140 (MQ=255)
                                 ggggATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGGCggga                                                                                                     >  2:244565/1‑138 (MQ=255)
                                  gggATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTAt                                                                                                    >  1:313832/1‑140 (MQ=255)
                                     aTTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCCCTTTTTTAGCGTGGTCGGGATTcc                                                                                                 >  2:43623/1‑140 (MQ=255)
                                            aaGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTTTTTGCGTGGTCGGGATTCCAATGGTc                                                                                          >  2:834838/1‑140 (MQ=255)
                                            aaGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTc                                                                                          >  1:175192/1‑140 (MQ=255)
                                            aaGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTc                                                                                          >  2:444827/1‑140 (MQ=255)
                                             aGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTcg                                                                                         >  1:906537/1‑140 (MQ=255)
                                              gcgcATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTcgc                                                                                        >  1:1327392/1‑140 (MQ=255)
                                                  aTCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCCATCCACCTTGTTTGCGGGGGCGGTATTCCAATGGTCGCGCCg                                                                                    >  2:961265/1‑140 (MQ=255)
                                                   tCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGt                                                                                   <  2:614789/140‑1 (MQ=255)
                                                     aGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATTCACTTTGTTAGCGTGGTCGGTATTCCAATGGGCACGCCGTTc                                                                                 >  1:1117949/1‑140 (MQ=255)
                                                            tATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTTTCTCTATCCACTTTGTTTGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTc                                                                          >  2:1396024/1‑140 (MQ=255)
                                                             tttGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTcg                                                                         <  1:177883/139‑1 (MQ=255)
                                                              ttGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTcgc                                                                        <  1:470152/140‑1 (MQ=255)
                                                                 ccgccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGGATTCCAATGGTCGCGCCGTTCCACGGTCGCGAc                                                                     >  2:141378/1‑140 (MQ=255)
                                                                    ccgTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTAGCGCCGTTCCACGGTCGCGACAGc                                                                  >  1:1212805/1‑140 (MQ=255)
                                                                              aCATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCa                                                        <  1:797339/140‑1 (MQ=255)
                                                                              aCATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCa                                                        <  2:957834/140‑1 (MQ=255)
                                                                              aCATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCa                                                        <  1:663449/140‑1 (MQ=255)
                                                                                   aCTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAAt                                                   <  1:834838/140‑1 (MQ=255)
                                                                                      cGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCg                                                >  1:1525402/1‑140 (MQ=255)
                                                                                        catcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                              <  1:43623/140‑1 (MQ=255)
                                                                                        catcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                              <  1:1477094/140‑1 (MQ=255)
                                                                                        catcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                              <  2:1028590/140‑1 (MQ=255)
                                                                                        catcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                              <  1:940953/140‑1 (MQ=255)
                                                                                        catcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGtt                                              <  2:411960/140‑1 (MQ=255)
                                                                                         atcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTc                                             >  2:1489847/1‑140 (MQ=255)
                                                                                         atcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTc                                             >  2:1510081/1‑140 (MQ=255)
                                                                                         atcatATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTCCAATGGTCGCGCCGTTCCCCGGTCGCGACAGCCGCTTTGGCACCCATCCGTTc                                             >  2:1557067/1‑140 (MQ=255)
                                                                                            atatCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTTCAATGGTCGCGCCGTTCCACGGTCGCGACCGCCGCTTTGGCCCCAATCCGTTCtgt                                          >  2:1193940/1‑140 (MQ=255)
                                                                                                  gCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTTGGGGGGGCGGTATTCCAATGGTCGCGCCGTCCCACGGGCGCGACAGCCGCTTTGGGGACAATCCCGTCTGTGTgggt                                    >  2:1160245/1‑138 (MQ=255)
                                                                                                     gTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTc                                 >  1:1340525/1‑140 (MQ=255)
                                                                                                      tATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTcc                                <  1:95329/140‑1 (MQ=255)
                                                                                                      tATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTAGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTcc                                <  2:1212805/140‑1 (MQ=255)
                                                                                                       aTCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTTGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTccc                               >  2:1454476/1‑140 (MQ=255)
                                                                                                       aTCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGGCGGGATTCCAATTGTCGCGCCGTTCCACGGTCGCGGCAGCCGCCTTGGCACCAATCCGTTCTGTGTGGGTTTtcc                               >  2:323943/1‑140 (MQ=255)
                                                                                                       aTCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACCTTGTTAGCGTGGTCGGGATTCCCATAGTCGCGCCGTTTCCCGGTCGCGACAGCCGCTTTGGCACCAATCCGGTTTGTGTGGGTTTtcc                               >  2:1287887/1‑140 (MQ=255)
                                                                                                        tCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGGTAGCGCGGGCGGGATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGGACCAATCCCTTCTGTGTGGGTTTTCCt                              >  2:872755/1‑140 (MQ=255)
                                                                                                                tGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCAAAGGGTCGCGCCGTTCCACGGTCGCGACAGCCGTTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGa                      <  2:313832/140‑1 (MQ=255)
                                                                                                                tGGGCGGAGCAGTGTGCAGCGGCGGGGGTTTTCTCTATCCCCTTTGTTTGCGTGGTCGGTATTCCCATGGTCGCGCCGTTCCCCGGTCGCGACAGCCGCTTTGGCACCAATCCGGTCTGTGTGGGTTTCCCTCGTAAAGa                      >  2:34113/1‑140 (MQ=255)
                                                                                                                  ggCGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATa                    >  1:145026/1‑140 (MQ=255)
                                                                                                                    cGGAGCAGTGTGCAGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGGATTCCAATGGGCGCGCCGGTCCACGGGCGCGACCGCCGCCTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGCAAAaataat                  >  2:682135/1‑140 (MQ=255)
                                                                                                                           gtgtGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGc           <  1:1307265/140‑1 (MQ=255)
                                                                                                                           gtgtGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGc           <  1:108973/140‑1 (MQ=255)
                                                                                                                              tgCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGt        >  1:287657/1‑140 (MQ=255)
                                                                                                                              tgCAGCGGCGGGGGTTTTCTCTTTCCACCTTGTTAGCGTGGTCGGGATTCCAATGGGCGCGCCGTTCCCCCGGCGCGACAGCCGCTTTGGCACCAATCCCCTCTGTGTGGGTTTCCCCCGCAAAGATAATTTCCCGCTGt        >  2:246070/1‑140 (MQ=255)
                                                                                                                                 aGCGGCGGGGGTTGTCTCTATCCACTTTGTTAGCGTGGTCGGGATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGc     >  1:686084/1‑140 (MQ=255)
                                                                                                                                    ggcggGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTg  <  2:1327392/140‑1 (MQ=255)
                                                                                                                                    ggcggGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTg  <  1:1430486/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                    
TGGTCAGGTCGCGGCACATGAAGCGATGGCGCTGGGGATTGAGAAAGCGCATCAGCACGGTATTGCCGCCGTGGCGCTACATAACTCGCATCATATCGGCCGTATCGGTTACTGGGCGGAGCAGTGTGCAGCGGCGGGGTTTGTCTCTATCCACTTTGTTAGCGTGGTCGGTATTCCAATGGTCGCGCCGTTCCACGGTCGCGACAGCCGCTTTGGCACCAATCCGTTCTGTGTGGTTTTCCCTCGTAAAGATAATTTCCCGCTGTTGCTTG  >  NC_000913/836617‑836888

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: