Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,065,747 (C)6→5 100% coding (18/546 nt) cobU ← cobinamide kinase and cobinamide phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,065,7420C.86.3% 263.0 / 41.6 80coding (23/546 nt)cobUcobinamide kinase and cobinamide phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base . (38/31);  total (45/35)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

GCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAA  >  NC_000913/2065643‑2065840
                                                                                                   |                                                                                                  
gCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTgcgc                                                                                                    <  1:913443/100‑1 (MQ=255)
gCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTgcgc                                                                                                    >  2:1531465/1‑100 (MQ=255)
  cATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGccc                                                                                                  >  2:660875/1‑100 (MQ=255)
    tctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGccccc                                                                                                >  1:2568147/1‑100 (MQ=255)
    tctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGccccc                                                                                                >  2:1373295/1‑100 (MQ=255)
     ctcATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccccg                                                                                               <  1:713513/100‑2 (MQ=255)
       catcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccccgtc                                                                                             <  1:2433095/100‑4 (MQ=255)
       catcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccccgtc                                                                                             <  1:2528418/100‑4 (MQ=255)
       catcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGcccccgtc                                                                                             <  2:296385/100‑4 (MQ=255)
       catcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGt                                                                                             >  1:471952/1‑100 (MQ=255)
         tcatcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCAc                                                                                          >  1:1451699/1‑100 (MQ=255)
          catcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCAcc                                                                                         >  2:894775/1‑100 (MQ=255)
           atcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCa                                                                                        >  1:1134254/1‑100 (MQ=255)
            tcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCa                                                                                        >  2:2138134/1‑100 (MQ=255)
            tcaAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCaa                                                                                       <  1:123896/100‑1 (MQ=255)
               aGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTGCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAat                                                                                    >  2:1809015/1‑100 (MQ=255)
                   tttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAatcatc                                                                                >  2:1678858/1‑100 (MQ=255)
                   tttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAatcatc                                                                                <  2:2158487/100‑1 (MQ=255)
                    ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAAtcatca                                                                               <  2:2512700/100‑1 (MQ=255)
                    ttGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCATCCGCATGGCGACCCCTCCCGCACCGTGCGCCCCCGTCACCAAAAtcatca                                                                               >  1:94100/1‑100 (MQ=255)
                     tGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATcatcat                                                                              >  2:1087737/1‑100 (MQ=255)
                     tGCGAGGTAGCGATATACAGAAACTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATcatcat                                                                              >  1:1503720/1‑100 (MQ=255)
                      gCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATa                                                                             <  2:2155481/100‑1 (MQ=255)
                        gAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATAc                                                                            >  1:1659161/1‑100 (MQ=255)
                        gAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATAcc                                                                           >  2:1560881/1‑100 (MQ=255)
                            tAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACt                                                                       <  2:527663/100‑1 (MQ=255)
                              gCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTcc                                                                     <  1:1775981/100‑1 (MQ=255)
                              gCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTca                                                                     <  1:335180/100‑2 (MQ=255)
                               cGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCt                                                                    >  1:1392740/1‑100 (MQ=255)
                                  tataCAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGAt                                                                 <  2:1393126/100‑1 (MQ=255)
                                    taCAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                               >  1:577538/1‑100 (MQ=255)
                                     aCAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTc                                                               <  1:2291452/100‑1 (MQ=255)
                                      cAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCaa                                                             >  1:472047/1‑100 (MQ=255)
                                       aGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCaaa                                                            <  2:1525319/100‑1 (MQ=255)
                                       aGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCaaa                                                            <  2:1424816/100‑1 (MQ=255)
                                        gAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAAt                                                           >  1:1454059/1‑100 (MQ=255)
                                        gAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAAt                                                           >  1:165887/1‑100 (MQ=255)
                                        gAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAAt                                                           <  2:1974420/100‑1 (MQ=255)
                                         aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATg                                                          <  2:2156342/100‑1 (MQ=255)
                                         aaCCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATg                                                          <  2:1552565/100‑1 (MQ=255)
                                           ccTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                        >  1:183954/1‑100 (MQ=255)
                                           ccTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTa                                                        <  1:168629/100‑1 (MQ=255)
                                             tgtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatg                                                      <  2:2529857/100‑1 (MQ=255)
                                              gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCCTACCGACTCCTGATTCAAATGt                                                         >  1:2265530/1‑96 (MQ=255)
                                              gtgAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatga                                                     >  2:1560029/1‑100 (MQ=255)
                                                gAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatgat                                                    >  2:1764324/1‑100 (MQ=255)
                                                gAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTatgatg                                                   >  2:1184564/1‑100 (MQ=255)
                                                 aaGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCTCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATgg                                                  <  1:1432469/100‑1 (MQ=255)
                                                 aaGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATgg                                                  <  2:1202353/100‑1 (MQ=255)
                                                 aaGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATgg                                                  <  2:2159188/100‑1 (MQ=255)
                                                 aaGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATgg                                                  <  2:1823411/100‑1 (MQ=255)
                                                  agagTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAGATCATTATACCGACTCCTGATTCAAATGTATGATGGt                                                 >  1:1745889/1‑100 (MQ=255)
                                                     gTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTAt                                               >  2:1384674/1‑100 (MQ=255)
                                                     ctCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCATGCGCCCCCGTCACAAAAATCATCATACAGACCCCAGATTCCAATGTATGATGGTAtt                                              >  1:1238937/2‑100 (MQ=255)
                                                      tCCCCAATAAGCGCCTCTGCATGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTagt                                              >  1:1576160/1‑98 (MQ=255)
                                                        cccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGCTTCAAATGTATGATGGTATTCAt                                           >  1:174652/1‑100 (MQ=255)
                                                        cccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCAt                                           >  2:628880/1‑100 (MQ=255)
                                                         ccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGCCTCCTGATTCAAATGTATGATGGTATTCAtt                                          >  1:376732/1‑100 (MQ=255)
                                                         ccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCAtt                                          >  2:1116897/1‑100 (MQ=255)
                                                         ccAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCAtt                                          <  2:587264/100‑1 (MQ=255)
                                                          cAATAAGCGCCTCTTCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCAtt                                          <  1:1182892/100‑1 (MQ=255)
                                                          cAATAAGCGCCTCTGCGTGGCGACTCTTCCTGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTg                                         <  2:568424/100‑1 (MQ=255)
                                                           aaTAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTga                                        <  1:1583588/100‑1 (MQ=255)
                                                              aaGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaaca                                     <  1:745299/100‑1 (MQ=255)
                                                              aaGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATAATGGTATTCATTtaaca                                     >  1:349993/1‑100 (MQ=255)
                                                               aGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCGTCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacag                                    >  2:262803/1‑100 (MQ=255)
                                                               aGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacag                                    >  1:3993/1‑100 (MQ=255)
                                                               aGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacag                                    >  2:2361583/1‑100 (MQ=255)
                                                                gcgcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAACTGTATGATGGTATTCATTgaacag                                    >  2:1686567/1‑100 (MQ=255)
                                                                gcgcCTCTGCGTGGCGACTATTCCGGCTCCGTGCGCCCCCGTCACTAAAATCATCGTACCGTCGCCTGATTCAAATGTATGATGGTATTCATTgaac                                      <  1:234782/97‑1 (MQ=255)
                                                                  gcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaac                                 <  2:1568409/100‑1 (MQ=255)
                                                                  gcCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaac                                 <  1:660872/100‑1 (MQ=255)
                                                                   ccTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaac                                 <  1:1995068/100‑1 (MQ=255)
                                                                   ccTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTgaacagaaca                                >  2:246973/1‑100 (MQ=255)
                                                                      ctGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATg                             <  2:1503729/100‑1 (MQ=255)
                                                                        gCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATgg                            <  1:1686558/100‑1 (MQ=255)
                                                                        gCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGa                           >  2:1053776/1‑100 (MQ=255)
                                                                             gCGACTCTTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATcccc                      >  2:843405/1‑100 (MQ=255)
                                                                                 ctctTCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCAcc                  <  2:2157702/100‑1 (MQ=255)
                                                                                    ttCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGtt               <  1:383819/100‑1 (MQ=255)
                                                                                    ttCCCGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGtt               >  2:2209644/1‑100 (MQ=255)
                                                                                      cccGCTCCGTGCGCCCCCGTCACCAAAATCATCATAGCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc             <  2:1817742/100‑1 (MQ=255)
                                                                                      cccGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc             >  1:1327868/1‑100 (MQ=255)
                                                                                      cccGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGc             >  1:281069/1‑100 (MQ=255)
                                                                                       ccGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCt            >  2:2263348/1‑100 (MQ=255)
                                                                                        cGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTg           >  1:1660127/1‑100 (MQ=255)
                                                                                        cGCTCCGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTg           <  1:2263333/100‑1 (MQ=255)
                                                                                             cGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGc      >  1:2080317/1‑100 (MQ=255)
                                                                                             cGTGCGCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGc      >  2:1684570/1‑100 (MQ=255)
                                                                                                  gcCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCaaaa  <  2:1908789/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
GCCATCTCATCATCAAGGATTTGCGAGGTAGCGATATACAGAACCTGTGAAGAGTCCCCAATAAGCGCCTCTGCGTGGCGACTCTTCCCGCTCCGTGCGCCCCCCGTCACCAAAATCATCATACCGACTCCTGATTCAAATGTATGATGGTATTCATTGAACAGAACAATGGATCCCCCACCGTTGCTGACCGCAAAA  >  NC_000913/2065643‑2065840

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: