Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,959 T→G 8.2% N110T (AAC→ACC)  hslV ← peptidase component of the HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,9590TG 8.2% 166.2 / 6.4 74N110T (AAC→ACC) hslVpeptidase component of the HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base G (4/2);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

CGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAA  >  NC_000913/4121816‑4122107
                                                                                                                                               |                                                                                                                                                    
cGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCg                                                                                                                                                 <  2:1171436/149‑1 (MQ=255)
 gCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTTGCTGCACCACGTCAACGtttcagt                                                                                                                                                >  2:534266/1‑144 (MQ=255)
  ccTTCAATATCCAACGCGTTTTCAGCAATTTAACGGGCGCTAAGTTCAGGGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                               <  2:917532/149‑1 (MQ=255)
  ccTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGt                                                                                                                                               <  2:1126538/149‑1 (MQ=255)
   cTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTTCCGGtg                                                                                                                                              >  2:449867/1‑149 (MQ=255)
     gCAATATCCAACGCCTTTTCAGCAATATCACGGGCGCTATGTTCAGTGTTTTCTAACATCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgat                                                                                                                                            <  1:1112484/149‑1 (MQ=255)
      cAATATCCAACGCCTTTTCAGCAATTTCACGGTCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTATGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGtgatg                                                                                                                                           <  2:677869/149‑1 (MQ=255)
      cAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGGTAACGGtgacg                                                                                                                                           <  2:887029/149‑3 (MQ=255)
        atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTACGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTgatgat                                                                                                                                         >  2:126568/1‑149 (MQ=255)
        atatCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTgatgat                                                                                                                                         <  1:1013732/149‑1 (MQ=255)
           tCCAACGCGTTTTCAGCTATTGGACGGGCGCTAAGTTCAGTGTTTTCTATCAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTTCACCACGTCACCGTTACCGGTGATGATAAg                                                                                                                                      <  2:976055/149‑1 (MQ=255)
           tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAg                                                                                                                                      <  1:1125035/149‑1 (MQ=255)
           tCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGAGTTTTTTAACAGCGCGCGCGCCGCAGCCTGGGCGAAAGCGCCGCCGGAGCCGAAAGCAATAAGATCGTATTCTGGCTGCACAACGTCACCGTTACCGGTGATGATAAg                                                                                                                                      >  2:252031/1‑149 (MQ=255)
           tCCAACGACTATTCAGCATTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTTATGATAAg                                                                                                                                      <  2:1146436/149‑1 (MQ=255)
             cAACGCCTTTTCAGCAATTTCCCGGGCGCTAAGTTCAGTGTATTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATAGTATTCTGCCTGCACCACGTCACCGATACCGGTGATGATAAGcc                                                                                                                                    >  1:716545/1‑148 (MQ=255)
             aaaagCCTTTTCTGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGAAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGTTAAGCg                                                                                                                                    <  2:835701/145‑1 (MQ=255)
                cGCCTTTTCTGCAATTTCCCGGGCGCTCAGTGCAGTGTTTTCTAACAGCTCGCGCGACGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATg                                                                                                                                 <  1:811651/149‑1 (MQ=255)
                cGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATg                                                                                                                                 <  1:895670/149‑1 (MQ=255)
                   cTTTTCAGCAATATCACGGGCGCTAAGTTCAGTGTTTACTAACAGAGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATCAGCGGTGCAg                                                                                                                              >  2:380208/1‑149 (MQ=255)
                    ttttCAGCAATTTCACGGGCGCTACGTTCAGTGTTTTCTCACCGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGCGCCGATAGCAATAAGAGCGTTTTCTGGCTGCAGCACGTCACCGTGACCGGTGATGATAAGCGCTGCAGt                                                                                                                             >  1:811077/1‑149 (MQ=255)
                    ttttCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGt                                                                                                                             >  1:674997/1‑149 (MQ=255)
                        cAGCAATTTCCCGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCa                                                                                                                         <  2:369861/149‑1 (MQ=255)
                           cAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATcc                                                                                                                      <  1:73459/149‑1 (MQ=255)
                            aaTTTCACGGGCGCTAAGTTCAGTCTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCcg                                                                                                                     >  1:399487/1‑149 (MQ=255)
                             aTTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCCTCcgc                                                                                                                    >  1:305270/1‑149 (MQ=255)
                                 taCGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACt                                                                                                                <  1:951393/148‑1 (MQ=255)
                                 cACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACt                                                                                                                >  2:995476/1‑149 (MQ=255)
                                          tAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGt                                                                                                       <  1:160928/149‑1 (MQ=255)
                                            aGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGc                                                                                                     <  1:351840/149‑1 (MQ=255)
                                              ttCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGCTGCAGTTTCATGCGCGACTGCCAGCAGTGCtt                                                                                                   >  2:103544/1‑149 (MQ=255)
                                              ttCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGcgt                                                                                                   >  2:390108/1‑147 (MQ=255)
                                                cAGTGTTTTCTACCCGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCCCCACGTCACCGTTACCGTTGATGATAAGCGATGCAGTTTCATCCGCGACTCCCAGCAGTGCTTCa                                                                                                 >  2:1109485/1‑149 (MQ=255)
                                                  gtgtTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAg                                                                                               <  1:498754/149‑1 (MQ=255)
                                                     ttttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCAACGTTACCGGTGATGTTAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGttt                                                                                            <  2:597712/149‑1 (MQ=255)
                                                      tttCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTg                                                                                           <  1:1097428/149‑1 (MQ=255)
                                                           aaCAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcat                                                                                      >  2:994952/1‑148 (MQ=255)
                                                            aCAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGCCTGCCAGCAGTGCTTCAAGTTTGCgcagc                                                                                     >  1:523654/1‑149 (MQ=255)
                                                             cAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCgcagca                                                                                    >  2:823054/1‑149 (MQ=255)
                                                              aGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGCTGCAGTTTCATCCGCGACTGCCAGCAGTGCGTCAAGTTTGCGCAGCAt                                                                                   >  1:286679/1‑149 (MQ=255)
                                                              aGCGCGCGCGCCGCAGCCTGGGCGTAACGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTTCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCAt                                                                                   <  1:217473/149‑1 (MQ=255)
                                                                cgcgcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCCCCACGTCACCGTTACCGGTGATGATAAGCGATGCTGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATgg                                                                                 >  2:1155331/1‑148 (MQ=255)
                                                                   gcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATTATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTATGCGCAGCATGCGaa                                                                              >  2:202224/1‑148 (MQ=255)
                                                                   gcgcgcCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGAt                                                                              <  1:863448/149‑1 (MQ=255)
                                                                           cAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGc                                                                      <  1:449282/149‑1 (MQ=255)
                                                                               cTGGGCGTAAGGGCCGCCGGAGCCTATAGCAATAAGATAGTTTTCTCTCTGCACCACGTCACCGTTACCGCTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCATTGCTTCAAGTTTGCGAAGCATCAGATCGGAACGCCAGt                                                                  >  2:663354/1‑149 (MQ=255)
                                                                               cTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGt                                                                  <  2:399487/149‑1 (MQ=255)
                                                                                  ggCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTTCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCtt                                                               <  1:820211/149‑1 (MQ=255)
                                                                                    cGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCtttt                                                             >  2:1034817/1‑149 (MQ=255)
                                                                                      gacggggCGCCGGAGCCGTTAGCAATAAGATCGTTTTCTGGCAGCACCTCGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGc                                                           <  1:252031/146‑1 (MQ=255)
                                                                                          ggCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGc                                                       >  2:52417/1‑149 (MQ=255)
                                                                                           gccgccGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCt                                                      >  1:553643/1‑149 (MQ=255)
                                                                                              gccgGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCTCGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCaa                                                   <  1:292578/149‑1 (MQ=255)
                                                                                                cgGAGCAGTTAGCAATAAGATCGTTTTCTGGCTGCACCACGGCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCGGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAAcg                                                 <  2:1061322/149‑1 (MQ=255)
                                                                                                      cgttAGCAATAAGATCGTTTTCTGTGTGCAGCTCGTCAGCGTTACCTGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGAGCAGCATGCGATCGGTTCGCCATTCTTTTGCCAGCTCAACGGCGGCt                                           <  1:51639/146‑1 (MQ=255)
                                                                                                          aGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGa                                       <  1:534266/149‑1 (MQ=255)
                                                                                                           gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGTTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTTCCAGCTCAACGGCGGCTTTGAc                                      <  2:243337/149‑1 (MQ=255)
                                                                                                           gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGc                                                                      <  1:778705/117‑1 (MQ=255)
                                                                                                           gCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGc                                                                      >  2:778705/1‑117 (MQ=255)
                                                                                                                   aTCGTTTTCTGGCTGCCCCACGTCACCGTTACCGGTGATGATAAGCGCAGCAGTTTCAACCGCGACTGCCCGCAGTGCGTAAAGTTTGCGCAGCATGCGGTCGGTACGCCCCGCTTTTGCCAGCTCCCCGGCGGCTTGGACCTGATGGc                              >  2:119876/1‑149 (MQ=255)
                                                                                                                   aTCGTATTCTGGCTGCACCACGTCACCGTTAACGGTGATGATAAGCGATTCAGTTTCATCCGCGACAGCCAGAAGTTCTTAAAGATTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGCCCAGATGGc                              >  2:846883/1‑149 (MQ=255)
                                                                                                                    tCGTTTTCGGGCTGCCCCACGTCCCCGGTACCGGTGATGATAAGCGATGCAGTTGCATCCGCGACTGACCGCAGGGCTTCACGGTTGCGCCGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGcc                             <  1:652803/149‑1 (MQ=255)
                                                                                                                     cGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGccc                            <  2:901536/149‑1 (MQ=255)
                                                                                                                      gTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGATTCATCCGCGACTGCCAGCAGTGCTTAAAGTTTGCGCAGCATGCGATCGGTACGCCACACTTTTGCCAGCTCAACGGCGGCTTTGACCAGAAGGccca                           >  2:1113954/1‑148 (MQ=255)
                                                                                                                      gTTTGCTGGCTGCACCACGTCACCGTTACCGGTGATGATACGAGCGTCAGTTTCAGCCGCGACTGCCAGCAGTGCTTCAATTTTGCGCAGCATGCGATCGTTACGCCAGTCGTTTGCCAGCGCAACGTCGGCTTTGACCAGAGGGCCCt                           <  2:847823/149‑1 (MQ=255)
                                                                                                                       gtgtcaggCTGCCCCACGTCACCGTTACCGGAGATGTTACGCGATGCAGTTCGATCCCCGACGGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGCTCGGTACGCCAGTGTTTTGCCAGCTCCACGGCGGCTTTGTCCAGATGGCCCTg                          <  2:1182512/143‑1 (MQ=255)
                                                                                                                         ttCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGAt                        <  2:1000454/149‑1 (MQ=255)
                                                                                                                              gCTGCACCACGGCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCAtt                   <  2:1158013/149‑1 (MQ=255)
                                                                                                                                 gCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTCCGGCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGAGGTCCCTGATGCATTTcc                >  2:667907/1‑149 (MQ=255)
                                                                                                                                     atcatCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTCCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTCATGCATTTCCAGtt            >  2:450229/5‑149 (MQ=255)
                                                                                                                                        gTCACCGTTACCGGTGATGATTAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTAc         <  2:485719/149‑1 (MQ=255)
                                                                                                                                        gTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGACAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCAAGTCTTTTGCCAGCGCAACGGAGGCTTTGACCAGATGGACCTGATGCAATTCCAGTTTAc         >  1:431721/1‑149 (MQ=255)
                                                                                                                                        atCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTAc         >  1:863209/2‑149 (MQ=255)
                                                                                                                                         tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACg        >  2:681966/1‑149 (MQ=255)
                                                                                                                                         tCACCGGTAACGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGAATACGATCGGTACGCCAGACTTTTGCCAGCACAACGGCGGCTTTGACAATATGGCCCTGATGCATAACCAGTTTACg        >  1:1054211/1‑149 (MQ=255)
                                                                                                                                           aCCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGCCCAGATGGCCCTGATGCATTTCCAGTTTACGtt      >  2:956033/1‑149 (MQ=255)
                                                                                                                                               ttACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCaaa  <  2:553643/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CGCCTGCAATATCCAACGCCTTTTCAGCAATTTCACGGGCGCTAAGTTCAGTGTTTTCTAACAGCGCGCGCGCCGCAGCCTGGGCGTAAGGGCCGCCGGAGCCGATAGCAATAAGATCGTTTTCTGGCTGCACCACGTCACCGTTACCGGTGATGATAAGCGATGCAGTTTCATCCGCGACTGCCAGCAGTGCTTCAAGTTTGCGCAGCATGCGATCGGTACGCCAGTCTTTTGCCAGCTCAACGGCGGCTTTGACCAGATGGCCCTGATGCATTTCCAGTTTACGTTCAAA  >  NC_000913/4121816‑4122107

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: