Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,7570TG42.6% 114.6 / 13.1 61G16G (GGT→GGGispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (1/34);  new base G (26/0);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 2.22e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333616‑3333906
                                                                                                                                              |                                                                                                                                                          
tttttCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGc                                                                                                                                                  <  1:383157/150‑1 (MQ=255)
     cAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCaa                                                                                                                                             <  2:608905/150‑1 (MQ=255)
       gTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGCGGTTTATACGGCAAtac                                                                                                                                          >  1:949351/1‑149 (MQ=255)
         aCAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCtt                                                                                                                                        <  2:97771/151‑1 (MQ=255)
         aCAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGCggga                                                                                                                                              >  1:549335/1‑143 (MQ=255)
          cAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTg                                                                                                                                       <  2:887319/151‑1 (MQ=255)
           aaTCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTAATgggg                                                                                                                                                >  2:130806/1‑141 (MQ=255)
              cACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGagc                                                                                                                                    <  1:826414/150‑1 (MQ=255)
                cccTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTAATTAGGCAATCCTTGagc                                                                                                                                    >  2:354128/1‑148 (MQ=255)
                    atatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCtt                                                                                                                             <  1:26262/149‑1 (MQ=255)
                    atatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCtt                                                                                                                             <  1:391749/149‑1 (MQ=255)
                     tatTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTGATTCGGCCAACCTTTagg                                                                                                                                    >  2:610459/1‑142 (MQ=255)
                       ttCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATACGGGAAACCTTTAGCAGCTTAAtt                                                                                                                         >  2:37150/1‑150 (MQ=255)
                         ccTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTcc                                                                                                                       <  1:809914/150‑1 (MQ=255)
                         ccTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATTAGGGAATCCTTTTGCAGCCTAAtttc                                                                                                                       >  2:1011616/1‑148 (MQ=255)
                          cTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATACGGGAAACCTTTAGGAGCGTTAATTcgc                                                                                                                      >  2:903939/1‑149 (MQ=255)
                            aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGa                                                                                                                     <  1:479773/149‑1 (MQ=255)
                            aaCCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACg                                                                                                                   <  2:751292/151‑1 (MQ=255)
                                   aaCCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACt                                                                                                           <  2:376939/151‑1 (MQ=255)
                                      ccTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGat                                                                                                        <  1:271893/151‑1 (MQ=255)
                                      ccTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGat                                                                                                        <  1:272468/151‑1 (MQ=255)
                                         aaGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGatcaa                                                                                                     <  2:468452/151‑1 (MQ=255)
                                         aaGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTAATGAGGCAAACCTTTAGCAGCATTATTTCGACGTccca                                                                                                             >  2:399382/1‑141 (MQ=255)
                                             tGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGAGGGAAACCTTTAGCAGCCTTATTTCGACGGCCAACTGatcaaaaa                                                                                                  >  2:543779/1‑147 (MQ=255)
                                              gCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGt                                                                                                <  1:346426/151‑1 (MQ=255)
                                              gCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGCGATCCTTTTTCAGCCTAATTCCGACGGCCAACTGat                                                                                                        >  1:379900/1‑143 (MQ=255)
                                                  ttGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAg                                                                                             <  1:845283/150‑1 (MQ=255)
                                                      tCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATTCGGGAAACCCTTTGCAACTTAATTTCCGCGGCCAACTGATTAAAtaggtagggt                                                                                          >  2:429481/1‑147 (MQ=255)
                                                         caGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTAc                                                                                     <  1:305668/151‑1 (MQ=255)
                                                         caGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGGGGCAATCCTTGAGCAGCCTTATTTCGACGTCCAACTGGTAAATAAATTAGGCTAtt                                                                                       >  2:464634/1‑149 (MQ=255)
                                                          aGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACa                                                                                    <  1:712671/151‑1 (MQ=255)
                                                            tAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTGATTCGGCAATCCTTTAGCAGATTAATATTCACGTCCAAATGCTCAAACAGTGATGCTAt                                                                                        >  1:788508/1‑145 (MQ=255)
                                                                gTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCa                                                                              <  1:250404/151‑1 (MQ=255)
                                                                  aaTCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgc                                                                            <  1:101141/151‑1 (MQ=255)
                                                                     cGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggc                                                                         <  2:860334/151‑1 (MQ=255)
                                                                      ggggCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATAGGGGAAAACCTTTGCAGCCTTAATTCCAACGCCAACCCGATCAAACAGTAGGGTAtt                                                                                       >  1:64158/1‑136 (MQ=255)
                                                                      ggggCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAgcggcg                                                                        <  1:524045/151‑1 (MQ=255)
                                                                            aaaaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTa                                                                  <  1:300257/151‑1 (MQ=255)
                                                                              aaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTaaa                                                                <  1:610363/151‑1 (MQ=255)
                                                                                gCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAAAGGGGGAAACCCTTGGGCGCCTAATTTCCGCGTCCCACTGATCAATTAGTTAGGGTATTACATCGTCAGCGGGGGCGGGaaa                                                                >  2:1008651/1‑149 (MQ=255)
                                                                                gCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACg                                                              <  1:9113/151‑1 (MQ=255)
                                                                                 cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTAATCGGGGATCCCTTAGCAGCCTAATTTCCACGGCCCACTTATCAAAAAGTTAAGGTCTTTCATCCTCCGcggcgggggggaa                                                                >  2:554764/1‑144 (MQ=255)
                                                                                 cccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTTATGAGGGAAACCTTTTGCAGCCTAATTTCCAAGTCCCACCGATCAAACAAATAGGATTTTTCATCCTCAGCGGCGGGGGGGAACgg                                                             >  1:928824/1‑150 (MQ=255)
                                                                                  ccGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATTTGGCAATCCTTTAGCAGCCTAATTTTGACGCCCCAATGATAAAACAATTAGGGTATTTAAATGTCAGCCGCGGGGGGaaaa                                                               >  2:273897/1‑147 (MQ=255)
                                                                                      cTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGCAAACCTTTAGCAGCCTAATTTCCACCTCCAACTGATCAATAAATTTGGCTTTTTAATTCTCAGcggggggggggaaagaatt                                                          >  2:165128/1‑139 (MQ=255)
                                                                                       ttttGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGCGGGAAACCTTTTGCAGCCTAATTCCCACGCTCCAATGATCAAAAAATTTTGCTATTTCATCTTCAGCCGCGGGGGGAAAAGGAtt                                                          >  2:226831/1‑148 (MQ=255)
                                                                                        tttGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGt                                                       <  1:261845/150‑1 (MQ=255)
                                                                                          tGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTcc                                                     <  2:493248/150‑1 (MQ=255)
                                                                                           gCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGGTTAATCGGCAAATCTTTAGCACCCTAATTTCCAACTCCCACCGATCAATCACATAGGGTATTTAATCGTTCGCGGCGGCGGGGAACGTATTTGt                                                       >  1:662777/1‑147 (MQ=255)
                                                                                           gCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCga                                                   <  1:369921/151‑1 (MQ=255)
                                                                                            cGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCgat                                                  <  2:299643/151‑1 (MQ=255)
                                                                                             gATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCgatg                                                 <  1:628497/151‑1 (MQ=255)
                                                                                               tGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTTATGCAGGAAAACTTTAGGAGCTTAATTTCGACGTCCCACTGATAAAACAGTTATGGTATTATATCGTCAAAGGCGGCGGGGAAAGTATTCGttccctg                                                 >  1:878551/1‑143 (MQ=255)
                                                                                                gAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAAAGAGGGAAACCTTTAGCAGCCTTATTTCCACGTCCAAATTATCAATAAAATATGCTATTATATTTTCAGCGGCGGCGGGAAAAGAATTTTTCCGATGAtt                                              >  2:353963/1‑151 (MQ=255)
                                                                                                   tttAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGTTACTATGGGGAAACCTTTGCAGCCTAATTCCGACGGCCAACTGATCAATCAGTTAGGGTCTTAACTCCTCAGCGGCGGCGGGAAAcacatttttcc                                                     >  2:966115/1‑131 (MQ=255)
                                                                                                    ttAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCt                                           <  2:679863/150‑1 (MQ=255)
                                                                                                      aGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAAATAGGCAAACACTTAGAAACCCAATTTCGTCCGCCCACCGATCAATCAGTTAGGCTATTACATCCTCCGCGGCCGCGGGGAACGAATTTGTCCCTTGATTGCTGt                                         >  1:79466/1‑150 (MQ=255)
                                                                                                             aTCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTg                                 <  1:197543/151‑1 (MQ=255)
                                                                                                             aTCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTg                                 <  2:50501/151‑1 (MQ=255)
                                                                                                             aTCAATGAGTTAACCGCGCAAGATATGGCGGGGGGTAATGCGGGAAAACTTTAGGAGCTTTATTTCCCCGTCCAACTCATCAACCAGATAGGGTATTTTATCGTCAGCGGGGGCGGGAAACGTATTTGTCCCATGATTGTTGTTCTggg                                   >  2:792438/1‑148 (MQ=255)
                                                                                                                   gAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGa                            <  2:579721/150‑1 (MQ=255)
                                                                                                                                             gtgtTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCa  <  1:609981/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                          
TTTTTCAGTACAATCACCCTATATTCCTAACCATAAACCCTAAGTTGCCTTTGTTCACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCA  >  NC_000913/3333616‑3333906

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