Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,745,543 A→G Y247C (TAT→TGT)  yiaN → 2,3‑diketo‑L‑gulonate TRAP transporter large permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,5430AG76.9% 175.9 / 44.0 91Y247C (TAT→TGT) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base G (34/36);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 8.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

AAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCTGGCTGA  >  NC_000913/3745409‑3745677
                                                                                                                                      |                                                                                                                                      
aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATgcg                                                                                                                                    <  1:759336/139‑1 (MQ=255)
        cGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            >  2:78780/1‑139 (MQ=255)
           tGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGc                                                                                                                         >  2:1118294/1‑139 (MQ=255)
             cAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  2:699507/139‑1 (MQ=255)
             cAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:1608016/139‑1 (MQ=255)
             cAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTGTACCCCTACCGAGGCCGGGGCCGTTGCTGCGTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:1622104/139‑1 (MQ=255)
             cAGGAAATCTGGCACTCTTTTCTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  2:1057342/139‑1 (MQ=255)
                    tCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTAt                                                                                                                <  2:519764/139‑1 (MQ=255)
                     cTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTAtt                                                                                                               >  1:1435915/1‑139 (MQ=255)
                              ttttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGACGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGaa                                                                                                      >  2:1227730/1‑139 (MQ=255)
                              ttttGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGaa                                                                                                      >  1:66402/1‑139 (MQ=255)
                                   tctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGAc                                                                                                 <  2:66402/139‑1 (MQ=255)
                                        ggTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCg                                                                                            >  2:117987/1‑139 (MQ=255)
                                        ggTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCca                                                                                                                       >  1:582534/1‑112 (MQ=255)
                                        ggTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCca                                                                                                                       <  2:582534/112‑1 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTgcg                                                                                                                                    >  2:267205/1‑90 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTgcg                                                                                                                                    <  1:267205/90‑1 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGc                                                                                                                         >  1:1466499/1‑101 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGc                                                                                                                         <  2:1466499/101‑1 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTAcc                                                                                                           <  1:940533/115‑1 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTAcc                                                                                                           >  2:940533/1‑115 (MQ=255)
                                                 gcgcTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACttt                                                                                   <  2:171358/139‑1 (MQ=255)
                                                    cTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                                                                <  1:1118294/139‑1 (MQ=255)
                                                    cTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                                                                <  1:117987/139‑1 (MQ=255)
                                                     tGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCTGCCTGTTTACCCCTACCGAATCCGGTGCTGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCa                                                                               <  2:1697729/139‑1 (MQ=255)
                                                      gTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            <  1:751342/93‑1 (MQ=255)
                                                      gTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            >  2:751342/1‑93 (MQ=255)
                                                      gTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCAt                                                                              >  2:31098/1‑139 (MQ=255)
                                                               aGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  1:1033405/1‑89 (MQ=255)
                                                               aGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  2:1033405/89‑1 (MQ=255)
                                                                 tcatcaTTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTgg                                                                   >  2:259827/1‑139 (MQ=255)
                                                                 tcatcaTTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTgg                                                                   >  2:1235302/1‑139 (MQ=255)
                                                                   atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:583701/85‑1 (MQ=255)
                                                                   atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  2:583701/1‑85 (MQ=255)
                                                                   atcatTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGgcg                                                                 >  1:1135242/1‑139 (MQ=255)
                                                                      attattGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCgcgg                                                              <  1:1552378/139‑1 (MQ=255)
                                                                          ttGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  2:1791588/1‑78 (MQ=255)
                                                                          ttGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:1791588/78‑1 (MQ=255)
                                                                          ttGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGaa                                                          <  2:1135242/139‑1 (MQ=255)
                                                                           tGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  2:1711756/1‑77 (MQ=255)
                                                                           tGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:1711756/77‑1 (MQ=255)
                                                                               gggTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca                                                     >  1:495434/1‑139 (MQ=255)
                                                                               gggTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAacca                                                     >  1:760286/1‑139 (MQ=255)
                                                                                ggTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCg                                                                                                                          >  2:17083/1‑69 (MQ=255)
                                                                                ggTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCg                                                                                                                          <  1:17083/69‑1 (MQ=255)
                                                                                  tttCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCt                                                  <  1:31098/139‑1 (MQ=255)
                                                                                    tCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCa                                                >  2:1464139/1‑139 (MQ=255)
                                                                                      gCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGt                                              <  2:117902/139‑1 (MQ=255)
                                                                                           ggCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            <  2:185956/56‑1 (MQ=255)
                                                                                           ggCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            >  1:185956/1‑56 (MQ=255)
                                                                                           ggCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGcc                                                                                          >  1:739361/1‑90 (MQ=255)
                                                                                           ggCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGcc                                                                                          <  2:739361/90‑1 (MQ=255)
                                                                                            gCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGc                                                                                                                         >  2:1672788/1‑58 (MQ=255)
                                                                                            gCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGc                                                                                                                         <  1:1672788/58‑1 (MQ=255)
                                                                                                    aCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGTGAAAACCACCTCAg                                               <  2:815077/124‑1 (MQ=255)
                                                                                                    aCCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGTGAAAACCACCTCAg                                               >  1:815077/1‑124 (MQ=255)
                                                                                                      cccTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCg                                                                                            >  2:1344377/1‑77 (MQ=255)
                                                                                                      cccTACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCg                                                                                            <  1:1344377/77‑1 (MQ=255)
                                                                                                      cccTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAac                                                       >  1:1685852/1‑114 (MQ=255)
                                                                                                      cccTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAac                                                       <  2:1685852/114‑1 (MQ=255)
                                                                                                         tACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            <  1:715235/42‑1 (MQ=255)
                                                                                                         tACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGt                                                                                                                            >  2:715235/1‑42 (MQ=255)
                                                                                                         tACCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGcc                           >  2:149976/1‑139 (MQ=255)
                                                                                                          aCCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  2:869494/46‑1 (MQ=255)
                                                                                                          aCCGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  1:869494/1‑46 (MQ=255)
                                                                                                           ccGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       <  1:543487/45‑1 (MQ=255)
                                                                                                           ccGAAGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCca                                                                                                                       >  2:543487/1‑45 (MQ=255)
                                                                                                           ccGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAg                                               >  2:1454897/1‑117 (MQ=255)
                                                                                                           ccGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAg                                               <  1:1454897/117‑1 (MQ=255)
                                                                                                              aaGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGtt                                          <  2:1674056/119‑1 (MQ=255)
                                                                                                              aaGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGtt                                          >  1:1674056/1‑119 (MQ=255)
                                                                                                               aGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                                                                     <  1:1309786/75‑1 (MQ=255)
                                                                                                               aGCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACt                                                                                     >  2:1309786/1‑75 (MQ=255)
                                                                                                                gCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                                                                <  1:1579718/79‑1 (MQ=255)
                                                                                                                gCCGGGGCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGc                                                                                >  2:1579718/1‑79 (MQ=255)
                                                                                                                   ggggCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAg                 >  2:1049607/1‑139 (MQ=255)
                                                                                                                     ggCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGtt               >  1:1440366/1‑139 (MQ=255)
                                                                                                                     ggCCGTTGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGtt               >  1:598277/1‑139 (MQ=255)
                                                                                                                     ggCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGtt               >  2:443365/1‑139 (MQ=255)
                                                                                                                      gCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGttt              >  2:337143/1‑139 (MQ=255)
                                                                                                                          ttGCTGCTTTCTGTGTGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGc          >  1:1757691/1‑139 (MQ=255)
                                                                                                                          ttGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCAt                                                                              >  2:425718/1‑71 (MQ=255)
                                                                                                                          ttGCTGCTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCAt                                                                              <  1:425718/71‑1 (MQ=255)
                                                                                                                           tgctgcTTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGcc         <  2:1440366/139‑1 (MQ=255)
                                                                                                                            tctggcttCTATTCGCTGGTTGTCCGCAAAGTTTTTTACCCTGAAATTAACCTCCCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCt        <  1:1464139/133‑1 (MQ=255)
                                                                                                                            gctgctTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCt        <  1:259827/139‑1 (MQ=255)
                                                                                                                            gctgctTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCt        <  2:495434/139‑1 (MQ=255)
                                                                                                                             ctgctTTCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCTg       <  1:1235302/139‑1 (MQ=255)
                                                                                                                              tgctTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCTgg      <  1:337143/139‑1 (MQ=255)
                                                                                                                                  ttCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATg                                                                             >  1:938519/1‑64 (MQ=255)
                                                                                                                                  ttCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATg                                                                             <  2:938519/64‑1 (MQ=255)
                                                                                                                                  ttCTGTGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCTGGCTGa  <  1:78780/139‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                      
AAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTACCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGGCACAAGTTTCCGCCTGGCTGA  >  NC_000913/3745409‑3745677

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