Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,2770AC47.4% 52.7 / 38.2 57V274G (GTA→GGA) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (24/6);  new base C (0/27);  total (24/33)
Fisher's exact test for biased strand distribution p-value = 1.49e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATA  >  NC_000913/3443133‑3443424
                                                                                                                                                    |                                                                                                                                                   
gCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCg                                                                                                                                               >  1:960307/1‑151 (MQ=255)
  aTAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGcc                                                                                                                                             <  2:404844/151‑1 (MQ=255)
     gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCata                                                                                                                                          >  2:996779/1‑151 (MQ=255)
     gagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAACGAAGATTGCCGGGATT‑ACCCCCGCCata                                                                                                                                          >  2:274920/1‑151 (MQ=255)
               tAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCa                                                                                                                                >  1:1008692/1‑151 (MQ=255)
                         ccccAGGCTGAAAATACAGCAAATTT‑G‑TTGCCCGCAAGTTCCAACCAGTACCCCCCCAGACCCAGG‑CCGCGGTGGCCGCCGGAAAAAAAATAATACTGGAAGCAAATATTGCCGGGATC‑ACCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                      <  2:340687/150‑1 (MQ=255)
                          cccgggcagaaaaaccAGCAAATTT‑G‑TGGCCACCCAGTTCCACCCAGTCCCCCCCCCGAACCATG‑CGGCGAGGGTCCCGGGAAACAAAATAATACGGAAAGAGAAGATGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAAtg                                                                                                                      <  2:758450/135‑1 (MQ=255)
                               ccTGAAAATCCAGCAAATTT‑T‑TTGTCACCCAGTCCCAACCATCCCCCCCCCCGAACCAGG‑CCGGGAGGGCCGCGGGGAACAAAAAAATCCTGGAGGCGAAGATTGCGGGTTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                 <  2:698155/149‑1 (MQ=255)
                                cTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                >  1:777462/1‑150 (MQ=255)
                                cTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                >  2:692355/1‑150 (MQ=255)
                                cTGAAAATACAGAGAAATT‑T‑TTGCCAGCCATTCCCAACCAGTCCCGCCCCCGACCCAGG‑CCGGGATGGCCCCGGGGAAAAGAATAATTCGGGAAGCTAAGATGGCGGGGTTT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctc                                                                                                               <  2:579705/151‑1 (MQ=255)
                                  gCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                >  2:845241/1‑148 (MQ=255)
                                            caaaaTG‑G‑TTGCCACCCAGTTCCAACCGGCCCCCCCCCCAAACCATG‑CCGGGAGGGTCGCGGGAAACAAAAAAATACGGGAAGCAAATATTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  2:671220/147‑1 (MQ=255)
                                            cGAAATT‑G‑TTGCCAGCCATTCCCAACAAGTCCCGCCCCCGACCCATG‑ACGCGATGGCCGCCGGGAACAAAATAATACTGAAAGCAAAGTTTGCGGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                       <  2:995545/147‑1 (MQ=255)
                                            agAATTT‑T‑TTGCCACCCAGTCCCAACAAGACCCGCCCCCAACCCAGG‑AGGGGAGGGTCGCGGGGAAAAAAATAATACTGGAAGC‑AAAATTTCCGGGGTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  2:923436/148‑1 (MQ=255)
                                            agAAGTT‑G‑TGGCCGGCCATTCCCACCCAGTCCCCCCCCCAAACAATG‑AGGCGTTGGCCGCCGGAAACAAAAAAATAGGGGAAGGGAAGTTTGCCGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  2:716025/148‑1 (MQ=255)
                                            aaaattt‑g‑ttgCCAGCCATTCCCACCCAGTCCCCCCCCCGACCCATG‑CCGGGATGGGCCCCGGGAACAAAATAAGACTGGAAGGGAAGATTGCGGGGTTT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  2:255496/144‑1 (MQ=255)
                                            aaaaaTT‑G‑TGCTCACCCAGTTCCAACCAGTCCCCCCCCCGACCCAGG‑CCGCGGTGGTCCCCGGGAACAAAATAATACGGGAGGAAAAATTGGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgcagc                                                                                                       <  1:48577/145‑1 (MQ=255)
                                             gAAATTGG‑TTGTCAGCCGGTTCAAACAAACCCCCCCCCGGAACCATG‑ACGCGAGGGCCGCGGGAAAAAAAAAAATTGTGGAAGAGAAGATGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  1:652666/149‑1 (MQ=255)
                                             ccaaattttttgcccgCCATTCCCACCCAGTCCCCCCCCAGACCCATG‑ACGGGAGGTCCGCGGGAAACAAAATAATCCTGGAAGCAAAGATTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  1:301738/135‑1 (MQ=255)
                                             aaaaTT‑GTTTGCCAGCAAGTCCCAACCAG‑ACCCCCCCCGGACCCGGAACGGGATGTGCCCCGGAAACAAAATAATCCGGAAGGCGAAGTTGGCCGGTATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  2:227234/148‑1 (MQ=255)
                                              aaaattttttgcCACCCAGTCCCAACCAGCACCGCCCCCAACCAATG‑ACGCGAGGGTCCCGGGAAAAAAAATAATTGGGGAAGAAAAGTTTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                     <  1:421179/141‑1 (MQ=255)
                                              aaaTT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                  >  2:534808/1‑150 (MQ=255)
                                              aaaTT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGa                                                                                                  >  2:162264/1‑150 (MQ=255)
                                               aaTT‑G‑TTGTCAGCCAGTCCCACCCAGTCCCCCCCCCGACCCAGG‑ACGGGTGGGTCGCGGGAAAAAAAATAATTCGGGAAGAGAAGTTGGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGaca                                                                                                <  1:605100/151‑1 (MQ=255)
                                                    g‑ttgtCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacga                                                                                             >  1:366650/1‑150 (MQ=255)
                                                      ttgtCAGCCAGTCCAAACCAGCCCCCCCCCAGACCCATG‑GCGCGTGGGTCGCCGGGAACAAAATAATAGGGGAAGAGAATTGGGCCGGTATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgacg                                                                                           <  1:512802/151‑1 (MQ=255)
                                                      tggccAGCCAGTCCCACACAGCCCCGCCCCAGAACCAGG‑CCGGGTTGGTCGCGGAAAACAAAAAAATCCGGGAGGCGAATTTTGCGGGTATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                            <  1:801651/146‑1 (MQ=255)
                                                        tcccgCCAGTCCCACCCAACCCCGCCCCCAAACCAGG‑CCGGGGGGGTCGCCGGAACAAAAAAAATCCTGGAGGCGAAGTTTCCGGGAATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                         <  2:517699/147‑1 (MQ=255)
                                                         tCAGCCAGTTCCACCCAGCCCCGCCCCCGACCCATG‑CCGCGATGGCCCCCGGGAAAAAAAAAATACTGGAAGCAAAGATTGCCGGAATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                        <  1:612049/151‑1 (MQ=255)
                                                         gcggCCGGTTCCACCAGGTCCCCCCCCCAAACAATG‑CCGCAAGGGCCGCCGGAACCAAAATAATACTGAAAGGGAAGTTGGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACgacgac                                                                                         <  2:97456/147‑1 (MQ=255)
                                                              cAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTg                                                                                     >  1:789257/1‑149 (MQ=255)
                                                                gTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGa                                                                                 >  1:252639/1‑151 (MQ=255)
                                                                   ccAACCAGTCCCCCCCCCAAACAATG‑AGGGGGGGGCCCCGGGAAACAAAATAATCCTGGAAGGGAAGATTGCGGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGt                                                                              <  1:938313/151‑1 (MQ=255)
                                                                    cAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGtt                                                                             >  1:956765/1‑151 (MQ=255)
                                                                      aCCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTACTGACGttt                                                                            >  1:690070/1‑150 (MQ=255)
                                                                               ccccccGACCCAGG‑CCGGGGTGCCCGCCGGAAAAAAAATAATATGGGAAGGGAAGTTTCCCGGAATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGttt                                                                  <  2:228149/150‑1 (MQ=255)
                                                                                cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                                 >  1:50463/1‑151 (MQ=255)
                                                                                cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                                 >  2:441161/1‑151 (MQ=255)
                                                                                cccccGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTa                                                                 >  2:547852/1‑151 (MQ=255)
                                                                                    caaaCCAGG‑ACGCGAGGGTCGCCGGAAAAAAAAAAATACTGGAAGCAAAGATTGCCGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcac                                                             <  2:868452/149‑1 (MQ=255)
                                                                                         cATG‑AGGGGATGTGCGCGGGAAACAAAATAATACGGAAAGAGAATTTTGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgc                                                        <  1:534574/151‑1 (MQ=255)
                                                                                              aCGCGATGGTCGCCGGAAACAAAATAATACTGGAAGCGAAGTTTGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggc                                                     <  1:363199/150‑1 (MQ=255)
                                                                                                    tgctcccggggAACAAAATAATACTGAACGAAAATTTGCCCGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGAc                                               <  1:171315/142‑1 (MQ=255)
                                                                                                               aacaaaaaaatcggggAAGGGAAATTTGCGGGGTTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                   <  2:218429/137‑1 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                >  2:95915/1‑150 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                >  2:631381/1‑150 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                >  1:307731/1‑150 (MQ=255)
                                                                                                                   gAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCCCAACAAATa                                >  2:172404/1‑150 (MQ=255)
                                                                                                                        ataCTGGAGGGAAAGAGTGCCGGGATC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATAcaacaa                          <  2:794381/151‑1 (MQ=255)
                                                                                                                            tggagggaaaGATGGCGGGATTT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGaa                      <  2:899456/143‑1 (MQ=255)
                                                                                                                                 gCGAAGTTGGCCGGTTTC‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTc                  <  1:928711/150‑1 (MQ=255)
                                                                                                                                     aGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg              >  1:718043/1‑150 (MQ=255)
                                                                                                                                       ttgcccGGTATC‑ACCCCCGCCATATTCNNTTNCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc             <  1:800127/145‑1 (MQ=255)
                                                                                                                                       aTTGCGGGGATT‑CCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGc             <  2:22840/149‑1 (MQ=255)
                                                                                                                                          gCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata        >  1:883645/1‑151 (MQ=255)
                                                                                                                                              ggATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaa    >  1:65131/1‑151 (MQ=255)
                                                                                                                                                    aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaata  >  2:619893/1‑148 (MQ=255)
                                                                                                                                                    aCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaata  >  1:979357/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATT‑G‑TTGTCAGCCAGTTCCAACCAGTACCGCCCCCGAACCATG‑ACGCGATGGTCGCCGGGAACAGAATAATACTGGAAGCGAAGATTGCCGGGATT‑ACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATA  >  NC_000913/3443133‑3443424

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