Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,966 C→T 6.3% A404V (GCG→GTG)  atoE → short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,325,9660CT 6.3% 205.2 / 3.5 63A404V (GCG→GTG) atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCGCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAATATAAAAT  >  NC_000913/2325819‑2326110
                                                                                                                                                   |                                                                                                                                                
tCGCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGgc                                                                                                                                                  >  2:276657/1‑148 (MQ=255)
tCGCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGgc                                                                                                                                                  >  8:221165/1‑148 (MQ=255)
tCGCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGgc                                                                                                                                                  >  7:254604/1‑148 (MQ=255)
  gCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGgcgc                                                                                                                                                <  3:239307/148‑1 (MQ=255)
  gCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGgcg                                                                                                                                                 <  1:66841/147‑1 (MQ=255)
       tCCGTCTGGCGGCGGTCACGGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCa                                                                                                                                           >  5:140452/1‑148 (MQ=255)
        ccGGCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGAAACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCa                                                                                                                                           <  6:100655/147‑1 (MQ=255)
         cGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGc                                                                                                                                         <  6:140452/148‑1 (MQ=255)
         agTCTGGCGGCGGTGACTGGGTTATTCAGGGTGCTTTGGTGATTGCCGGAGGTCAGGGGCTTGGCGCTGTTCTCGGTAAATCTGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACATCCATTCTGGGCGCTGCCAgg                                                                                                                                         <  4:240277/147‑2 (MQ=255)
          gTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCa                                                                                                                                        <  4:78134/148‑1 (MQ=255)
           tCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCAc                                                                                                                                       >  1:110906/1‑148 (MQ=255)
            cTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGTGCTGCCAGCACt                                                                                                                                      >  3:99112/1‑148 (MQ=255)
                         cTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGa                                                                                                                         >  1:57212/1‑148 (MQ=255)
                              ttATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCgg                                                                                                                    <  8:117168/148‑1 (MQ=255)
                                 ttCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGcc                                                                                                                            >  6:74797/1‑137 (MQ=255)
                                 ttCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGcc                                                                                                                            <  5:74797/137‑1 (MQ=255)
                                 ttCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGgtg                                                                                                                  <  7:68374/147‑1 (MQ=255)
                                   cAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTc                                                                                                                <  7:221165/147‑1 (MQ=255)
                                            ttCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGAcatca                                                                                                      >  1:165553/1‑148 (MQ=255)
                                             tCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGAcatca                                                                                                      <  1:178341/147‑1 (MQ=255)
                                              cGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGAcatcat                                                                                                     <  5:3173/147‑1 (MQ=255)
                                              cGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGAcatcat                                                                                                     <  5:112407/147‑1 (MQ=255)
                                               gTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATg                                                                                                    <  2:165553/147‑1 (MQ=255)
                                                 gATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGc                                                                                                 >  2:144362/1‑148 (MQ=255)
                                                 gATACCCGCAGCCAAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGc                                                                                                 >  2:143894/1‑148 (MQ=255)
                                                    aCCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTAc                                                                                              >  3:75213/1‑148 (MQ=255)
                                                     cccGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTAc                                                                                              <  4:253540/147‑1 (MQ=255)
                                                     cccGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACt                                                                                             <  8:254604/148‑1 (MQ=255)
                                                          aGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATc                                                                                        <  4:75213/148‑1 (MQ=255)
                                                            cccAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCAc                                                                                      >  8:66470/1‑148 (MQ=255)
                                                              cAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTg                                                                                    <  1:144362/148‑1 (MQ=255)
                                                               aGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGc                                                                                   >  8:62366/1‑148 (MQ=255)
                                                                    ctgggcgctgATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGTATGAACATGGCACAACCCTTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGc                                                                              <  4:12009/148‑1 (MQ=255)
                                                                                     tAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCgg                                                                    >  4:199407/1‑141 (MQ=255)
                                                                                     tAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCgg                                                                    <  3:199407/141‑1 (MQ=255)
                                                                                     tAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCAtt                                                             >  1:302098/1‑148 (MQ=255)
                                                                                        aTCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTc                                                          >  4:207872/1‑148 (MQ=255)
                                                                                        aTCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTc                                                          >  1:272128/1‑148 (MQ=255)
                                                                                            gTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCa                                                      <  2:272128/148‑1 (MQ=255)
                                                                                             tAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCAt                                                     <  2:302098/148‑1 (MQ=255)
                                                                                               aTGGCGTTCGACTACGGCGAGCAATGGATGAACATTGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCAtt                                                    <  5:124019/147‑1 (MQ=255)
                                                                                               aTGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCAtt                                                    <  1:276657/147‑1 (MQ=255)
                                                                                                tGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTgg                                                  >  1:250420/1‑148 (MQ=255)
                                                                                                tGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTgg                                                  <  7:62366/148‑1 (MQ=255)
                                                                                                tGGAGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTCCTGCATCACTGACCTGCTCTTCTCCGGTGTCATTTTCGTCATTgg                                                  >  2:239149/1‑148 (MQ=255)
                                                                                                 ggCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGTGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCAttt                                                            <  6:176567/137‑1 (MQ=255)
                                                                                                 ggCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGTGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCAttt                                                            >  5:176567/1‑137 (MQ=255)
                                                                                                 ggCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGt                                                 >  3:93062/1‑148 (MQ=255)
                                                                                                  gCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCTCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGATACTGCATCACTGCCATGCTCTTCTCCGGTTTCATTTTCGTCATTTGtt                                                >  4:133895/1‑148 (MQ=255)
                                                                                                  gCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGt                                                 <  2:110906/147‑1 (MQ=255)
                                                                                                    gATCGCCTACGGCGAGAAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTa                                              <  3:133895/148‑1 (MQ=255)
                                                                                                               gCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTg                                   <  7:66470/148‑1 (MQ=255)
                                                                                                                              aCATGGCACAACCATTCTGGGTGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACaa                    <  4:99112/148‑1 (MQ=255)
                                                                                                                                  ggCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCCGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACaaa                   >  8:322655/1‑145 (MQ=255)
                                                                                                                                   gCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCCGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGa               <  7:322655/148‑1 (MQ=255)
                                                                                                                                    cacaACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAg                >  5:173126/1‑146 (MQ=255)
                                                                                                                                      caACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATGTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAg             >  1:63781/1‑147 (MQ=255)
                                                                                                                                       aaCCATTCTGGGCGCTGCCAGCACTGGTAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGa           <  4:277506/148‑1 (MQ=255)
                                                                                                                                       aaCCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCGCCCCTACTAACAGAAgg            >  6:159563/1‑147 (MQ=255)
                                                                                                                                        aCCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGaa          >  2:36785/1‑148 (MQ=255)
                                                                                                                                          cATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAata        <  1:36785/148‑1 (MQ=255)
                                                                                                                                          cATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAata        >  8:206416/1‑148 (MQ=255)
                                                                                                                                          cATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAata        >  5:125806/1‑148 (MQ=255)
                                                                                                                                           aTTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTAGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAatat       <  7:206416/148‑1 (MQ=255)
                                                                                                                                             tCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGtt                                                >  3:212652/1‑105 (MQ=255)
                                                                                                                                             tCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGtt                                                <  4:212652/105‑1 (MQ=255)
                                                                                                                                                gggCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAATATAAAAt  >  2:71166/1‑148 (MQ=255)
                                                                                                                                                gggCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAATATAAAAt  >  1:292762/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TCGCCGTTCCGTCTGGCGGCGGTCACTGGGTTATTCAGGGACCTTTCGTGATACCCGCAGCCCAGGCGCTGGGCGCTGATCTCGGTAAATCGGTAATGGCGATCGCCTACGGCGAGCAATGGATGAACATGGCACAACCATTCTGGGCGCTGCCAGCACTGGCAATCGCCGGACTCGGTGTCCGCGACATCATGGGCTACTGCATCACTGCCCTGCTCTTCTCCGGTGTCATTTTCGTCATTGGTTTAACGCTGTTCTGACGGCACCCCTACAAACAGAAGGAATATAAAAT  >  NC_000913/2325819‑2326110

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: