Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,263,502 T→C 7.7% A53A (GCA→GCG)  tdcD ← propionate kinase/acetate kinase C, anaerobic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,263,5020TC 7.7% 323.0 / 14.9 118A53A (GCA→GCG) tdcDpropionate kinase/acetate kinase C, anaerobic
Reads supporting (aligned to +/- strand):  ref base T (46/60);  new base C (7/2);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 7.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.42e-01

Only 100 of 121 total aligned reads displayed.

TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAA  >  NC_000913/3263353‑3263651
                                                                                                                                                     |                                                                                                                                                       
tCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGt                                                                                                                                                         <  1:1336306/150‑1 (MQ=255)
 cATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTg                                                                                                                                                        <  2:161864/150‑1 (MQ=255)
  agcGGAAATAATGGCTGACTCTGTAAAAATACTGCCGCCGTGATCGAGTCGGTGGCCAATTTAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGc                                                                                                                                                       <  2:301753/148‑1 (MQ=255)
    cGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTg                                                                                                                                                     >  1:646632/1‑150 (MQ=255)
     ggTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTgg                                                                                                                                                    >  1:1121505/1‑150 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGCTGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCCATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                   <  1:1485690/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                   <  1:16134/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                   <  2:646632/150‑1 (MQ=255)
      gTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGGGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGc                                                                                                                                                   <  1:529403/150‑1 (MQ=255)
           aatGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTcc                                                                                                                                              <  1:1549465/150‑1 (MQ=255)
            atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCTCACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTccc                                                                                                                                             >  1:302854/1‑150 (MQ=255)
            atGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTccc                                                                                                                                             <  1:1435953/150‑1 (MQ=255)
             tGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCTAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAACGGTGCTGGatcaaccc                                                                                                                                            >  2:1456856/1‑145 (MQ=255)
                cGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGACAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACAGTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑At                                                                                                                                        >  2:1110782/1‑150 (MQ=255)
                 ggACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGATGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑Att                                                                                                                                       <  1:1537311/150‑1 (MQ=255)
                  gACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCATTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCCGTTGTGGGACTGCCCC‑Attt                                                                                                                                      >  2:1586679/1‑150 (MQ=255)
                    cTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTAc                                                                                                                                    <  1:212578/150‑1 (MQ=255)
                      cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCACCCGTGGTGGGGCTCCCCC‑AATTAAgg                                                                                                                                  >  2:1448372/1‑150 (MQ=255)
                      cGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACgg                                                                                                                                  <  2:1552305/150‑1 (MQ=255)
                        gTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGATCCCGCGGTGCTGGAACTCCCCCAATTTACgtt                                                                                                                                 >  2:568247/1‑148 (MQ=255)
                            aaaTACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGa                                                                                                                            >  2:552565/1‑150 (MQ=255)
                            aaaTAATGCCGCCGTGAGCGATGCGTGGGCCAATTAGGGCCACCATGGCTTTTAAATGCGTTTTTGCCATTTAAATTGAAATTGCATTCAATGCCCCTTCTTAGCGTTGGTGAGCAAGCGGTCATGGATCTCCCCC‑ATTTACGGATacta                                                                                                                            >  2:1216430/1‑147 (MQ=255)
                               tACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGTATg                                                                                                                         <  2:1564986/150‑1 (MQ=255)
                                aCTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGGTGGGTCGCCCCC‑AATTACGGATAAGAATgg                                                                                                                        >  2:374799/1‑149 (MQ=255)
                                aCTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGc                                                                                                                        <  1:1208795/150‑1 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCCCGGGGGGGGGGGTCCCCC‑ATTTTAGGATTAGAATGGAt‑tt                                                                                                                   >  1:1063222/1‑150 (MQ=255)
                                    ccgccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGATCCCCC‑CATTAACGTTAATAATGGAt‑tt                                                                                                                   >  2:597037/1‑150 (MQ=255)
                                     cgccgTGCGCGATGCGGTGGCCAATTAAGGCCACACTGTCAATTAAATTCCGTTGTGCCAGTTCAACTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCGCCCCC‑ATTTACGGATAAGAATGCAt‑ttt                                                                                                                  <  1:287197/150‑1 (MQ=255)
                                      gccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTc                                                                                                                 <  1:990538/150‑1 (MQ=255)
                                       ccgTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTcc                                                                                                                >  2:57764/1‑150 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGGGGTGCTGGGGCTCCCCC‑AATTACGGGTTAGAATGCAT‑TTTGCgg                                                                                                              >  1:414324/1‑149 (MQ=255)
                                         gTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGa                                                                                                              <  2:162955/150‑1 (MQ=255)
                                           gAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGTATGCAT‑TTTCCGAGt                                                                                                            >  1:714026/1‑150 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCTTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGtt                                                                                                           <  1:1453424/150‑1 (MQ=255)
                                            aGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGtt                                                                                                           <  1:1111166/150‑1 (MQ=255)
                                                   cGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCg                                                                                                    <  1:1019000/150‑1 (MQ=255)
                                                    ggTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGGGGTGGCTCTACCCC‑ATTTTCGGATTATAATGCAT‑TTTACGAGTTAATAACgg                                                                                                   >  1:1703449/1‑149 (MQ=255)
                                                      tGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCg                                                                                                 >  1:1015107/1‑150 (MQ=255)
                                                       ggCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCgg                                                                                                >  1:920587/1‑150 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGc                                                                                               <  1:1415423/150‑1 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGc                                                                                               <  1:825709/150‑1 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGc                                                                                               <  2:195188/150‑1 (MQ=255)
                                                        gCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCACTACCCC‑ATTTACGGATAAGAAAGCAT‑TTTACGATTTAATACCTTCgga                                                                                               >  2:297380/1‑149 (MQ=255)
                                                           aaTTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAAt                                                                                            >  1:36153/1‑150 (MQ=255)
                                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGCTGGGTCTCCCCC‑ATTTAAGGATAAAAATTCAT‑TTTACGAGTTTATAACGTAGGCAATACCAGa                                                                                      >  2:491064/1‑150 (MQ=255)
                                                                 ggCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGa                                                                                      <  1:1448372/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                     <  2:1453338/150‑1 (MQ=255)
                                                                  gCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGAc                                                                                     <  1:606572/150‑1 (MQ=255)
                                                                   ccACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACa                                                                                    <  1:584359/150‑1 (MQ=255)
                                                                     acacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCTGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAAAGCAT‑TTTCCGAGTTAATACCGTCGGCAATTCCTGACAtt                                                                                  >  1:242259/1‑150 (MQ=255)
                                                                      cacTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCACCGGTGGGTGCGGACCCCC‑AGTTAAAGAAAAGAATGCATGTTAACAGGTTAAAACCG‑GGGGAATACCTGACagtt                                                                                 >  2:871129/1‑147 (MQ=255)
                                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                               <  2:1121505/150‑1 (MQ=255)
                                                                        cTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAAAGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAAt                                                                               <  1:459365/150‑1 (MQ=255)
                                                                            cATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACtt                                                                           >  1:1020699/1‑150 (MQ=255)
                                                                              tttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGAAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTca                                                                         <  2:723895/150‑1 (MQ=255)
                                                                               ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                        <  2:63964/150‑1 (MQ=255)
                                                                               ttAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTcac                                                                        >  2:236613/1‑150 (MQ=255)
                                                                                tAAATTCCGTTTTTCCAGTTCAAATGAAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCGCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCAGCAATACCTGACATTAATACTTcaca                                                                       <  2:319279/150‑1 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTATACCCC‑AATTAAGGATAAGAATGAAT‑TTTACGAGTTAATACCGTCGGCAAAACATGACAATAATAAATAACAGTCGCTg                                                                >  1:1034137/1‑150 (MQ=255)
                                                                                       cGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTg                                                                >  1:986154/1‑150 (MQ=255)
                                                                                        gTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTgg                                                               <  1:980228/150‑1 (MQ=255)
                                                                                           tttCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCAt                                                            <  1:1271494/150‑1 (MQ=255)
                                                                                              ccAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGa                                                         <  2:217589/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCTTATCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTGGCATCGAGc                                                       <  1:297380/150‑1 (MQ=255)
                                                                                                aGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCACCGGGGCTGGCTCTCCCCC‑CATTAAGGATAATAAATCAA‑TTTAAGAGTTAATACCGTCGGCCATACCGGACAATAATAAATAACAATAGCAGGCATCGAGc                                                       >  2:385737/1‑150 (MQ=255)
                                                                                                  ttCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTTGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAAAACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCAc                                                     <  2:748133/150‑1 (MQ=255)
                                                                                                   tCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGCGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAAAACCGTCGGCAATACCAGACCTTAAAACAGCACAGTCGCAGGCATCGAGCaag                                                    >  2:894179/1‑148 (MQ=255)
                                                                                                       aTGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑AATTACGGATAAGAATGCAT‑TTTACGAGTTAATACCGTCGGCAATAACTGACATTAATAAATAAAACTCGCTGGCATCGAGCACGGaaa                                                >  2:1609291/1‑150 (MQ=255)
                                                                                                         gCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCCGCGGTGGTGGCTCTCCCCC‑ATTTAAGGATAAGAATGCAT‑TTTCCGAGTTAAAACCGTGGGCAATTCCTGACATTAAAACGTCACAGTAGCTGGCCTAGAGCCCGGAAAAc                                              >  2:141996/1‑150 (MQ=255)
                                                                                                          cAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑AATTAAGGATAAGAAAGCAT‑TTTACGAGTTAAAACCCTCGGCAAAACCTTAAAATAATAAATAACAGTAGCTTGCCTCGAGCACGGAAAAca                                             >  2:453589/1‑149 (MQ=255)
                                                                                                           aaTTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCGCTCCCCC‑ATTTACGGATAAGAATGCAA‑TTTACGAGTTTATACCGGCGGGAATACCTGGCATTAATAAGTCACAGTCGCTGGCAGCGAGAACGGAAAAcgt                                            >  2:833930/1‑148 (MQ=255)
                                                                                                              tGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAAt                                         <  1:1456856/150‑1 (MQ=255)
                                                                                                               gCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATc                                        <  2:665156/150‑1 (MQ=255)
                                                                                                                 cTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGa                                      <  1:236613/150‑1 (MQ=255)
                                                                                                                   tCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAg                                    <  1:788308/150‑1 (MQ=255)
                                                                                                                       tGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGa                                <  1:22533/150‑1 (MQ=255)
                                                                                                                            cTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACcac                           >  1:1124063/1‑150 (MQ=255)
                                                                                                                                  aGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACGTCAAAGTCGCTGGCATCGAGCACGGAAAACTTAATAGAAGTCGCACCACACTtaa                     >  1:909727/1‑150 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                 <  1:552565/150‑1 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                 >  1:362307/1‑150 (MQ=255)
                                                                                                                                      gtggtgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                 <  1:453589/150‑1 (MQ=255)
                                                                                                                                      ggggagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCGGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                 >  2:278721/6‑150 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                                   >  2:1273074/5‑131 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGa                                   <  1:1273074/127‑1 (MQ=255)
                                                                                                                                       gggagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAcc                >  1:635309/5‑150 (MQ=255)
                                                                                                                                         gtgCGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa              <  2:1501750/150‑1 (MQ=255)
                                                                                                                                         gtgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa              >  1:107929/1‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa              >  1:633492/3‑150 (MQ=255)
                                                                                                                                         gagAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaa              >  1:730257/3‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGCATGCAT‑TTTCCGAGTTAATACCGTCGGCAATTCCTGCCATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa             >  1:77555/1‑150 (MQ=255)
                                                                                                                                          tgAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa             <  1:491064/150‑1 (MQ=255)
                                                                                                                                          agAGCCAGCACCGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCaaa             >  1:122473/2‑150 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca          <  2:1703449/150‑1 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca          <  2:134618/150‑1 (MQ=255)
                                                                                                                                             gCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaaca          <  2:886878/150‑1 (MQ=255)
                                                                                                                                              ccAGCGGTGCTGGATCTCCCCC‑ATTTACGTATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAaacaa         <  1:1110782/150‑1 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAAc                 >  2:867777/1‑140 (MQ=255)
                                                                                                                                                aGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATAAATACTTCACAGTCGCTCGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAAcc       <  1:57764/150‑1 (MQ=255)
                                                                                                                                                     tGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACCGTTAATAACCAACACAACCGGaaa  >  1:556419/1‑150 (MQ=255)
                                                                                                                                                     tGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGaaa  >  1:563192/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                       
TCATCGGTAATAATGGCGGACTCGGTAAAAATACTGCCGCCGTGAGCGATGCGGTGGCCAATTAAGGCCACACTGTCATTTAAATTCCGTTTTTCCAGTTCAAATGCAATTGCCTTCAATGCACCTTCGTAGCTGTGGTGAGCCAGCGGTGCTGGCTCTCCCCC‑ATTTACGGATAAGAATGCAT‑TTTCCGAGTTAATACCGTCGGCAATACCTGACATTAATACTTCACAGTCGCTGGCATCGAGCACGGAAAACTTAATCGAAGACGAACCACAGTTAATAACCAAAACAACCGGAAA  >  NC_000913/3263353‑3263651

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