Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,443,2770AC44.9% 106.9 / 25.4 69V274G (GTA→GGA) secYpreprotein translocase membrane subunit
Reads supporting (aligned to +/- strand):  ref base A (35/3);  new base C (0/31);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.50e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443127‑3443426
                                                                                                                                                        |                                                                                                                                                     
gCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTa                                                                                                                                                       <  2:75505/151‑1 (MQ=255)
 cAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTAc                                                                                                                                                      <  2:399503/151‑1 (MQ=255)
       cATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGc                                                                                                                                                >  2:118734/1‑151 (MQ=255)
         tAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGcc                                                                                                                                               >  1:157007/1‑150 (MQ=255)
          agagTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCa                                                                                                                                              <  1:393461/150‑1 (MQ=255)
              tAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTc                                                                                                                                         >  2:584153/1‑151 (MQ=255)
                         ggggTTGCCCAGGCTGCAAATACAGCGAATTTGTTCCCACCCGGTCCCAA‑CCGGTACCGCCCCCGACCCATGACGGGAGGTTCGCCGGGAAAAAAATAATACGGG‑AAGCAAAGATTGCCGGAATCCCCCCCGCCATATTCACTTTCAGCgg                                                                                                                              <  2:387636/149‑1 (MQ=255)
                           gtttccccAGGCTGAAAATACAGCAAATTTTTTGCCAGCCGTTCCCAC‑CCAGTCCCGCCCCCGAACCATGACGCGTTGGTCGCCGGAAAAAAAATAATACTGG‑AAGAGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                            <  2:65457/146‑1 (MQ=255)
                           ggTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCgg                                                                                                                              >  2:107430/1‑149 (MQ=255)
                           ggTTCCCCGGGCTGCAAATACAGCAAAATTTTTGCCAGCCAGTTCCAA‑CCAGTCCCGCCCCCGACCCATGACGCGATGGTCGCGGGGAACAAAATAATACTGG‑AAGCGAAGATTGCGGGGATCCCCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                            <  2:459665/151‑1 (MQ=255)
                            gTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTa                                                                                                                            >  1:400347/1‑150 (MQ=255)
                                cccAGGCGGAAAAAACAGGAAATTTGTTGTCGGCCATTCCCAC‑CCAGTCCCGCCCCCGAACCATGCCGCGATGGTCGCGGGGAAAAGAAAAATACGGG‑AGGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGt                                                                                                                             <  2:811496/145‑1 (MQ=255)
                                   aGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAAtgtgtg                                                                                                                    >  1:570298/1‑151 (MQ=255)
                                    ggCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGc                                                                                                                   >  1:1041959/1‑151 (MQ=255)
                                     gCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGct                                                                                                                  >  1:869545/1‑151 (MQ=255)
                                          aaaTACAGCAAATTTTTTGCCAGCCA‑TTCCAACCCAGTCCCCCCCCCGACCCAGGCCGCGATGGCCGCGGGAAACAAAATAATACTGG‑AAGCGAAGATGGCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGctct                                                                                                                <  2:736617/148‑1 (MQ=255)
                                           aaaccagaaaagttgttgCCAGCCAGTTCCAC‑CCACCCCCCCCCCCAAACCATGCCGCGAGGGCCGCGGGAAAAAAAATAATACTGG‑AAGGGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCtg                                                                                                               <  1:781275/136‑1 (MQ=255)
                                           aaaaCAGCAAATTTTTTGTCCCCCA‑TTCCAACCCAGTCCCCCCCCCAACCCATGACGCGAGGTTCGCGGGAAACAAAATAATACTGG‑AAGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTgc                                                                                                            <  2:797664/148‑1 (MQ=255)
                                                   gAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCAt                                                                                                       >  2:68061/1‑148 (MQ=255)
                                                     aaTTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGac                                                                                                   >  2:758769/1‑150 (MQ=255)
                                                      aTTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                 >  1:367027/1‑151 (MQ=255)
                                                       tttttgtCACCTATTCCCAC‑CAAGTGCCGCCCCCGACCCAGAACGGGATGGGCGCGGGAAAAAAAAAAATAGTGG‑AAGGGAAGTGGGCGGGTATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGacac                                                                                                 <  1:897439/147‑1 (MQ=255)
                                                          ttgtCCGCCATTCCCAC‑CCAGTCCCCCCCCCGAACCATGCCGCGATGGTCGCCGGAAACAAAATAATACTGG‑AAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                              <  2:999163/150‑1 (MQ=255)
                                                           tgccacccAGTCCAAC‑CCAGTCCCCCCCCCAACCCATGCCGCGTTGGTCGCCGGGAACAAAATAATCCTGA‑AAGCGAAGATGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACacgac                                                                                              <  1:464760/143‑1 (MQ=255)
                                                             tCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGAcc                                                                                          >  2:117687/1‑151 (MQ=255)
                                                              cAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCt                                                                                         >  2:400901/1‑151 (MQ=255)
                                                               aGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        >  1:702362/1‑151 (MQ=255)
                                                                 ccattcccAC‑CCAGTCCCCCCCCCAAACCATGCCGCGAGGGTCGCGGGAAACAAAATAATACGGG‑AAGCGAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCtt                                                                                        <  2:875512/143‑1 (MQ=255)
                                                                  cAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCt                                                                                     >  1:384523/1‑151 (MQ=255)
                                                                         aa‑CCATTCCCCCCCCCAACCAAGGCCGCGATGGTCGCGGGAAAAAAAATAATACTGG‑AAGCAAAGATTCCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                              <  2:213489/151‑1 (MQ=255)
                                                                         aa‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGttt                                                                              >  2:690715/1‑151 (MQ=255)
                                                                            ccAGTCCCCCCCCCGACCCTTCACGCGAGGCTCGCGGGGAAAAAAATAATCCTGG‑AAGAAAATATTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                             <  2:482898/150‑1 (MQ=255)
                                                                              aGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTcgcg                                                                          >  1:777069/1‑151 (MQ=255)
                                                                                tACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGt                                                                         >  1:606508/1‑150 (MQ=255)
                                                                                 ccccccccatacccATGAGGCGATGGTCGCGGGGAACAAAATAATACTGG‑AAGCAAAGATTGCCGGTTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAg                                                                       <  1:612150/139‑1 (MQ=255)
                                                                                  ccccccccAAACCAGGCCGCGTTGGCCCCCGGAACAAAAATAAT‑CTGGAAAGCAAATTTTCCCGGATTTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTc                                                                             <  1:340207/141‑1 (MQ=255)
                                                                                   cccccccAAACCAGGACGCGATGGTCGCCGGGAACAAAAAAATACTGG‑AAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGtt                                                                     <  2:58626/149‑1 (MQ=255)
                                                                                      ccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATTTTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                  >  2:592861/1‑151 (MQ=255)
                                                                                      ccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                  >  2:705420/1‑151 (MQ=255)
                                                                                      ccccGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTAc                                                                  >  2:427484/1‑151 (MQ=255)
                                                                                            aCCATGACGCGATGGGCGCCGGAAACAAAAAAATACGGG‑AAGCGAAGTTTGCCGGATTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACcaca                                                              <  1:510898/149‑1 (MQ=255)
                                                                                             ccatggcggggggggcgcgGGGAACAGAATAATAGTGG‑AGGCGAAGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATg                                                           <  1:493823/136‑1 (MQ=255)
                                                                                               aTGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcg                                                         >  1:656761/1‑151 (MQ=255)
                                                                                                 gACGGGATGGTCGCGGGAAAAAAAAAAATACGGG‑AAGGGAAGTTTGCCGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAAt                                                            <  2:808046/146‑1 (MQ=255)
                                                                                                  aCGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                      >  1:863100/1‑151 (MQ=255)
                                                                                                   cgcgATGGTCGCCGGGAACAAAATAATACTGG‑AAGCGAAGATTGCCGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                      <  2:672725/150‑1 (MQ=255)
                                                                                                    gcgATGGCCGCGGGAAAAAAAATAATACGGG‑AAGAGAAGTTGGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATgcggcg                                                      <  1:778645/149‑1 (MQ=255)
                                                                                                     cgATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                   >  1:122725/1‑151 (MQ=255)
                                                                                                         ggTCGCGGGAAAAAGAAAAATACTGG‑AAGCGAAGTTTCCGGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                   <  1:721201/147‑1 (MQ=255)
                                                                                                         ggTCGCGGGAAAAAAAATAATGCTGG‑AAGCGAAGATGGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                   <  1:124394/147‑1 (MQ=255)
                                                                                                         ggTCGCCGGAAAAAAAATAATACTGG‑AGGCGATGATGGCCGGTTTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTg                                                   <  1:893590/147‑1 (MQ=255)
                                                                                                          gTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCAc                                              >  1:120958/1‑151 (MQ=255)
                                                                                                                ggaaaaaaaataatGCTGA‑AAGCGAGGATTGCCGGGTTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTca                                         <  1:962913/142‑1 (MQ=255)
                                                                                                                  gAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaac                                       >  1:516854/1‑150 (MQ=255)
                                                                                                                     aaaaataatACCGA‑AAGCAAATATTGCCGGCCTCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTcaacaa                                     <  1:781860/146‑1 (MQ=255)
                                                                                                                      aaaataatACTGG‑AGGCGAAGATTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACaaa                                    <  1:108334/147‑1 (MQ=255)
                                                                                                                      aGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                   >  2:828486/1‑150 (MQ=255)
                                                                                                                       aaaaaaggattgaaaGCGAAGTTTGCGGGAATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATa                                  <  2:684190/138‑1 (MQ=255)
                                                                                                                        aataatACGGG‑AGGCAAAGTTTGCGGGGATTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAAt                                   <  2:301204/148‑1 (MQ=255)
                                                                                                                           aataCTGA‑AAGCAAAGTTTGCGTGTATACCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaac                              <  1:946674/150‑1 (MQ=255)
                                                                                                                            ataCTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATacaaca                             >  1:536228/1‑150 (MQ=255)
                                                                                                                                         aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg                >  2:556779/1‑151 (MQ=255)
                                                                                                                                         aaGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTg                >  2:808434/1‑151 (MQ=255)
                                                                                                                                             ttGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaa             >  1:776139/1‑150 (MQ=255)
                                                                                                                                                ccGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac         >  1:278807/1‑151 (MQ=255)
                                                                                                                                                 ggttaTTCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatac         <  2:985198/146‑1 (MQ=255)
                                                                                                                                                 cGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatact        >  2:606601/1‑151 (MQ=255)
                                                                                                                                                  ggtataaCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCaaa            <  2:157028/140‑1 (MQ=255)
                                                                                                                                                  gggATCCCCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaata          <  2:639160/148‑1 (MQ=255)
                                                                                                                                                       tACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatact  >  1:988236/1‑151 (MQ=255)
                                                                                                                                                       tACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAaatactaatact  >  1:662931/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                     
GCAGACGCATAGAGTAACACATAAAGCGGTTGCCCAGGCTGCAAATACAGCGAAATTGTTGTCAGCCAGTTCCAA‑CCAGTACCGCCCCCGAACCATGACGCGATGGTCGCCGGGAACAGAATAATACTGG‑AAGCGAAGATTGCCGGGATTACCCCCGCCATATTCACTTTCAGCGGTAAATGTGTGCTCTGTGCAGCATAGACACGACGACCTTGCTGACGTTTCGCGTAGTTTACCACAATGCGGCGTTGACCACGCTCAACAAATACAACAAAGAACGTCACTGCAAATACTAATACT  >  NC_000913/3443127‑3443426

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