Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,806,042 (C)6→5 100% coding (42/759 nt) ydiY ← acid‑inducible putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,806,0370C.77.1% 195.4 / 79.3 83coding (47/759 nt)ydiYacid‑inducible putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base . (35/29);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.08e-01

CGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACT  >  NC_000913/1805939‑1806134
                                                                                                  |                                                                                                 
cGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcc                                                                                                  >  1:1338494/1‑100 (MQ=255)
 gTTCAGATTACCTTCGAAAGGTTTTTTTTTTTTTTTAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGTCGCAAACATGccc                                                                                                 >  1:1883955/1‑100 (MQ=255)
 gTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAATAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGccc                                                                                                 >  1:1869794/1‑100 (MQ=255)
  ttCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccc                                                                                                <  1:2478216/100‑1 (MQ=255)
      gATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGcccccgcc                                                                                            <  1:1781760/100‑4 (MQ=255)
        ttACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCaa                                                                                         <  2:804104/100‑1 (MQ=255)
        ttACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCaa                                                                                         <  2:243496/100‑1 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACAt                                                                                      >  2:538344/1‑100 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACAt                                                                                      >  2:1275929/1‑100 (MQ=255)
           ccTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACAt                                                                                      >  1:655084/1‑100 (MQ=255)
               cGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACg                                                                                  >  2:2092073/1‑100 (MQ=255)
                gAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGa                                                                                 >  2:1662493/1‑100 (MQ=255)
                gAAAGGCTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGa                                                                                 >  1:1129851/1‑100 (MQ=255)
                 aaaGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGAt                                                                                <  2:1207700/100‑1 (MQ=255)
                   aGGTTTTTTAGCGCTGGCAAGTTCATCCATGACAGTAAAACCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGAtn                                                                               <  2:1093710/99‑2 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGc                                                                             <  1:1857838/100‑1 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGc                                                                             >  2:628677/1‑100 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGc                                                                             >  1:1418590/1‑100 (MQ=255)
                    ggTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGc                                                                             >  1:2036999/1‑100 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCa                                                                            <  1:2225085/100‑1 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCa                                                                            >  1:2380251/1‑100 (MQ=255)
                     gTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCa                                                                            <  1:81558/100‑1 (MQ=255)
                      ttttttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAg                                                                           <  2:1324700/100‑1 (MQ=255)
                      gatcttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAg                                                                           <  1:211345/96‑1 (MQ=255)
                        ttttAGCGCTGCCAGGATCATCCGTTACAGTAAAAACGGAACCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCa                                                                            <  1:125667/97‑1 (MQ=255)
                          ttAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAAc                                                                       <  2:2052459/100‑1 (MQ=255)
                           tAGCGCTGGCAGGATCCACCATGACAGTCAAAACGGAATTAGCGACAGCCGCACTCAGTGACGCAAACATGCCCTCGCCAACATAACGAGTGCAg                                                                           >  1:2424647/1‑95 (MQ=255)
                           tAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAgagatc                                                                      >  2:211347/1‑95 (MQ=255)
                           tAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACt                                                                      <  1:1345460/100‑1 (MQ=255)
                             gcgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTg                                                                     <  2:1433002/100‑1 (MQ=255)
                              cgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTc                                                                   >  1:135412/1‑100 (MQ=255)
                              cgcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTc                                                                   <  1:1275924/100‑1 (MQ=255)
                               gcTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCt                                                                  <  2:1914597/100‑1 (MQ=255)
                                cTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCt                                                                  <  2:2582590/100‑1 (MQ=255)
                                cTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCtt                                                                 <  1:2411119/100‑1 (MQ=255)
                                 tGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTc                                                                <  2:301978/100‑1 (MQ=255)
                                  ggCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCCACATAACGATGGCAGGAACTGTCTTc                                                                >  1:2245505/1‑100 (MQ=255)
                                  ggCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTc                                                                <  2:2337089/100‑1 (MQ=255)
                                   gCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCa                                                               <  1:1242755/100‑1 (MQ=255)
                                     aGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGTCAACATAACGATGGCAGGAACTGTCTTCaaa                                                             <  1:555521/100‑1 (MQ=255)
                                     aGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCaaaa                                                            <  2:1686423/100‑1 (MQ=255)
                                       gATCATCCATTACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                          <  1:1631973/100‑1 (MQ=255)
                                        atcatcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGc                                                          >  1:283200/1‑100 (MQ=255)
                                          catcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTc                                                       <  2:1349098/100‑1 (MQ=255)
                                            tcCCTGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGTCCCCAGCCAACATAACGATGGCAGGAACAGTCGTCAAAAGCTTCa                                                      >  1:17707/1‑100 (MQ=255)
                                            tcCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAt                                                     <  1:672020/100‑1 (MQ=255)
                                             ccATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAt                                                     <  2:1238538/100‑1 (MQ=255)
                                              cATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAttt                                                   >  1:240234/1‑100 (MQ=255)
                                              cATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCAttt                                                   <  1:1870873/100‑1 (MQ=255)
                                                tGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTa                                                  >  2:959367/1‑100 (MQ=255)
                                                tGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTAt                                                 >  1:1889351/1‑100 (MQ=255)
                                                tGACAGTAAAAACGGAACCATTGGCAGCTGCATTCACTGACGCAAACAAGCCCACGCCAACATAACGACGGCAGGAACTGTCTTCAAAAGCTTCATTTa                                                  <  1:2374477/99‑1 (MQ=255)
                                                  aCAGTAAATACGGAATTATCGGCAGCAGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGGACTGTCTTCAAAAGCTTCATTTATCa                                               >  1:561202/1‑100 (MQ=255)
                                                   cAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCaa                                              >  1:399988/1‑100 (MQ=255)
                                                    aGTAAAAACGGAATCATCGGCAGCTGTATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAACAGCTTCATTTATCAag                                             <  1:389745/100‑1 (MQ=255)
                                                     gTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAag                                             >  1:2576317/1‑100 (MQ=255)
                                                     gTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAaga                                            >  1:37123/1‑100 (MQ=255)
                                                       aaaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAagag                                           >  2:844722/1‑99 (MQ=255)
                                                        aaaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTc                                         <  1:64180/100‑1 (MQ=255)
                                                         aaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTc                                         >  2:1621758/1‑100 (MQ=255)
                                                         aaaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTcc                                        >  1:1616429/1‑100 (MQ=255)
                                                          aaCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCg                                       >  2:674332/1‑100 (MQ=255)
                                                           aCGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGt                                      >  2:936631/1‑100 (MQ=255)
                                                              gAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                    <  2:825915/100‑1 (MQ=255)
                                                              gAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTaca                                   >  2:2245018/1‑100 (MQ=255)
                                                                 tcatcGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTacaaca                                >  1:22512/1‑100 (MQ=255)
                                                                  catcGGCAGCGGCAGTCAGTGACGCAAACATGACCCCGCCAACAAAACGAAGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTac                                    <  1:6341/95‑1 (MQ=255)
                                                                   atcGGCAGCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACaaa                              >  1:684816/1‑100 (MQ=255)
                                                                          gCTGCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAagag                       >  2:2389588/1‑100 (MQ=255)
                                                                             gCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGAcc                    >  2:1930879/1‑100 (MQ=255)
                                                                             gCATTCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGAcc                    <  2:2364186/100‑1 (MQ=255)
                                                                                 tCAGTGACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcg                >  2:150962/1‑100 (MQ=255)
                                                                                   aGTGACGCAAACATGCCCCCCGCCAACATAACGAGGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATcgc               <  2:2162889/100‑1 (MQ=255)
                                                                                    ttGAAGCAAACATGCCCCCGTCAACATATCGATGGCAGGAACTATCTTCAAAAGCTTCATTTATCAAGAGCCCGTACAACAAAAAAAGAGACCATCGCgg             <  1:2156823/99‑1 (MQ=255)
                                                                                      gACGCAAACATGCCCCCGCCAACCTAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGACCGTACAACAAAAAAAGAGACCATAGCGGTc           >  2:2058338/1‑100 (MQ=255)
                                                                                      gACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTc           >  1:2359333/1‑100 (MQ=255)
                                                                                      gACGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGACACCATCGCGGTc           >  1:958239/1‑100 (MQ=255)
                                                                                       aCGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          <  2:70968/100‑1 (MQ=255)
                                                                                        cGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          >  2:473183/1‑100 (MQ=255)
                                                                                        cGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTcc          >  2:1331197/1‑100 (MQ=255)
                                                                                        cGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:1128960/1‑100 (MQ=255)
                                                                                        cGCAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:630792/1‑100 (MQ=255)
                                                                                        cGCAAACATGCCCCCGCCAACATAACGATGGCAGGAAATGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTccc         >  1:220529/1‑100 (MQ=255)
                                                                                          cAAACATGCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCgg       <  1:94623/100‑1 (MQ=255)
                                                                                           aaaCATGCCCCCGCCAACATAACGATGGCAGGAAATGTCTTCAAAAGCTTGATTTATCAGGAGTCCGTACATCAAAAAAAGAGACCATGGCGGTCCCGGa      >  1:1155902/1‑100 (MQ=255)
                                                                                              cATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaaa    >  1:2099634/1‑100 (MQ=255)
                                                                                              cATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGaaa    >  2:1406243/1‑100 (MQ=255)
                                                                                                tGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACt  >  2:1610918/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
CGTTCAGATTACCTTCGAAAGGTTTTTTAGCGCTGGCAGGATCATCCATGACAGTAAAAACGGAATCATCGGCAGCTGCATTCAGTGACGCAAACATGCCCCCCGCCAACATAACGATGGCAGGAACTGTCTTCAAAAGCTTCATTTATCAAGAGTCCGTACAACAAAAAAAGAGACCATCGCGGTCCCGGAAACT  >  NC_000913/1805939‑1806134

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