Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,840,547 T→C 100% intergenic (+39/+2) yicS → / ← nepI putative periplasmic protein/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,840,5470TC99.0% 325.9 / ‑4.7 97intergenic (+39/+2)yicS/nepIputative periplasmic protein/putative transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (50/45);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATAT  >  NC_000913/3840451‑3840643
                                                                                                |                                                                                                
taataaTCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATc                                                                                               >  2:2527775/1‑100 (MQ=255)
  ataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAg                                                                                             <  2:1439408/100‑1 (MQ=255)
  ataatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAg                                                                                             >  1:1562530/1‑100 (MQ=255)
   taatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAgg                                                                                            <  2:2417531/100‑1 (MQ=255)
    aatCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGa                                                                                           >  1:120341/1‑100 (MQ=255)
      tCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGAtt                                                                                         <  1:1352162/100‑1 (MQ=255)
       cAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGCACAGCAACAACTACTGACCAACGCTTATTTACCTGACATCAGGAtt                                                                                         <  2:711751/99‑1 (MQ=255)
        agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                       <  1:2288820/100‑1 (MQ=255)
        agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                       <  2:1922114/100‑1 (MQ=255)
        agagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttc                                                                                       >  1:2139559/1‑100 (MQ=255)
         gagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttct                                                                                      <  2:692846/100‑1 (MQ=255)
         gagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttct                                                                                      <  1:2527790/100‑1 (MQ=255)
         gagCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACAAACGCTTCTTTACCTGACATCAGGATttct                                                                                      >  2:2255472/1‑100 (MQ=255)
         gagCAAGTATTGGGAAGCGGTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAAAGCTTCTTTACCTGACATCAGGATttc                                                                                       >  2:1206657/1‑99 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                    <  1:2414776/100‑1 (MQ=255)
           gCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATttcttc                                                                                    <  2:775232/100‑1 (MQ=255)
               ttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCAttt                                                                                >  2:614011/1‑100 (MQ=255)
               ttattGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGAATTTCTTCAttt                                                                                >  2:1493037/1‑100 (MQ=255)
                  ttGAGAAGCGTTGGGTAAGATAGAATGCCCGGATATGTAAGTACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTc                                                                              <  2:1796862/99‑1 (MQ=255)
                   tGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAc                                                                            >  1:1295609/1‑100 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                           <  1:1979833/100‑1 (MQ=255)
                    gggAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCAcc                                                                           >  2:2104711/1‑100 (MQ=255)
                      gAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCtt                                                                         <  1:456233/100‑1 (MQ=255)
                         gCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  1:2141999/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTTTTCATTTTCACCTTTGCAg                                                                    >  2:2068485/1‑100 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGc                                                                      <  2:1597602/98‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                    <  1:1918352/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                    <  1:2513360/100‑1 (MQ=255)
                           gTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAg                                                                    <  1:1428432/100‑1 (MQ=255)
                            ttGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGt                                                                   <  1:1450197/100‑1 (MQ=255)
                               ggTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                >  2:253914/1‑100 (MQ=255)
                                gTACGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAc                                                                <  1:144299/99‑1 (MQ=255)
                                gTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAAcc                                                               <  1:2521420/100‑1 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                            <  2:2261514/100‑1 (MQ=255)
                                   aGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACATGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaac                                                            >  2:1174572/1‑100 (MQ=255)
                                    ggTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAgatcggaaga                                                           >  1:2063971/1‑91 (MQ=255)
                                     gTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                          >  2:1610178/1‑100 (MQ=255)
                                      tAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACcaacaa                                                          <  2:2481946/99‑1 (MQ=255)
                                       aGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATg                                                        <  2:1952929/100‑1 (MQ=255)
                                       aGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTATTCATTTTCACCTTTGCAGTAACCAACAATg                                                        >  1:2032860/1‑100 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                       <  1:2108849/100‑1 (MQ=255)
                                        gAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                       <  1:1656728/100‑1 (MQ=255)
                                         aaTGCCCGGATATGTAAGAACAGCAACATCTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGc                                                       <  2:89847/99‑1 (MQ=255)
                                           tGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGt                                                    >  2:2457125/1‑100 (MQ=255)
                                            gCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTc                                                   >  2:2537917/1‑100 (MQ=255)
                                                ggATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCa                                               <  1:22602/100‑1 (MQ=255)
                                                   tatGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACa                                            <  1:1742337/100‑1 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                         >  2:2318537/1‑100 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCa                                         >  2:2207578/1‑100 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTCCAGTAACCAACAATGCAGTCAGCAACATc                                          >  2:356300/1‑99 (MQ=255)
                                                      gTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTACATTTTCACCTTTGCACTAACCAACAATGCAATCAGCAACCTCa                                         >  2:203359/1‑100 (MQ=255)
                                                       tAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAAAATCaa                                        >  2:1376645/1‑100 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                      <  1:1346972/100‑1 (MQ=255)
                                                         aGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                      >  2:2400937/1‑100 (MQ=255)
                                                           aaCAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAAtg                                      <  2:1719546/98‑1 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                  <  1:218697/100‑1 (MQ=255)
                                                             cAGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGcc                                  >  2:120129/1‑100 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTTTTTTTTTTTTTAACCAACAATGCAGTCAACAACATccttatgcc                                  >  2:551261/1‑91 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCg                                 >  2:289505/1‑100 (MQ=255)
                                                              aGCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCg                                 >  1:1067348/1‑100 (MQ=255)
                                                               gCAACAACTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCgg                                >  1:592906/1‑100 (MQ=255)
                                                                caacaaCTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGa                               <  1:2437472/100‑1 (MQ=255)
                                                                 aacaacTACTGACCAATGTTTCTTTAACTTCAAATAGGATATCTTCTTATTAACCTTTGTAGACACCAACAATATAATTCGCATCATTATTGTGACGGac                              >  2:992367/1‑100 (MQ=255)
                                                                 aacaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGac                              <  2:2530024/100‑1 (MQ=255)
                                                                  acaacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                             <  1:1930265/100‑1 (MQ=255)
                                                                   caacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGacaa                            >  1:34605/1‑100 (MQ=255)
                                                                   caacTACTGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGaca                             >  2:1135952/1‑99 (MQ=255)
                                                                      ctactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACAt                         <  1:1135851/100‑1 (MQ=255)
                                                                      cgactGACCAACGCTTCTTTACCTTACATCAGGATTTCTTCATTTTCACCTTTGCAGTAATCAACAATGCAGTCAGCAACATCAATGTGCCGGacaaca                          <  2:1221275/97‑1 (MQ=255)
                                                                       tactGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATc                        >  1:148767/1‑100 (MQ=255)
                                                                        actGACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                       <  2:145993/100‑1 (MQ=255)
                                                                        actGACCAACGCTTCTTTACCTGACATCAGAATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCa                       <  2:1782082/100‑1 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCaa                      >  2:1778680/1‑98 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                    >  1:1293949/1‑100 (MQ=255)
                                                                           gACCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                    >  2:1488995/1‑100 (MQ=255)
                                                                            aCCAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACgg                   >  2:1708492/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACg                    >  2:664827/1‑98 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  2:2254348/1‑100 (MQ=255)
                                                                             ccAACGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGc                  >  1:1424768/1‑100 (MQ=255)
                                                                             ccAAAGCATCTTTACCAGACATCAGGATTTCTTCATTATCAACTTTGCAGTAACCAACAATCCAGTCAGCAACAGTATTGTGGCTGACAACATCAActtc                  >  2:473944/1‑97 (MQ=255)
                                                                               aaCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGa                >  2:1952086/1‑100 (MQ=255)
                                                                                aCGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGaa               <  1:352112/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:1316075/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:2104725/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              >  2:23585/1‑100 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:1610189/100‑1 (MQ=255)
                                                                                 cGCTTCTTTACCTGACATCAGGATTTCTTCATTCTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAg              <  1:1248738/100‑1 (MQ=255)
                                                                                     tcttTACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACAGCAACGGCGAAGTCAg          >  1:2515563/1‑100 (MQ=255)
                                                                                       tttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAc        >  2:1322842/1‑100 (MQ=255)
                                                                                        ttACCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGAcc       <  2:419079/100‑1 (MQ=255)
                                                                                          aCCTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCaa     >  1:2363511/1‑100 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat    >  1:2018777/1‑100 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat    >  2:1153625/1‑100 (MQ=255)
                                                                                           ccTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAat    >  2:1777179/1‑100 (MQ=255)
                                                                                            cTGACATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAata   >  1:2231513/1‑100 (MQ=255)
                                                                                             tGACATCATGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAatat  <  1:2394078/100‑1 (MQ=255)
                                                                                             tGACATCAGGATTGCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAatat  >  1:1742679/1‑100 (MQ=255)
                                                                                                |                                                                                                
TAATAATCAGAGCAATTATTGGGAAGCGTTGGGTAAGGTAGAATGCCCGGATATGTAAGAACAGCAACAACTACTGACCAACGCTTCTTTACCTGATATCAGGATTTCTTCATTTTCACCTTTGCAGTAACCAACAATGCAGTCAGCAACATCAATGTGCCGGACAACATCAACGGCGAAGTCAGACCAATAT  >  NC_000913/3840451‑3840643

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: