Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,206,4350CT50.0% 36.9 / 47.8 58intergenic (‑179/‑28)ttdR/ttdAtranscriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (3/26);  new base T (29/0);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 3.30e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAA  >  NC_000913/3206287‑3206583
                                                                                                                                                     |                                                                                                                                                     
gtgtATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCa                                                                                                                                                     <  1:664350/151‑1 (MQ=255)
    aTTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTCTAGa                                                                                                                                                   >  2:887497/1‑149 (MQ=255)
         gCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAc                                                                                                                                            <  1:432883/151‑1 (MQ=255)
           aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTcaaacaaaaa                                                                                                                                             >  2:735696/1‑146 (MQ=255)
           aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATACa                                                                                                                                           >  2:285185/1‑150 (MQ=255)
            aaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACaaa                                                                                                                                         <  1:487103/151‑1 (MQ=255)
              aCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACaaaa                                                                                                                                        >  2:816014/1‑150 (MQ=255)
                aTTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACt                                                                                                                                     <  1:697415/151‑1 (MQ=255)
                     cATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGt                                                                                                                               <  2:398888/151‑1 (MQ=255)
                         cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGcc                                                                                                                           <  1:174821/151‑1 (MQ=255)
                           cTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCCGACCAAAACAAAAAAAGGAGGTGCCa                                                                                                                          >  2:865729/1‑150 (MQ=255)
                                ctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGGCatg                                                                                                                        >  2:408212/1‑147 (MQ=255)
                                 tGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCTGAAAAATAAAAAAAATGGAGTTGCCatgatgtg                                                                                                                   >  1:879333/1‑149 (MQ=255)
                                   gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                                 <  1:345909/151‑1 (MQ=255)
                                   gcgcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCg                                                                                                                 <  2:463123/151‑1 (MQ=255)
                                     gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAAAACAAAAACTGGAGTTTGCatggtggggg                                                                                                                 >  1:201924/1‑145 (MQ=255)
                                     gcAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGaa                                                                                                               <  2:567657/151‑1 (MQ=255)
                                             cattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGtaataa                                                                                                       <  2:495739/151‑1 (MQ=255)
                                                tcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGc                                                                                                     <  2:201953/150‑1 (MQ=255)
                                                   tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAAATTGAGGTGCCATGATGAGCGAAAGTAATAAACaa                                                                                                   >  2:179597/1‑149 (MQ=255)
                                                    aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACa                                                                                                 <  2:499105/150‑1 (MQ=255)
                                                    aTTTTCGCGTTGATTCACGAACTATTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTAGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAg                                                                                                <  2:597233/151‑1 (MQ=255)
                                                        tCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCa                                                                                             <  1:643543/150‑1 (MQ=255)
                                                            gTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGaa                                                                                        <  2:425675/151‑1 (MQ=255)
                                                                aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAAATGGAGTTGCCATGATGAGGGAAAGGAATAAACAAAAGGGAGTGAATaa                                                                                     >  2:780374/1‑150 (MQ=255)
                                                                aTTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAAACAAAAAAAAAAAAAGGAGGTTCCATGATGAGGGAAAGGAAAAAAAAAAAGGCAGTGAaaaa                                                                                     >  1:202426/1‑150 (MQ=255)
                                                                   cACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTg                                                                                 <  1:132624/151‑1 (MQ=255)
                                                                     cGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGAc                                                                               <  1:159456/151‑1 (MQ=255)
                                                                      gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAAATGGAGGTGCCATGATTAGGGAAAGTAATAAAAAACAAGCAGGGGATAAGGTGAc                                                                               >  2:99539/1‑150 (MQ=255)
                                                                       aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTCTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCCTGATGAGCGGAAGTAATAAACAACAAGCAGGGGATAAAtt                                                                                  >  2:393676/1‑146 (MQ=255)
                                                                             tAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAtt                                                                        <  2:310439/150‑1 (MQ=255)
                                                                                  aCATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAACTTGAGTTGCCATGATTAGCGGAAGTAATAAACAACAGGCAGTGAATAAGGTGACAGAGAATGTCGc                                                                   >  2:225081/1‑150 (MQ=255)
                                                                                           cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                                           <  1:874920/149‑1 (MQ=255)
                                                                                           cGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTAcc                                                         <  2:882181/151‑1 (MQ=255)
                                                                                              tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGAGGGAAAGAAATAAACAAAAAGCAGGGAATAAGTTGACAGAGAATGTCGCTAACTTTACCGcc                                                      >  2:481130/1‑151 (MQ=255)
                                                                                               ggCAAAAAATATAAAACCACCTTTTTAGTGGTAGTTTGTGGAGGCCAATTTTTTCAAACAAATACAAAAAATGGCAGTGGTAATATTAGCGACAGGAATAAAAAACAAGGAATGGATAAGAAGAAAGAGACTGTAGATaaacttaattctc                                                     >  2:825225/1‑140 (MQ=255)
                                                                                                gCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAACAAAAAATGGAGGTGCCATGATGAGGGGAAGTAATAAACAACAAGCAGGGAAAAAAATGACAGAGATTGTCGCTAAATTTaacccc                                                      >  1:75366/1‑144 (MQ=255)
                                                                                                   aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGGGGTGGCATGAGGGGGGGAAGAAATAAACAACAGGCAGTGAATAAAATGAAAAAGATTGTCGCTAACATTACAGCCATg                                                   >  1:317702/1‑149 (MQ=255)
                                                                                                   aaaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATAAAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAACCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAt                                                 <  1:34725/151‑1 (MQ=255)
                                                                                                    aaaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGGTGAGCGGAAGGAATAAACAACAAGCAGTGAATAAGTTGACAGAGATTTTTGCTAACTTTACCGCCATGAtt                                                >  2:806749/1‑151 (MQ=255)
                                                                                                        tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTTCAAAAAAAAAAAAAAATGGGGGTGGCATGATGGGGGAAAAAAAAAAACAAAAAGCAGGGAAAAAATTGGCAGAGAATTTTGCTTAATTTTTCGCCCTGAttt                                               >  1:416360/1‑148 (MQ=255)
                                                                                                        tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCAGACAAAAACAAAAAAATGGGTTGCCCTGATGAGCGGAAAAAATAAACAACAAGCAGAGAAAAAGTTTACAGAGAATGTTGCTTACTTTACCCCCATGaataaa                                             >  2:575609/1‑145 (MQ=255)
                                                                                                         ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTTCCAACCAATTCAAAAAATGGGGTTTCCTTGATGAGCGAAAGGAATAAAGCACAGGGAGGGAATTAATTGGAACGGATTGTTGCTAAATTTTcccgcccgattttt                                             >  1:118229/1‑137 (MQ=255)
                                                                                                         ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGGAAGGAATAAACAAAAGGCAGTGAAAAAATTGACAGAGAATGTCGCTAAATTTAACGCCATGATTTCt                                             >  2:216005/1‑149 (MQ=255)
                                                                                                         ataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTAAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGAtt                                                <  1:575567/146‑1 (MQ=255)
                                                                                                            aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAACAAAAACAAAAAATTGAGTTGGCATGGTGGGGGGAAGTAAAAAACAAAAAGGAGGGAAAAAAGTGACAAAGATTGTCGCTTAATTTAACGCCCTGAttttttacc                                         >  1:11331/1‑144 (MQ=255)
                                                                                                            aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGTTGCCATGATGGGGGGAAGGAATAAGGAACAGGGAGTGGATAAGATGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTAccc                                         >  2:187663/1‑150 (MQ=255)
                                                                                                            aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGAGGTGCCATGATGAGCGAAAGGAATAAACAAAAAGCAGTGAATAAATTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTTTTccc                                         >  2:780153/1‑150 (MQ=255)
                                                                                                            aaaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAGAAAACAAAAAAGTGAGTGGTGCTGATGTGGGAAAGGAATAACCAAAAAGCGGTGAATAAGGTGAGAGAGAGTTTCGCTAAAATTAGCGCCATGAtttattcacc                                        >  2:494247/1‑144 (MQ=255)
                                                                                                             aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTCCGACAAAAAAAAAAAATGGAGTTGCCCTTATGAGGGGAAGTAATAAGCAAAAGGCAGTGAAAAAGGTTAAAGAGATTGTCGCTAACTTTACCGCCCTGGTTTTTAacc                                         >  2:801671/1‑148 (MQ=255)
                                                                                                              aCCACATTTTTAGTGGTAGTTTGTGGCGGGGAATTTTTTCATACAAAAACAAAAACTGGAGGTTGCATGATGGGCGGAAGGAAAAAACAAAAAGCAGTGAATAAAATGACAGAGAATTTCGCTAAATTTTTCCCCCTGATTTTTatc                                          >  1:568401/1‑145 (MQ=255)
                                                                                                               ccACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCAAGAAAAATAAAAAAACTGGAGTTGCCATGATGAGCAAAATTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAAt                                     <  2:276822/151‑1 (MQ=255)
                                                                                                                cacaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGGCAAAAACAAAAAATGGAGTTGGCCTGGTGGGGGAAAGTAATAAGCAACAAGCAGTGAATAAATTTACAGAGATTTTCGCTAACTTTTACGCCCTTATTTTTAcccccatg                                    >  1:471570/1‑146 (MQ=255)
                                                                                                                 acaTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAAATGGGGTTGCCATGATGAGCGAAAGGAATAAACAACAGGCAGGGAATAAATTGACAGAGATTGTTGCTAACTTTACCGCCATGATTTTTACCCGaa                                      >  2:84166/1‑148 (MQ=255)
                                                                                                                  caTTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGCCGCTAACTTTACCGCCATGATTTCTACCCGAATGc                                   <  2:394080/150‑1 (MQ=255)
                                                                                                                      tttAGTGGTAGTTTGTGGCGGGGAATTTTTTCCGAAAAAAAAAAAAACTGGGGTTGCCATGATGAGGGAAAGTAATAAGCAACAAGGAGTGAATAAATTGTAAGAGAATGTCGGTAAATTTACCCCCCTGAttttttacccaa                                      >  2:216573/1‑134 (MQ=255)
                                                                                                                            ggTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACgtg                        <  2:485576/151‑1 (MQ=255)
                                                                                                                                                   ttCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTa   <  2:207599/150‑1 (MQ=255)
                                                                                                                                                     ccAGACAAATACAAAAACGGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGTTTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTaa  <  2:885837/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGTATTTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGAGTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAA  >  NC_000913/3206287‑3206583

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: