Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 924,183 T→C 5.5% I307T (ATC→ACC)  clpA → ATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913924,1830TC 5.5% 236.0 / 2.2 73I307T (ATC→ACC) clpAATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity
Reads supporting (aligned to +/- strand):  ref base T (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

TCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGTTCAAA  >  NC_000913/924040‑924332
                                                                                                                                               |                                                                                                                                                     
tCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAAc                                                                                                                                                  <  2:77517/149‑1 (MQ=255)
    ggcgACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGtt                                                                                                                                              >  2:69142/1‑149 (MQ=255)
     gcgACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTg                                                                                                                                             <  3:542290/149‑1 (MQ=255)
          tttGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCtctc                                                                                                                                        >  3:7593/1‑149 (MQ=255)
          tttGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCtctc                                                                                                                                        >  1:175152/1‑149 (MQ=255)
                     ttttAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTaaaa                                                                                                                             <  3:168536/149‑1 (MQ=255)
                          aaGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTcc                                                                                                                                       >  2:546025/1‑134 (MQ=255)
                          aaGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTcc                                                                                                                                       <  1:546025/134‑1 (MQ=255)
                          aaGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGTGGTACAATTCgt                                                                                                                        <  4:584972/149‑1 (MQ=255)
                          aaGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCgt                                                                                                                        <  3:261668/149‑1 (MQ=255)
                          aaGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCgt                                                                                                                        <  2:130257/149‑1 (MQ=255)
                                tGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAAtt                                                                                                                  <  3:256571/149‑1 (MQ=255)
                                  cTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATggg                                                                                                                >  1:174354/1‑149 (MQ=255)
                                     aaGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTc                                                                                                             <  2:203346/149‑1 (MQ=255)
                                     aaGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTc                                                                                                             >  1:316519/1‑149 (MQ=255)
                                      agcagcTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCg                                                                                                            <  3:378456/149‑1 (MQ=255)
                                         agcTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACa                                                                                                         >  3:194043/1‑149 (MQ=255)
                                          gcTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACaa                                                                                                        >  3:294603/1‑149 (MQ=255)
                                                   ggACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAgg                                                                                               >  4:310627/1‑149 (MQ=255)
                                                    gACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGTGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGa                                                                                              <  2:149654/149‑1 (MQ=255)
                                                    gACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGa                                                                                              >  4:375762/1‑149 (MQ=255)
                                                    gACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAACCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGa                                                                                              >  4:4449/1‑149 (MQ=255)
                                                    gACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAACCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGa                                                                                              >  4:6348/1‑149 (MQ=255)
                                                        cTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTc                                                                                          <  2:116181/149‑1 (MQ=255)
                                                         tAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCa                                                                                         >  2:316910/1‑149 (MQ=255)
                                                         tAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCa                                                                                         <  2:476834/149‑1 (MQ=255)
                                                          aaCAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAg                                                                                        >  4:106982/1‑149 (MQ=255)
                                                               cATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAACCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACa                                                                                   <  3:4449/149‑1 (MQ=255)
                                                               cATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAACCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACa                                                                                   <  3:6348/149‑1 (MQ=255)
                                                                    tGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTc                                                                              <  2:175152/149‑1 (MQ=255)
                                                                    tGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTc                                                                              <  1:643528/149‑1 (MQ=255)
                                                                          ttGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGaaa                                                                        <  3:304998/149‑1 (MQ=255)
                                                                            gatgaGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCa                                                                                         >  3:363462/1‑130 (MQ=255)
                                                                            gatgaGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCa                                                                                         <  4:363462/130‑1 (MQ=255)
                                                                              tgaGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGAcc                                                                    >  3:179989/1‑149 (MQ=255)
                                                                               gagaTCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCg                                                                   >  2:517880/1‑149 (MQ=255)
                                                                                agaTCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGt                                                                  <  1:298436/149‑1 (MQ=255)
                                                                                      acacCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTg                                                            >  4:322775/1‑149 (MQ=255)
                                                                                        acCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGgc                                                          <  3:86665/149‑1 (MQ=255)
                                                                                          cATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGgcgc                                                        >  2:460598/1‑149 (MQ=255)
                                                                                          cATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGgcgc                                                        >  4:523674/1‑149 (MQ=255)
                                                                                          cATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGgcgc                                                        >  2:427840/1‑149 (MQ=255)
                                                                                            ttATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGgcgcat                                                      >  1:459586/1‑147 (MQ=255)
                                                                                            ttATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGt                                                      >  2:122225/1‑149 (MQ=255)
                                                                                             tATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTc                                                     <  1:460598/149‑1 (MQ=255)
                                                                                                         tGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGaaaa                                         <  4:130014/149‑1 (MQ=255)
                                                                                                            aGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTg                                      >  1:689508/1‑149 (MQ=255)
                                                                                                             gCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGa                                     <  1:56513/149‑1 (MQ=255)
                                                                                                                     tggtCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTg                             >  3:6716/1‑149 (MQ=255)
                                                                                                                       gtCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTc                                                                              >  3:359711/1‑98 (MQ=255)
                                                                                                                       gtCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTc                                                                              <  4:359711/98‑1 (MQ=255)
                                                                                                                       gtCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGaa                           >  1:171279/1‑149 (MQ=255)
                                                                                                                         cAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAAcc                         <  1:427840/149‑1 (MQ=255)
                                                                                                                          aGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCg                        >  1:88714/1‑149 (MQ=255)
                                                                                                                           ggTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGt                       <  4:7593/149‑1 (MQ=255)
                                                                                                                                   cGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               <  2:316519/149‑1 (MQ=255)
                                                                                                                                   cGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCATCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               <  2:459586/149‑1 (MQ=255)
                                                                                                                                     gCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaaga             >  1:337501/1‑149 (MQ=255)
                                                                                                                                         aCCTAATCGAACCGTTGCTCTCCAGCGGTAAAATTCGGGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               >  3:636776/1‑143 (MQ=255)
                                                                                                                                         aCCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               <  4:636776/143‑1 (MQ=255)
                                                                                                                                          ccTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               >  3:625659/1‑142 (MQ=255)
                                                                                                                                          ccTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               <  4:625659/142‑1 (MQ=255)
                                                                                                                                           cTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGt                                                                  >  3:244377/1‑90 (MQ=255)
                                                                                                                                           cTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGt                                                                  <  4:244377/90‑1 (MQ=255)
                                                                                                                                           cTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGt       >  1:117592/1‑149 (MQ=255)
                                                                                                                                           cTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGt       >  1:600928/1‑149 (MQ=255)
                                                                                                                                           cTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGt       >  3:169227/1‑149 (MQ=255)
                                                                                                                                            tAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAACGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               >  4:66452/1‑140 (MQ=255)
                                                                                                                                            tAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCgaa               <  3:66452/140‑1 (MQ=255)
                                                                                                                                            tAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGtt      <  3:103582/149‑1 (MQ=255)
                                                                                                                                            tAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGtt      <  3:523674/149‑1 (MQ=255)
                                                                                                                                            tAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGtt      >  4:103582/1‑149 (MQ=255)
                                                                                                                                             aaTCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGTTc     >  1:415439/1‑149 (MQ=255)
                                                                                                                                               tCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGTTCaa   <  2:174354/149‑1 (MQ=255)
                                                                                                                                               tCAAACCGTTGCTC‑CCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGTTCaaa  >  2:326628/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
TCGCGGCGACTTTGAAAAACGTTTTAAAGCGTTGCTCAAGCAGCTGGAGCAGGACACTAACAGCATCCTGTTTATTGATGAGATCCACACCATTATCGGTGCGGGTGCAGCGTCTGGTGGTCAGGTCGATGCGGCTAACCTAATCAAACCGTTGCTCTCCAGCGGTAAAATTCGTGTAATTGGTTCGACAACCTATCAGGAGTTCAGCAACATTTTCGAGAAAGACCGTGCTCTGGCGCGTCGCTTCCAGAAAATTGATATTACTGAACCGTCGATCGAAGAAACTGTTCAAA  >  NC_000913/924040‑924332

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