Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 116,035 G→A 5.2% C355C (TGC→TGT hofB ← T2SE secretion family protein; P‑loop ATPase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913116,0350GA 5.2% 242.6 / 2.1 78C355C (TGC→TGThofBT2SE secretion family protein; P‑loop ATPase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (39/33);  new base A (2/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

ATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTG  >  NC_000913/115889‑116181
                                                                                                                                                  |                                                                                                                                                  
aTGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGc                                                                                                                                                   >  5:211420/1‑148 (MQ=255)
   ggCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCaaa                                                                                                                                                <  2:323249/148‑1 (MQ=255)
     cAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGtt                                                                                                                                             >  4:72393/1‑149 (MQ=255)
       gAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTg                                                                                                                                           >  2:5300/1‑149 (MQ=255)
         aCTTCAAATAAGGCCGTACGTACATAAACACCGTGGTAGCAATGTACACAACAGGGTGCCTGCCAGTTGTGCAGCGGCGATGGCCATACATTGTCTGTAATGTTGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCg                                                                                                                                         <  2:268625/149‑1 (MQ=255)
          cTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGt                                                                                                                                        <  3:72393/149‑1 (MQ=255)
           ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGtc                                                                                                                                       >  7:186540/1‑148 (MQ=255)
           ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTa                                                                                                                                       >  3:250848/1‑149 (MQ=255)
           ttCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTa                                                                                                                                       >  2:351404/1‑149 (MQ=255)
           ttCAAATAAGGCCGTACGACCATAAAAACCGGGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTa                                                                                                                                       >  1:224042/1‑149 (MQ=255)
             cAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAc                                                                                                                                      >  7:199352/1‑148 (MQ=255)
                     gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTg                                                                                                                             <  2:224042/149‑1 (MQ=255)
                     gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTg                                                                                                                             <  7:352827/149‑1 (MQ=255)
                     gCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTg                                                                                                                             <  7:48386/149‑1 (MQ=255)
                       cGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTggg                                                                                                                           <  1:5300/149‑1 (MQ=255)
                         tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGc                                                                                                                                          >  4:309828/1‑132 (MQ=255)
                         tACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGc                                                                                                                                          <  3:309828/132‑1 (MQ=255)
                             aCCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAAAGAGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTa                                                                                                                     >  3:222687/1‑149 (MQ=255)
                                aTAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGTCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTTGGATATTACCa                                                                                                                  <  8:194654/149‑1 (MQ=255)
                                aTAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTAc                                                                                                                    >  2:348593/1‑147 (MQ=255)
                                 tAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAg                                                                                                                 >  4:37367/1‑149 (MQ=255)
                                  aaaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCTGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGc                                                                                                                >  3:216061/1‑149 (MQ=255)
                                   aaaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGc                                                                                                                <  6:259583/148‑1 (MQ=255)
                                    aaaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCg                                                                                                               >  5:150118/1‑148 (MQ=255)
                                     aaCCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGt                                                                                                              >  3:256046/1‑148 (MQ=255)
                                             tAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCga                                                                                                     <  1:348593/149‑1 (MQ=255)
                                             tAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCga                                                                                                     >  3:230063/1‑149 (MQ=255)
                                                             cGGGTGCCTGCCAGGGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTTGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCgg                                                                                     <  5:276676/149‑1 (MQ=255)
                                                                    ctgccAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGAccc                                                                              >  4:193454/1‑149 (MQ=255)
                                                                       ccAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGCCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCa                                                                           >  6:51694/1‑149 (MQ=255)
                                                                           tGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTg                                                                       <  6:150118/149‑1 (MQ=255)
                                                                             gggCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACTCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCt                                                                     <  1:351404/149‑1 (MQ=255)
                                                                               gCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTg                                                                    >  2:162276/1‑148 (MQ=255)
                                                                                  gcggcgATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGACAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAg                                                                                        >  5:337515/1‑125 (MQ=255)
                                                                                  gcggcgATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGACAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAg                                                                                        >  5:342869/1‑125 (MQ=255)
                                                                                  gcggcgATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGACAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAg                                                                                        <  6:337515/125‑1 (MQ=255)
                                                                                  gcggcgATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGACAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAg                                                                                        <  6:342869/125‑1 (MQ=255)
                                                                                     gcgATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACg                                                              >  5:148979/1‑148 (MQ=255)
                                                                                       gATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTAc                                                           <  3:37367/149‑1 (MQ=255)
                                                                                       gATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTAc                                                           >  4:174896/1‑149 (MQ=255)
                                                                                       gATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTAc                                                           >  4:60922/1‑149 (MQ=255)
                                                                                          ggCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAg                                                        >  3:60864/1‑149 (MQ=255)
                                                                                           gCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGc                                                       >  1:357489/1‑149 (MQ=255)
                                                                                               tACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGttt                                                   <  4:216061/149‑1 (MQ=255)
                                                                                               tACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACACGCGttt                                                   >  3:9722/1‑149 (MQ=255)
                                                                                                 gaTTGTCTGGAATTTGGATGGTCTCCCCTTTCTGCCGGCGACTATGTGGGCAAAGTTTGCGTACCAGTCGCTGGGCTATTACCAGCTTAAGCGGCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTc                                                  <  4:51724/147‑1 (MQ=255)
                                                                                                  aTTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGc                                                >  2:89399/1‑149 (MQ=255)
                                                                                                       cTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTg                                           >  3:306489/1‑149 (MQ=255)
                                                                                                       cTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTg                                           >  7:156909/1‑149 (MQ=255)
                                                                                                         ggAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTgg                                          >  4:322969/1‑148 (MQ=255)
                                                                                                          gAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTAGCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTg                                           <  6:211420/146‑1 (MQ=255)
                                                                                                          gAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGaa                                        <  4:250848/149‑1 (MQ=255)
                                                                                                             tgtgGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTa                                     <  8:199352/149‑1 (MQ=255)
                                                                                                               tgGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGaa                                        <  4:90840/144‑1 (MQ=255)
                                                                                                               tgGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGaa                                        >  3:90840/1‑144 (MQ=255)
                                                                                                               tgGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTATCGCAGGTGCAATAAgt                                   >  4:32303/1‑149 (MQ=255)
                                                                                                                ggATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGCTGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAgtg                                  >  3:228475/1‑149 (MQ=255)
                                                                                                                ggATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAgtg                                  >  2:37014/1‑149 (MQ=255)
                                                                                                                ggATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAgtg                                  >  4:115695/1‑149 (MQ=255)
                                                                                                                     ggCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAg                             <  3:115695/149‑1 (MQ=255)
                                                                                                                      gCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAgg                            >  2:161565/1‑149 (MQ=255)
                                                                                                                      gCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAgg                            <  3:193454/149‑1 (MQ=255)
                                                                                                                          cccTTGCTGCCGGCGACAATTTGGGCAAAGTTTGCTTACCAGACGATGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGt                          <  8:128980/147‑1 (MQ=255)
                                                                                                                          cccTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAggg                           >  4:288417/1‑146 (MQ=255)
                                                                                                                           ccTTGCTGCCGGCGACAATGTGGGCAAATGTTGAGTACCAGACTCTGGGCTATTACCAGCGTAAGCGCCGATGTTAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGa                       <  4:9722/149‑1 (MQ=255)
                                                                                                                           ccTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCGCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGa                       <  3:60922/149‑1 (MQ=255)
                                                                                                                           ccTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGa                       <  3:322969/149‑1 (MQ=255)
                                                                                                                            cTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCa                                               >  8:378486/1‑124 (MQ=255)
                                                                                                                            cTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCa                                               <  7:378486/124‑1 (MQ=255)
                                                                                                                               gctgcCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGacaa                    <  4:60864/148‑1 (MQ=255)
                                                                                                                                 tgcCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCg                                                                                                      >  1:236499/1‑64 (MQ=255)
                                                                                                                                 tgcCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCg                                                                                                      <  2:236499/64‑1 (MQ=255)
                                                                                                                                 tgcCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACTTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGTTAGACAacac                 >  4:24344/1‑149 (MQ=255)
                                                                                                                                   ccGGCGACAATGTGGTCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCa               <  5:51694/149‑1 (MQ=255)
                                                                                                                                   ccGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAacac                 <  4:222687/147‑1 (MQ=255)
                                                                                                                                    cGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAg              <  4:306489/149‑1 (MQ=255)
                                                                                                                                      gCGACAATGTGGGCAAAGTTTGCGTACCAGACTCTGGGCTATTACCTGCGTAAGCGCCGATTATAGCATCCAGCGGGCTACCCCCATTTTCTTTAAACGTACCAGAGTTTCGCATGTGGAATTATTGTGTAGGGTAGACAACACCAGGt            <  1:170293/149‑1 (MQ=255)
                                                                                                                                          cAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGAc         >  3:220680/1‑148 (MQ=255)
                                                                                                                                                gggCAAAGTTTGCGTACCAGACGATGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGAGGGCGACACCCATTTGATGTAAAAGTAACAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACAAGGTGAACATTTTg  >  2:66067/1‑149 (MQ=255)
                                                                                                                                                  gCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGttt   >  6:361158/1‑146 (MQ=255)
                                                                                                                                                  gCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGttt   <  5:361158/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGGGCAGAACTTCAAATAAGGCCGTACGACCATAAAAACCGTGGTAGCAATGTACACAACCGGGTGCCTGCCAGTGGGGCAGCGGCGATGGCCATACATTGTCTGGAATGTGGATGGGCTCCCCTTGCTGCCGGCGACAATGTGGGCAAAGTTTGCGTACCAGACGCTGGGCTATTACCAGCGTAAGCGCCGATGATAGCATCCAGCGGGCGACCCCCATTTGCTGTAAACGTACCAGCGTTTCGCAGGTGGAATTAGTGTGTAGGGTAGACAACACCAGGTGACCAGTTTG  >  NC_000913/115889‑116181

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