Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,063,343 T→C A148A (GCA→GCG cobU ← cobinamide kinase/cobinamide phosphate guanylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,063,3430TC99.4% 493.7 / ‑4.9 155A148A (GCA→GCGcobUcobinamide kinase/cobinamide phosphate guanylyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (84/70);  total (85/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 161 total aligned reads displayed.

TACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAA  >  NC_000913/2063253‑2063432
                                                                                          |                                                                                         
tACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCgcgc                                                                                         >  1:4010272/1‑93 (MQ=255)
  cgccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       <  1:386042/91‑1 (MQ=255)
  cTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       >  1:1039963/1‑93 (MQ=255)
  cTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCa                                                                                       <  1:3983862/93‑1 (MQ=255)
   tCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAg                                                                                      >  2:192682/1‑93 (MQ=255)
   tCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAg                                                                                      >  1:245181/1‑93 (MQ=255)
    ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAga                                                                                     >  2:812638/1‑93 (MQ=255)
    ccAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAga                                                                                     >  2:2404602/1‑93 (MQ=255)
     cAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgac                                                                                    <  2:3799308/93‑1 (MQ=255)
      aaTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCAGGAAAATGTCGCGCCAgacg                                                                                   >  2:1169464/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  2:2828106/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  1:1620923/1‑93 (MQ=255)
       aTACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAgacga                                                                                  >  2:2135736/1‑93 (MQ=255)
         aCCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGAct                                                                                >  1:3445220/1‑93 (MQ=255)
          cccGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGActc                                                                               >  1:1681394/1‑93 (MQ=255)
          cccGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCACACGActc                                                                               >  1:2565478/1‑93 (MQ=255)
             gAACCCCCCAGCCATACTTCATTTGCCGCAGCGGCCGACTGCTGATTTGCCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            <  2:1260843/93‑1 (MQ=255)
             gAAACCCCCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            <  1:3876332/93‑1 (MQ=255)
             gAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACtctca                                                                            >  2:1788490/1‑92 (MQ=255)
             gAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTcc                                                                            >  2:3040528/1‑93 (MQ=255)
                accaccAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGc                                                                         >  2:4630950/1‑93 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  1:3100301/92‑1 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  2:1229704/92‑1 (MQ=255)
                   cccaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      <  2:1229702/92‑1 (MQ=255)
                   accaGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGcaca                                                                      >  2:2471473/1‑93 (MQ=255)
                     caGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAAt                                                                    <  2:4453986/93‑1 (MQ=255)
                      aGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATc                                                                   <  2:2044425/93‑1 (MQ=255)
                      aGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATc                                                                   >  2:1779497/1‑93 (MQ=255)
                       gCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATcc                                                                  >  1:2610679/1‑93 (MQ=255)
                        ccATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATccc                                                                 >  1:2206810/1‑93 (MQ=255)
                        ccATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATccc                                                                 <  1:4576138/93‑1 (MQ=255)
                           tACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCAt                                                              <  2:3066031/93‑1 (MQ=255)
                            aCTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCAtc                                                             >  2:4019906/1‑92 (MQ=255)
                                  tttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACt                                                       >  1:4409739/1‑93 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      <  1:2730754/93‑1 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      <  2:1500536/93‑1 (MQ=255)
                                   ttGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACtt                                                      >  2:3602751/1‑93 (MQ=255)
                                           gCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTc                                              <  1:3533109/93‑1 (MQ=255)
                                             ggCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCAc                                            <  1:3010892/93‑1 (MQ=255)
                                              gCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACt                                           >  1:1775847/1‑93 (MQ=255)
                                                cacctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTaa                                         <  2:4213634/90‑1 (MQ=255)
                                                cAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTaa                                         <  1:2159262/93‑1 (MQ=255)
                                                 aaCTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAAt                                        >  2:16202/1‑93 (MQ=255)
                                                   ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAc                                      >  1:2125759/1‑93 (MQ=255)
                                                   ctgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATAc                                      <  2:4643175/93‑1 (MQ=255)
                                                    tgctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACa                                     >  1:763429/1‑93 (MQ=255)
                                                     gctgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACaa                                    >  1:4704950/1‑93 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCCCGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  <  2:2583710/93‑1 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  <  1:4567926/93‑1 (MQ=255)
                                                       tgATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAAcc                                  >  2:3012717/1‑93 (MQ=255)
                                                        gATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCt                                 >  2:3524311/1‑93 (MQ=255)
                                                        gATTTACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCt                                 >  2:2986742/1‑93 (MQ=255)
                                                         aTTTCCCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCtt                                <  2:3180356/93‑1 (MQ=255)
                                                          tttACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCttt                               >  1:2764866/1‑93 (MQ=255)
                                                           ttACCCGCCCGGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTg                              <  2:596030/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:1900202/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  2:1096927/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  2:1480110/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:3400807/93‑1 (MQ=255)
                                                              cccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgg                           <  1:3164283/93‑1 (MQ=255)
                                                               ccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                          <  2:1972551/93‑1 (MQ=255)
                                                               ccgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCggg                          <  2:4738482/93‑1 (MQ=255)
                                                                cgcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCgggg                         <  2:1531628/93‑1 (MQ=255)
                                                                 gcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTCGTCACTAATACCACCTTTGCGGGGc                        >  1:2302329/1‑93 (MQ=255)
                                                                 gcccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGc                        <  2:1282051/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      >  2:3383343/1‑93 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      <  1:280502/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      <  2:3771810/93‑1 (MQ=255)
                                                                   ccgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCaa                      >  2:1625370/1‑93 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     >  2:2135591/1‑93 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  2:248038/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  1:711122/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     <  2:3441621/93‑1 (MQ=255)
                                                                    cgGCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAAc                     >  2:1288349/1‑93 (MQ=255)
                                                                     ggCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACg                    >  2:4527374/1‑93 (MQ=255)
                                                                      gCAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGt                   <  2:420924/93‑1 (MQ=255)
                                                                       cAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACttt                  <  2:1503666/93‑1 (MQ=255)
                                                                       cAATATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGtt                  <  2:563466/93‑1 (MQ=255)
                                                                        aaTATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                 >  2:2303049/1‑93 (MQ=255)
                                                                        aaTATCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTg                 >  2:1382297/1‑93 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGCCGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                <  2:899754/93‑1 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                <  2:2316081/93‑1 (MQ=255)
                                                                         atatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTgg                >  1:376229/1‑93 (MQ=255)
                                                                          tatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc               >  2:1707089/1‑93 (MQ=255)
                                                                          aatCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggc               <  1:4079551/92‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              >  2:2005214/1‑93 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              <  1:4338537/93‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTggca              <  1:2469048/93‑1 (MQ=255)
                                                                           atCACGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACATTggca              >  2:1765626/1‑93 (MQ=255)
                                                                               cGAAAATGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTAACTAATACAACCTTTGCGGGGCAACGTTggcaggca          >  1:1872990/1‑93 (MQ=255)
                                                                                 aaaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagc        >  2:3337902/1‑93 (MQ=255)
                                                                                 aaaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagc        >  1:3534923/1‑93 (MQ=255)
                                                                                  aaaTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGgcagca       >  1:505281/1‑93 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCCCTAATACAACCTTTGCGGTGCAACGTTGGCAGGCAGCAAt     <  2:770811/93‑1 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAt     >  1:3119109/1‑93 (MQ=255)
                                                                                    aTGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCATTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAAt     >  1:1098692/1‑93 (MQ=255)
                                                                                     tGTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATc    >  2:4593668/1‑93 (MQ=255)
                                                                                      gTCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCa   <  2:1999887/93‑1 (MQ=255)
                                                                                       tCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCaa  <  1:2814400/93‑1 (MQ=255)
                                                                                       tCGCGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCaa  <  2:2596175/93‑1 (MQ=255)
                                                                                          |                                                                                         
TACTCCAATACCCGAAACCACCAGCCATACTTCATTTGCCGCAGCGGCCAACTGCTGATTTACCCGCCCGGCAATATCACGAAAATGTCGTGCCAGACGACTCTCCGGCACAATCCCCATTCCCACTTCGTTAGTCACTAATACAACCTTTGCGGGGCAACGTTGGCAGGCAGCAATCAA  >  NC_000913/2063253‑2063432

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