Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,016:1 +C 17.2% intergenic (+52/+32) yfbP → / ← nuoN TPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,390,0161.C17.2% 579.3 / 130.0 169intergenic (+52/+32)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base . (55/85);  new base C (13/16);  total (68/101)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01

Only 100 of 176 total aligned reads displayed.

CAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTA  >  NC_000913/2389917‑2390115
                                                                                                    |                                                                                                     
 aGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑t                                                                                                      >  2:2373427/1‑100 (MQ=255)
 aGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑t                                                                                                      >  1:140865/1‑100 (MQ=255)
 aGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑t                                                                                                      <  2:1752119/100‑1 (MQ=255)
  gTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tc                                                                                                     <  1:2143041/100‑1 (MQ=255)
    tATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcgg                                                                                                   <  1:995295/100‑1 (MQ=255)
      tGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcggcg                                                                                                 >  1:2687177/1‑100 (MQ=255)
       gTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcggcgg                                                                                                <  2:1298261/100‑1 (MQ=255)
       gTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcggcgg                                                                                                >  2:1187156/1‑100 (MQ=255)
        tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTcggcggt                                                                                               >  1:2778278/1‑99 (MQ=255)
        tttAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTcggcggt                                                                                               <  2:1517176/100‑2 (MQ=255)
           atcttATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑ttt                                                                                          <  1:3278517/96‑1 (MQ=255)
           agagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTCGGCGG‑t‑ttt                                                                                          <  2:2778301/100‑1 (MQ=255)
           agagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACGATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTCGGCGG‑t‑ttt                                                                                          <  2:370688/100‑1 (MQ=255)
            gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTc                                                                                         <  2:2910623/100‑1 (MQ=255)
            gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttc                                                                                        <  2:778357/100‑3 (MQ=25)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTg                                                                                       <  1:575879/100‑1 (MQ=255)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTg                                                                                       <  1:254568/100‑1 (MQ=255)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTg                                                                                       <  2:619976/100‑1 (MQ=255)
              gTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTCGGCGG‑t‑ttttct                                                                                       <  1:844752/100‑4 (MQ=38)
               tATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTATGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTg                                                                                       <  1:439367/99‑1 (MQ=255)
                 ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACt                                                                                    >  2:1753007/1‑100 (MQ=255)
                 ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACt                                                                                    <  1:3073083/100‑1 (MQ=255)
                 ttttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaa                                                                                   <  2:3403868/100‑8 (MQ=21)
                  tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaaa                                                                                  <  1:818564/100‑9 (MQ=18)
                   ttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaaaa                                                                                 <  1:739070/100‑10 (MQ=18)
                   ttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaaaa                                                                                 <  2:331743/100‑10 (MQ=18)
                   ttGGCGATAGTCGTACACAGTTCACGTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaaaa                                                                                 <  2:889770/100‑10 (MQ=18)
                    tGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAACCGCCGA‑CATGCTCGGCGG‑t‑ttttctgaaaaag                                                                                <  2:1822985/100‑11 (MQ=21)
                     ggCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACa                                                                                >  2:3272113/1‑100 (MQ=255)
                       cGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACaaa                                                                              >  2:3869757/1‑100 (MQ=255)
                        gATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAg                                                                             <  2:2159516/100‑1 (MQ=255)
                         aTAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGt                                                                            <  1:2494966/100‑1 (MQ=255)
                         aTAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGt                                                                            <  2:2224136/100‑1 (MQ=255)
                             tCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTTATTTTATCCGCCGAGTATG‑TCGGCGG‑C‑TTTCTGACTTACAAAATAACa                                                                        >  1:3027005/1‑100 (MQ=255)
                             tCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACa                                                                        <  2:459855/100‑1 (MQ=255)
                              cGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAg                                                                       <  2:1952451/100‑1 (MQ=255)
                              cGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAg                                                                       >  1:2287232/1‑100 (MQ=255)
                              cGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAg                                                                       <  1:1775399/100‑1 (MQ=255)
                                  cacaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTa                                                                   <  1:1629691/100‑1 (MQ=255)
                                  cacaGTTCACTTAACGCTCTGCGTCTATCTACTAGCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCTGCGG‑C‑TTTCTGACTTACTAAGTAACAGATTa                                                                   <  2:2614061/100‑1 (MQ=255)
                                   acaGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTAc                                                                  <  1:2976296/100‑1 (MQ=255)
                                    caGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACa                                                                 >  1:1500586/1‑100 (MQ=255)
                                         cACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGc                                                            >  2:136407/1‑100 (MQ=255)
                                          aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCg                                                           >  2:3833808/1‑100 (MQ=255)
                                          aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCg                                                           <  1:245316/100‑1 (MQ=255)
                                          aCTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCg                                                           <  2:140725/100‑1 (MQ=255)
                                             tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                        >  2:54023/1‑100 (MQ=255)
                                             tAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                        >  2:1915584/1‑100 (MQ=255)
                                                 gCTCTGCGTCTATCTACTATCACTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGc                                                    <  1:811728/100‑1 (MQ=255)
                                                   tctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCa                                                  <  1:2683292/100‑1 (MQ=255)
                                                   tctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCa                                                  >  1:397926/1‑100 (MQ=255)
                                                   tctGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGAATTACAAACTAACAGAATACATCAGCGGCATTGCCa                                                  >  1:1085556/1‑100 (MQ=255)
                                                     tGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaaa                                                <  2:1494488/100‑1 (MQ=255)
                                                     tGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaaa                                                <  2:642866/100‑1 (MQ=255)
                                                         tctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGca                                            <  1:2517247/100‑1 (MQ=255)
                                                            atctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACaa                                         >  2:1660484/1‑100 (MQ=255)
                                                             tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAAt                                        <  2:3569443/100‑1 (MQ=255)
                                                             tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAAt                                        <  1:2120003/100‑1 (MQ=255)
                                                             tctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAAt                                        <  2:352355/100‑1 (MQ=255)
                                                              ctactaTCCCTTCTTTTTCATTTTATCTTCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATg                                       >  1:1188516/1‑100 (MQ=255)
                                                              ctactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATg                                       >  2:1008503/1‑100 (MQ=255)
                                                               tactaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                      <  1:1483374/100‑1 (MQ=255)
                                                                actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCt                                     <  2:3342364/100‑1 (MQ=255)
                                                                   aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAAt                                  <  2:3286307/100‑1 (MQ=255)
                                                                     cccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCa                                <  1:3619043/100‑1 (MQ=255)
                                                                      ccTGCTTGTTCAGAACAAGCGCCTAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTTACAGATTAAATCAGCGGCATTGCCAAACGCAGAATGCTAATCa                                <  1:3819869/99‑1 (MQ=255)
                                                                        ttcttTTTCAGAAAAACTGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCg                             <  1:89374/100‑1 (MQ=255)
                                                                        ttcttTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCg                             <  2:3150136/100‑1 (MQ=255)
                                                                         tcttTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCgg                            >  1:537941/1‑100 (MQ=255)
                                                                            ttttCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtg                         >  2:3164206/1‑100 (MQ=255)
                                                                             tttCAGACAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTCCATCAGCTGCATTGCCAAACGCACAATGCTAATCAGCgg                            >  1:3291940/1‑96 (MQ=255)
                                                                             tttCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  1:1348732/1‑100 (MQ=255)
                                                                              agtcagaaagccgccgA‑CATGCTCGGCGG‑TTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  2:2074872/11‑100 (MQ=21)
                                                                              agtcagaaagCCGCCGA‑CATGCTCGGCGGTT‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                        >  2:1144159/7‑100 (MQ=21)
                                                                               tCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      <  2:1989958/100‑1 (MQ=255)
                                                                               tCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      <  2:1312256/100‑1 (MQ=255)
                                                                                tcagaaagCCGCCGA‑CATGCTCGGCGGTT‑TTTCTGACTTACAAAGTAGCAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                      >  1:2898594/5‑100 (MQ=21)
                                                                                cAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  2:3684114/1‑100 (MQ=255)
                                                                                 cagaaagccgccgA‑CATGCTCGGCGG‑TTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGc                     >  2:98527/8‑100 (MQ=25)
                                                                                 aGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGcc                    >  2:257812/1‑100 (MQ=255)
                                                                                  gAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCa                   <  2:1465028/100‑1 (MQ=255)
                                                                                   gaaagCCGCCGA‑CATGCTCGGCGGTT‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCa                   >  2:211215/2‑100 (MQ=25)
                                                                                   aaaaaCCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCAt                  >  1:2994026/1‑100 (MQ=255)
                                                                                    aaaaCCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                 >  1:2015718/1‑100 (MQ=255)
                                                                                    aaaaCCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATa                 >  2:3595561/1‑100 (MQ=255)
                                                                                       gccgccgA‑CATGCTCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac              >  2:1281569/2‑100 (MQ=255)
                                                                                       aCCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac              >  2:3473970/1‑100 (MQ=255)
                                                                                          gccgTGCATG‑TCGGCGT‑G‑TTTCTGACTTCCCAAGTAACAGATTACATCAGCGGCATTGCCAAACGCAAAATGCTAATCAGCGGTTGTGGCCATACACCCa           <  1:1417943/100‑1 (MQ=255)
                                                                                          gccgAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCa           >  1:628692/1‑100 (MQ=255)
                                                                                          gccgAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCa           >  1:2601712/1‑100 (MQ=255)
                                                                                            cgAGCATG‑TCGGCGG‑C‑TTTCTGACTCATAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc         >  1:1435002/1‑100 (MQ=255)
                                                                                             gAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCa        <  1:3724830/100‑1 (MQ=255)
                                                                                             gAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCa        <  2:3457660/100‑1 (MQ=255)
                                                                                              aGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAc       <  1:3799021/100‑1 (MQ=255)
                                                                                              aGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAc       <  2:3875745/100‑1 (MQ=255)
                                                                                               gCATG‑TCGGCGT‑C‑TGGCTGACTTAAATGGGGAGGGATTACGTCAGCGGCAGGGCCAAACGCACAATGCTAATCAGCGGGTGTGGCGATACACCCAGCAcc      <  1:3115577/100‑1 (MQ=255)
                                                                                               gCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGc         <  2:2422795/97‑1 (MQ=255)
                                                                                               gCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAcc      <  2:2753929/100‑1 (MQ=255)
                                                                                               gCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAcc      >  2:3421716/1‑100 (MQ=255)
                                                                                                   g‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTa  <  1:1150374/100‑1 (MQ=255)
                                                                                                    |                                                                                                     
CAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑C‑TTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTA  >  NC_000913/2389917‑2390115

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