Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,487,142 G→A 5.1% A118T (GCA→ACA)  yhfK → putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,487,1420GA 5.1% 249.0 / 2.4 78A118T (GCA→ACA) yhfKputative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base A (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CGTCGCTGTTTGCCACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGAACAACCGCTGC  >  NC_000913/3487002‑3487286
                                                                                                                                            |                                                                                                                                                
cGTCGCTGTTTGCCACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCgctg                                                                                                                                          <  2:259326/149‑1 (MQ=255)
        tttGCCACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTcgccgc                                                                                                                                  >  1:2140/1‑149 (MQ=255)
        tttGCCACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTcgcc                                                                                                                                    <  6:355861/147‑1 (MQ=255)
            ccACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAtt                                                                                                                              >  7:378345/1‑149 (MQ=255)
             cACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAttt                                                                                                                             <  1:422192/149‑1 (MQ=255)
              aCCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAtttt                                                                                                                            <  6:313025/149‑1 (MQ=255)
               ccTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAttttt                                                                                                                           <  3:330489/149‑1 (MQ=255)
                             acacAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTTGTACTTGGCGTCACTGCGGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCgc                                                                                                                                            <  7:118778/118‑1 (MQ=255)
                             acacAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTTCTTCATGCATAgc                                                                                                                                            >  8:118778/1‑118 (MQ=255)
                             acacAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTgg                                                                                                              <  1:144741/148‑1 (MQ=255)
                             acacAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGc                                                                                                             <  2:404203/149‑1 (MQ=255)
                               acaGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTa                                                                                                                          >  1:70938/1‑134 (MQ=255)
                               acaGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTa                                                                                                                          <  2:70938/134‑1 (MQ=255)
                               acaGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCgg                                                                                                           <  6:441775/149‑1 (MQ=255)
                               acaGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCgg                                                                                                           <  2:2140/149‑1 (MQ=255)
                                  gCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCTGATTAACGCTGGTACTTGGCGTCACTGCTGATCTGGGGCCATTGCACTCAAAATTGCTTCCTGAATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGAt                                                                                                        <  3:292829/149‑1 (MQ=255)
                                   ctactaCTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATa                                                                                                       <  3:238384/149‑1 (MQ=255)
                                   ctactaCTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCACGCCGCCATTTTTACCCTCAGTTTGGCGGGATa                                                                                                       >  4:112893/1‑149 (MQ=255)
                                    tactacTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATa                                                                                                       >  7:434220/1‑148 (MQ=255)
                                         ctGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTg                                                                                                               <  3:397950/135‑1 (MQ=255)
                                         ctGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTg                                                                                                               >  4:397950/1‑135 (MQ=255)
                                          tGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCg                                                                                                >  4:267193/1‑149 (MQ=255)
                                          tGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCg                                                                                                >  1:745/1‑149 (MQ=255)
                                          tGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGATGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCg                                                                                                >  4:122238/1‑149 (MQ=255)
                                            gCAAAAGATGTTCCCCTGCCCTTTTTGCTGCCCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGt                                                                                              >  6:119519/1‑149 (MQ=255)
                                            gCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTCACGCTGGGACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCg                                                                                                <  1:419375/147‑1 (MQ=255)
                                               aaaGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTg                                                                                           >  2:11352/1‑149 (MQ=255)
                                                    tGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAAc                                                                                      <  1:410801/149‑1 (MQ=255)
                                                         cccTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCTCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTg                                                                                 >  3:27308/1‑149 (MQ=255)
                                                         cccTGCCCTTTTTGCTGACCGGATAAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTg                                                                                 <  1:377771/149‑1 (MQ=255)
                                                          ccTGCCCTTTTTGCTGACCGTATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGc                                                                                >  2:365311/1‑149 (MQ=255)
                                                          ccTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGc                                                                                <  5:17008/149‑1 (MQ=255)
                                                            tGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGTGTCACTGCTGAGCTGGGGCCATTGCACTCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATACCGGTCTGGGAACCGTTGCt                                                                               >  6:220063/1‑148 (MQ=255)
                                                             gcacTTTTTCCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAATTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTc                                                                              >  2:208813/4‑148 (MQ=255)
                                                             gCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCa                                                                             >  2:128225/1‑149 (MQ=255)
                                                                 tttttGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTa                                                                         >  7:236667/1‑149 (MQ=255)
                                                                   tttGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTAt                                                                        <  8:378345/148‑1 (MQ=255)
                                                                   tttGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCACGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGGTGCTCATCTAt                                                                        <  3:112893/148‑1 (MQ=255)
                                                                        tGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTg                                                                  <  8:236667/149‑1 (MQ=255)
                                                                          aCCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTggg                                                                >  4:47333/1‑149 (MQ=255)
                                                                          aCCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTgg                                                                 <  2:745/148‑1 (MQ=255)
                                                                          aCCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTACATCGCTGCTcgc                                                                                                                                     >  5:433626/1‑80 (MQ=255)
                                                                          aCCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTACATCGCTGCTcgc                                                                                                                                     <  6:433626/80‑1 (MQ=255)
                                                                              gATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAttt                                                                                                                             <  4:403490/84‑1 (MQ=255)
                                                                              gATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCAttt                                                                                                                             >  3:403490/1‑84 (MQ=255)
                                                                              gATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCAct                                                            >  5:363376/1‑149 (MQ=255)
                                                                                   aCGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTg                                                       <  6:363376/149‑1 (MQ=255)
                                                                                   aCGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTg                                                       <  5:119519/149‑1 (MQ=255)
                                                                                     gCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGc                                                                                                             <  8:437124/93‑1 (MQ=255)
                                                                                     gCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGc                                                                                                             >  7:437124/1‑93 (MQ=255)
                                                                                      cTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTa                                                    <  3:267193/149‑1 (MQ=255)
                                                                                             ttGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTg                                             >  3:243743/1‑149 (MQ=255)
                                                                                               ggCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGtt                                           >  8:224718/1‑149 (MQ=255)
                                                                                                gCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGttt                                          <  4:243743/149‑1 (MQ=255)
                                                                                                   tCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGc                                                                                                  >  4:225438/1‑90 (MQ=255)
                                                                                                   tCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGc                                                                                                  <  3:225438/90‑1 (MQ=255)
                                                                                                     aCTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTg                                     <  1:365311/149‑1 (MQ=255)
                                                                                                       tgctgAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGt                                   >  5:133458/1‑149 (MQ=255)
                                                                                                         ctgAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGttt                                 <  1:128225/149‑1 (MQ=255)
                                                                                                         ctgAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGttt                                 <  3:122238/149‑1 (MQ=255)
                                                                                                         cggcgcgggggCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTGTTACCCTCAGTTTGGCGGGATACATGCCGGGCTGGGAACCGTTGCTCCTCTATGCGTTGGGCACTCTCTGGTACGGAGTGTTTAACTGGtt                                  <  5:220063/141‑1 (MQ=255)
                                                                                                          tgAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGtttt                                >  1:403350/1‑149 (MQ=255)
                                                                                                           gAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCt                                                                                            >  1:180941/1‑88 (MQ=255)
                                                                                                           gAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCt                                                                                            <  2:180941/88‑1 (MQ=255)
                                                                                                                      aTTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCt                    <  6:133458/149‑1 (MQ=255)
                                                                                                                              cAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGaa            <  1:11352/149‑1 (MQ=255)
                                                                                                                              cAAAATTGCTTCCTACATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGTTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCt                          <  6:188568/135‑1 (MQ=255)
                                                                                                                              cAAAATTGCTTCCTACATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCt                          >  5:188568/1‑135 (MQ=255)
                                                                                                                               aaaaTTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACg                                                  >  8:469101/1‑110 (MQ=255)
                                                                                                                               aaaaTTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACg                                                  <  7:469101/110‑1 (MQ=255)
                                                                                                                                  aTTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTATGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGaacaac        <  3:47333/149‑1 (MQ=255)
                                                                                                                                       ttCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGAACAACCGCTg   >  4:134515/1‑149 (MQ=255)
                                                                                                                                       ttCCTGCATCGCTGCTCGCCGCCATATTTACCCTCAGTTAGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGATTAACTGCTTTTGCTTCTGGATCTGGCGCGAACAACCGCt    >  5:332412/1‑148 (MQ=255)
                                                                                                                                         ccTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGAACAACCGcg    <  5:311872/146‑2 (MQ=255)
                                                                                                                                         ccTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGAACAACCGCt    >  6:311872/1‑146 (MQ=255)
                                                                                                                                         ccTGCATCGCTGCTCGCAGCCATTTTTACCCTCAGTTTGGAGGGATACATGCGGGTATGGGAACCGTTGATCATGTATGCGTTGGTGAATCTCTGGTACGGATTGTTTAACTGGTTTTTGTTCTGGATCTTGCTCTAACAACCGCTgc  <  4:27308/148‑1 (MQ=255)
                                                                                                                                          cTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGa                                                <  6:384703/101‑1 (MQ=255)
                                                                                                                                          cTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGa                                                >  5:384703/1‑101 (MQ=255)
                                                                                                                                            |                                                                                                                                                
CGTCGCTGTTTGCCACCTGTAGCTTGCTGACACAGCTACTACTGGCAAAAGATGTTCCCCTGCCCTTTTTGCTGACCGGATTAACGCTGGTACTTGGCGTCACTGCTGAGCTGGGGCCATTGCACGCAAAATTGCTTCCTGCATCGCTGCTCGCCGCCATTTTTACCCTCAGTTTGGCGGGATACATGCCGGTCTGGGAACCGTTGCTCATCTATGCGTTGGGCACTCTCTGGTACGGATTGTTTAACTGGTTTTGGTTCTGGATCTGGCGCGAACAACCGCTGC  >  NC_000913/3487002‑3487286

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