Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,390,016:1 +C 5.1% 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.C 5.1% 321.8 / 12.5 79intergenic (+52/+32)yfbP/nuoNTPR‑like repeats‑containing protein/NADH:ubiquinone oxidoreductase, membrane subunit N
Reads supporting (aligned to +/- strand):  ref base . (37/38);  new base C (2/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

GGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTACTGCC  >  NC_000913/2389871‑2390163
                                                                                                                                                  |                                                                                                                                                    
ggAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcgg                                                                                                                                                  <  1:896863/150‑1 (MQ=255)
ggAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTACACTTAACGCTCTGCGTATATCTAATATACCTTCTTTTTCAGAAAAACCGCCGAGCAAG‑Tcgg                                                                                                                                                  >  2:195648/1‑150 (MQ=255)
 gAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑Tcggc                                                                                                                                                 <  1:375057/150‑1 (MQ=255)
    cGCTATGATTTTAAAATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGA‑CATGCTcggcggt                                                                                                                                              <  2:1177469/150‑2 (MQ=255)
    agCTATTATTTTAATCTCGATGAATGGCGGGCAAAATTCAGGCAGGTATTTTTTGAGTAGGTTGGCGATAGTCCTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGTGCATG‑TCggcgg‑c                                                                                                                                             <  2:801153/149‑1 (MQ=255)
      cTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑Ctt                                                                                                                                           <  2:1198148/150‑1 (MQ=255)
      cTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑Ctt                                                                                                                                           <  1:1314909/150‑1 (MQ=255)
        aTGATTTTAATATAGATGAATGGCGTGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGGCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTc                                                                                                                                         <  1:109458/150‑1 (MQ=255)
         tGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCGTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCt                                                                                                                                        >  1:806522/1‑150 (MQ=255)
         tGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCt                                                                                                                                        >  2:484953/1‑150 (MQ=255)
           aTTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAG‑AAAGCCGCCGA‑CATGCTCGGCGG‑tttttctga                                                                                                                                     <  1:860736/150‑6 (MQ=38)
                     cGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAACACCGCCGCGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGt                                                                                                                            >  1:501202/1‑150 (MQ=255)
                        tGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAAc                                                                                                                         <  2:686724/150‑1 (MQ=255)
                         gAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACa                                                                                                                        <  2:955795/150‑1 (MQ=255)
                              gctggCATAATTGACGCTGTTATGTTTAGAGTATTTTTGCGATAGTCGTACACAGTTCACTTAATGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGAATG‑TGGGCGG‑CTTTCTGACTTACAAAGTAACAGATTa                                                                                                                   <  2:574352/147‑1 (MQ=255)
                               cGGGCAAAATTCAGGCAGTTATGTGTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCAAG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATGAc                                                                                                                  >  2:207044/1‑150 (MQ=255)
                                 ggCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAACCCGCCGAGCATG‑GCGGAGG‑CTTTCTGACTTACAAAGTACCAGATTACAt                                                                                                                >  2:490078/1‑150 (MQ=255)
                                 ggCAAAATTCAGGCAGTTAGGTTTAGAGTATTTTGGCGATAGTCGTGCACAGGTCACTTAACGCTCTGCGTCTATCTAGTAGCCCTTCTTTTTCAGAAAAACCGCGGGGCTTG‑TCGGGGG‑CTTTCTGGCTGAGAAAGGAACGGGTTAcct                                                                                                                >  2:851578/1‑148 (MQ=255)
                                  gCAATATTAAGGCAGTGATGTTTAGAGTATTTTGGAGATAGTCGTACACATTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATc                                                                                                               <  2:1170106/150‑1 (MQ=255)
                                  gCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATc                                                                                                               <  2:1051371/150‑1 (MQ=255)
                                         tCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                                                                        <  2:856226/150‑1 (MQ=255)
                                         tCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCa                                                                                                        <  2:194909/150‑1 (MQ=255)
                                             gCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTTACAGATTACATCAGCGGCATTGc                                                                                                    <  1:195887/150‑1 (MQ=255)
                                             gCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCAGAGCATA‑TCGGCGG‑CTTTCTGACATAAAAAGTAACAGATTACATCAGCGGCAATGc                                                                                                    >  2:1189227/1‑150 (MQ=255)
                                             gCAGCTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGc                                                                                                    <  1:451158/150‑1 (MQ=255)
                                               aGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCa                                                                                                  >  1:772052/1‑150 (MQ=255)
                                                gTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCaa                                                                                                 >  2:775656/1‑150 (MQ=255)
                                                   aTGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACg                                                                                              <  1:1160301/150‑1 (MQ=255)
                                                    tGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGc                                                                                             >  1:1161203/1‑150 (MQ=255)
                                                        tAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACaa                                                                                         >  2:1239757/1‑150 (MQ=255)
                                                          gagTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATg                                                                                       <  1:661004/150‑1 (MQ=255)
                                                           agTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGc                                                                                      <  1:261585/150‑1 (MQ=255)
                                                              aTTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTaa                                                                                   >  2:554907/1‑150 (MQ=255)
                                                              aTTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATATACTATCCATTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGAACAATGCTaa                                                                                   >  1:1245596/1‑150 (MQ=255)
                                                               ttttGGCGATAGTCGTACACAGTTCACTTAAGGCGCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAAt                                                                                  <  2:483583/150‑1 (MQ=255)
                                                                tttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATc                                                                                 >  2:286577/1‑150 (MQ=255)
                                                                tttGGCGATAGTCGTACACAGTTCACATAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCTCCGAGCATG‑TCGGCGG‑CATTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATc                                                                                 <  1:286375/150‑1 (MQ=255)
                                                                 ttGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCa                                                                                >  1:1026407/1‑150 (MQ=255)
                                                                  tGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCAGG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCCGCGGCATTGCCAAACGCACAATGCTAATCAg                                                                               >  2:798596/1‑150 (MQ=255)
                                                                    gCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTAAATCAGCGGCATTGCCAAACGCACGATGCTAATCAGCg                                                                             >  1:194166/1‑150 (MQ=255)
                                                                      gATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGt                                                                           <  2:1267351/150‑1 (MQ=255)
                                                                         aGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTtgt                                                                        >  1:1189901/1‑150 (MQ=255)
                                                                           tCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTgg                                                                      <  2:1037936/150‑1 (MQ=255)
                                                                                  caGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacc                                                               >  1:524793/1‑149 (MQ=255)
                                                                                  caGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCAAG‑TCGGCGG‑CATTCTGACTTACAAAGTAACAGATTACATCAGCGTCATTGCCAAACGCACAATACTAATCAGCGGTTGTGGCCAAaca                                                               >  2:62198/1‑150 (MQ=255)
                                                                                   aGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac                                                              <  1:361766/150‑1 (MQ=255)
                                                                                   aGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATacac                                                              >  1:414372/1‑150 (MQ=255)
                                                                                      tCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCa                                                           <  2:1012920/150‑1 (MQ=255)
                                                                                       cACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAAAAGATTACATCAGCGGCATTGCCAAACGCACAATTCTAATCAGCGGTTGTGGCCATACACCCAg                                                          >  1:753961/1‑150 (MQ=255)
                                                                                        agttAACTCTCTGGGTGTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAAAAGATTACATCAGCGGCATTGCCAAACGCACAATTCTAATCAGCGGTTGTGGCCATACACCCAgg                                                         <  2:543046/148‑2 (MQ=255)
                                                                                           tAACGCTCTGCGTCTATCTACTATCCCTTCTATTTCAGAAAAACCGCCGAGCATG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCAcc                                                      <  2:124421/150‑1 (MQ=255)
                                                                                              cGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCTGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGt                                                   <  1:712664/150‑1 (MQ=255)
                                                                                                        ctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACATAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATCCACCCAGCACCAGTACCCACAGTg                                         >  2:550200/1‑150 (MQ=255)
                                                                                                        ctatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTg                                         <  2:438653/150‑1 (MQ=255)
                                                                                                         tatctaCTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACAACAGCGGCATTGCCAAACGCACAATGCTAATCAGCGTTTGTGGCAATACACCCAGCACCAGTACCAACAGTGc                                        >  1:619253/1‑150 (MQ=255)
                                                                                                              actaTCCCTTCTTTTTCAGAACAACCGCCGAGCATG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAgaga                                   <  2:459139/150‑1 (MQ=255)
                                                                                                              actaTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCTG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGAGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGAGGCCATACACCCAGAACCAGTACAAACAGAGCAgaga                                   >  1:198386/1‑150 (MQ=255)
                                                                                                                 aTCCCTTCTTTTTCAGAAAAACCGCCGCGCATG‑GCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCa                                <  1:697583/150‑1 (MQ=255)
                                                                                                                 aTCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCa                                >  2:456688/1‑150 (MQ=255)
                                                                                                                  tCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAg                               >  1:26524/1‑150 (MQ=255)
                                                                                                                   cccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGc                              >  1:579153/1‑150 (MQ=255)
                                                                                                                    ccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGTTTACATCAGCGGCATTGCCAATCGCACAATGCTAATCAGCGGTTGTGGCCATACAGCCAGCACCAGTACCAACAGTGCAGAGATCAGca                             >  1:63911/1‑150 (MQ=255)
                                                                                                                    ccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCTTTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGca                             <  2:392164/150‑1 (MQ=255)
                                                                                                                    ccTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGca                             <  2:1210263/150‑1 (MQ=255)
                                                                                                                      ttcttTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGcaca                           <  1:574806/150‑1 (MQ=255)
                                                                                                                            agtcagaaagCCGCCGA‑CATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGAt                      >  1:677833/7‑150 (MQ=38)
                                                                                                                            agtcagaaagCCGCCGA‑CATGCTCGGCGGTTTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGAt                      >  1:783312/7‑150 (MQ=38)
                                                                                                                             tCAGAAAAACGGCCGAGGATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCAGAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAc                    <  1:80365/150‑1 (MQ=255)
                                                                                                                             tCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGTCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAc                    <  1:83390/150‑1 (MQ=255)
                                                                                                                              cAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAcc                   >  1:1099758/1‑150 (MQ=255)
                                                                                                                               aGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg                  <  2:1161203/150‑1 (MQ=255)
                                                                                                                               aGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccg                  >  1:328563/1‑150 (MQ=255)
                                                                                                                                gAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACAGCAGCGGCATTGCCACACGCACATTGCTAATCCGCGGTTTTGGCCATACACCCCGCACCAGTACCACCAGTGCAGAGATCAGAACAACGATTccgc                 >  1:865171/1‑150 (MQ=255)
                                                                                                                                  aaaaCCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccgccc               >  1:863141/1‑150 (MQ=255)
                                                                                                                                   aaaCCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATAccgcccg              >  1:896718/1‑150 (MQ=255)
                                                                                                                                        gccgAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACCACGATACCGCCCGCGCTg         >  2:747932/1‑150 (MQ=255)
                                                                                                                                        gccgAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTg         <  1:1132746/150‑1 (MQ=255)
                                                                                                                                              cATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTACTGc   >  1:818816/1‑150 (MQ=255)
                                                                                                                                               aTG‑TCGGCGG‑CTTTCTGACTTGCAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTACTCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTACTGcc  <  1:125086/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GGAACGCTATGATTTTAATATCGATGAATGGCGGGCAAAATTCAGGCAGTTATGTTTAGAGTATTTTGGCGATAGTCGTACACAGTTCACTTAACGCTCTGCGTCTATCTACTATCCCTTCTTTTTCAGAAAAACCGCCGAGCATG‑TCGGCGG‑CTTTCTGACTTACAAAGTAACAGATTACATCAGCGGCATTGCCAAACGCACAATGCTAATCAGCGGTTGTGGCCATACACCCAGCACCAGTACCAACAGTGCAGAGATCAGCACAACGATACCGCCCGCGCTGTACTGCC  >  NC_000913/2389871‑2390163

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