Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,543 T→A 5.3% V83D (GTC→GAC)  ycdZ → DUF1097 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,5430TA 5.3% 325.1 / 7.2 95V83D (GTC→GAC) ycdZDUF1097 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/57);  new base A (3/2);  total (36/59)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

CTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGC  >  NC_000913/1100395‑1100689
                                                                                                                                                    |                                                                                                                                                  
cagggCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTACGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGt                                                                                                                                                    <  4:900160/147‑1 (MQ=255)
cTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGt                                                                                                                                                    <  2:864124/149‑1 (MQ=255)
   ggCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGACAt                                                                                                                                                 >  1:557741/1‑149 (MQ=255)
         aCGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCgg                                                                                                                                           <  2:553291/149‑1 (MQ=255)
         aCGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGACATAACCga                                                                                                                                           <  2:557741/149‑2 (MQ=255)
            gCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGCGGGGTGTGGGCGATGGGCATTATTTACGGTAGTTCACTGGCACCCCCTCTGGAAATTCTCGGTTATGTCATAACCGGTAtt                                                                                                                                       >  1:489916/1‑149 (MQ=255)
             ccTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTAtt                                                                                                                                       >  4:246573/1‑149 (MQ=255)
             ccTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTAGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTAtt                                                                                                                                       <  4:147021/149‑1 (MQ=255)
                aCTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTc                                                                                                                                    <  3:132892/149‑1 (MQ=255)
                aCTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTc                                                                                                                                    <  3:560710/149‑1 (MQ=255)
                  tttGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTcgc                                                                                                                                  >  3:626356/1‑149 (MQ=255)
                  tttGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTcgc                                                                                                                                  <  1:541133/149‑1 (MQ=255)
                  tttGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTcgc                                                                                                                                  <  2:674315/149‑1 (MQ=255)
                    tgcctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGt                                                                                                                                <  1:851876/149‑1 (MQ=255)
                     gcctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGtt                                                                                                                               <  3:246573/149‑1 (MQ=255)
                     gcctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGtt                                                                                                                               <  1:754701/149‑1 (MQ=255)
                      cctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGt                                                                                                                                     <  1:24391/142‑1 (MQ=255)
                      cctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGt                                                                                                                                     >  2:24391/1‑142 (MQ=255)
                       ctgccGGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTAtt                                                                                                                                       <  4:421165/139‑1 (MQ=255)
                       ctgccCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTAtt                                                                                                                                       >  3:421165/1‑139 (MQ=255)
                          cccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          <  2:265518/149‑1 (MQ=255)
                          cccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          <  4:163997/149‑1 (MQ=255)
                          cccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          <  4:406323/149‑1 (MQ=255)
                           ccGCAAGGTGGGCTGAAAGTGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAt                                                                                                                         <  4:717480/149‑1 (MQ=255)
                           ccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAt                                                                                                                         <  1:574758/149‑1 (MQ=255)
                           ccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAt                                                                                                                         <  1:528815/149‑1 (MQ=255)
                           ccGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAt                                                                                                                         >  1:486578/1‑149 (MQ=255)
                             gCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAtgt                                                                                                                       <  3:194990/149‑1 (MQ=255)
                             gCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAtgt                                                                                                                       <  1:430901/149‑1 (MQ=255)
                                     ggCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAg                                                                                                               <  1:416755/149‑1 (MQ=255)
                                     ggCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAg                                                                                                               <  2:545441/149‑1 (MQ=255)
                                      gCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGt                                                                                                                                          <  3:122015/121‑1 (MQ=255)
                                      gCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGt                                                                                                                                          >  4:122015/1‑121 (MQ=255)
                                         gAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCa                                                                                                           <  4:261717/149‑1 (MQ=255)
                                            aGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAAc                                                                                                        >  3:330377/1‑149 (MQ=255)
                                                 tGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAgctg                                                                                                   <  2:799993/149‑1 (MQ=255)
                                                   gCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAgc                                                                                                     <  4:720436/145‑1 (MQ=255)
                                                   gCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAgc                                                                                                     >  3:720436/1‑145 (MQ=255)
                                                     gATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAgctgctgc                                                                                               >  3:333895/1‑149 (MQ=255)
                                                        ctcCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCttt                                                                                            <  1:193954/149‑1 (MQ=255)
                                                         tcCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTc                                                                                           >  3:377055/1‑149 (MQ=255)
                                                         tcCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTc                                                                                           >  2:730516/1‑149 (MQ=255)
                                                          ccGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCa                                                                                          <  1:791784/149‑1 (MQ=255)
                                                           cGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCAt                                                                                         <  1:326254/149‑1 (MQ=255)
                                                             cTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCAttt                                                                                       <  1:217915/149‑1 (MQ=255)
                                                             cTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCAttt                                                                                       <  3:485717/149‑1 (MQ=255)
                                                             cTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCAttt                                                                                       <  2:150722/149‑1 (MQ=255)
                                                              tGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTg                                                                                      >  1:804007/1‑149 (MQ=255)
                                                              tGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTg                                                                                      >  4:474565/1‑149 (MQ=255)
                                                                 aaCCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTc                                                                                   >  4:124362/1‑149 (MQ=255)
                                                                  aCCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTcc                                                                                  >  4:43724/1‑149 (MQ=255)
                                                                   cccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCg                                                                                 >  3:169368/1‑149 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          >  3:680026/1‑107 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          >  1:1003786/1‑107 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          <  4:680026/107‑1 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGa                                                                                                                          <  2:1003786/107‑1 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCgg                                                                                >  1:804260/1‑149 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCgg                                                                                <  2:371826/149‑1 (MQ=255)
                                                                    ccTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCggg                                                                               <  3:560474/149‑1 (MQ=255)
                                                                     cTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCggg                                                                               <  1:182650/149‑1 (MQ=255)
                                                                     cTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCggg                                                                               <  3:446881/149‑1 (MQ=255)
                                                                      tGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCgggg                                                                              >  1:706025/1‑149 (MQ=255)
                                                                       gCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGa                                                                             >  4:224487/1‑149 (MQ=255)
                                                                          tAGTGGCGTGTTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGt                                                                          <  3:387430/149‑1 (MQ=255)
                                                                            gTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGttt                                                                        >  2:778388/1‑149 (MQ=255)
                                                                               gCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTtata                                                                     <  2:804007/149‑1 (MQ=255)
                                                                               gCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTtata                                                                     <  4:333895/149‑1 (MQ=255)
                                                                                  tggtgTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGgc                                                                  <  1:778388/149‑1 (MQ=255)
                                                                                      gtgGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGTCGCAt                                                              <  2:397215/149‑1 (MQ=255)
                                                                                       tgGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCgg                                                                                <  1:1063039/130‑1 (MQ=255)
                                                                                       tgGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCgg                                                                                >  2:1063039/1‑130 (MQ=255)
                                                                                        gggCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATgc                                                            <  3:158418/149‑1 (MQ=255)
                                                                                         ggCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAtgt                                                                                                                       >  4:173252/1‑89 (MQ=255)
                                                                                         ggCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGAtgt                                                                                                                       <  3:173252/89‑1 (MQ=255)
                                                                                             aTGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGAc                                                       <  1:1001558/149‑1 (MQ=255)
                                                                                              tGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACg                                                      >  3:498847/1‑149 (MQ=255)
                                                                                                gTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGtt                                                    >  4:758433/1‑149 (MQ=255)
                                                                                                 tCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGttt                                                   <  3:43724/149‑1 (MQ=255)
                                                                                                    ttattTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGcc                                                <  1:783756/149‑1 (MQ=255)
                                                                                                           cGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAgg                                         <  4:498847/149‑1 (MQ=255)
                                                                                                           cGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAgg                                         <  2:681536/149‑1 (MQ=255)
                                                                                                              tAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTATCATTTGTTCCGGGGACGTTTATAGACGCATGCGCGACGTTTGCCGGGCAGGGTg                                      >  1:707148/1‑149 (MQ=255)
                                                                                                               aGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGa                                     <  2:486578/149‑1 (MQ=255)
                                                                                                               aGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGa                                     <  4:621961/149‑1 (MQ=255)
                                                                                                                    aCTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCggg                                                                               >  4:310339/1‑102 (MQ=255)
                                                                                                                    aCTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCggg                                                                               <  3:310339/102‑1 (MQ=255)
                                                                                                                        gCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACt                            >  1:708233/1‑149 (MQ=255)
                                                                                                                             acaTCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGgtgt                       <  1:180695/149‑1 (MQ=255)
                                                                                                                                    gAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTc                <  1:888598/149‑1 (MQ=255)
                                                                                                                                    gAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTc                >  4:795033/1‑149 (MQ=255)
                                                                                                                                        ttCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTg            >  2:70146/1‑149 (MQ=255)
                                                                                                                                         tctcGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTgg           >  1:683192/1‑149 (MQ=255)
                                                                                                                                           ccGGTTATGACATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGgcg         >  4:746074/2‑149 (MQ=255)
                                                                                                                                            cGGTTATGACATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTAc                    <  3:757428/137‑1 (MQ=255)
                                                                                                                                            cGGTTATGACATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTAc                    >  4:757428/1‑137 (MQ=255)
                                                                                                                                                tATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTgg    <  1:283297/149‑1 (MQ=255)
                                                                                                                                                  tGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAATCAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGc  <  3:188941/149‑1 (MQ=255)
                                                                                                                                                    tCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGgt                         <  1:655324/124‑1 (MQ=255)
                                                                                                                                                    tCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGgt                         >  2:655324/1‑124 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CTGGGCTGTACGGCCTACTTTGCCTGCCCGCAAGGTGGGCTGAAAGGGCTGGCGATCTCCGCTGCAACCCTGCTTAGTGGCGTGGTGTGGGCGATGGTCATTATTTACGGTAGTGCACTGGCACCACATCTGGAAATTCTCGGTTATGTCATAACCGGTATTGTCGCGTTTCTGATGTGTATTCAGGCCAAACAGCTGCTGCTTTCATTTGTTCCGGGGACGTTTATAGGCGCATGCGCGACGTTTGCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGC  >  NC_000913/1100395‑1100689

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