Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,870,876 Δ1 bp coding (365/1764 nt) dgoA ← 2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,870,8750G.99.2% 1206.9 / ‑0.3 262coding (366/1764 nt)dgoA2‑oxo‑3‑deoxygalactonate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (164/96);  total (166/96)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 288 total aligned reads displayed.

CAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCA  >  NC_000913/3870783‑3870963
                                                                                            |                                                                                        
cAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcg                                                                                          >  2:876818/1‑93 (MQ=255)
cAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcg                                                                                          >  1:819642/1‑93 (MQ=255)
cAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcg                                                                                          <  1:3927047/93‑1 (MQ=255)
 aaaGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcgc                                                                                         >  1:3690983/1‑92 (MQ=255)
  aaGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcgcg                                                                                        >  2:1961090/1‑91 (MQ=255)
    gACTGCGATGTCCGATGGCAATACCGCTTTTAGCGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcgcgcc                                                                                      <  2:500253/93‑5 (MQ=255)
    gACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTgcgcgcc                                                                                      >  1:728384/1‑89 (MQ=255)
     aCTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGccc                                                                                    >  1:4870284/1‑93 (MQ=255)
      cTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCg                                                                                   >  1:5047129/1‑93 (MQ=255)
       tGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGc                                                                                  <  2:137028/93‑1 (MQ=255)
       tGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGc                                                                                  >  1:3376733/1‑93 (MQ=255)
        gCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCt                                                                                 >  2:72143/1‑93 (MQ=255)
         cGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCtt                                                                                <  2:5629546/93‑1 (MQ=255)
          gATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTc                                                                               >  2:1652014/1‑93 (MQ=255)
           aTGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCg                                                                              <  2:4926823/93‑1 (MQ=255)
            tGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGa                                                                             <  2:1560651/93‑1 (MQ=255)
             gTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAg                                                                            <  2:1253635/93‑1 (MQ=255)
              tCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAgc                                                                           >  2:5448514/1‑93 (MQ=255)
              tCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAgc                                                                           >  2:478851/1‑93 (MQ=255)
                  aTGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCgg                                                                       >  2:5534208/1‑93 (MQ=255)
                   tGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGt                                                                      >  2:4731382/1‑93 (MQ=255)
                      cAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTaaa                                                                   <  1:603585/93‑1 (MQ=255)
                       aaTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGACCGGTAAAg                                                                  >  2:1389753/1‑93 (MQ=255)
                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAgg                                                                 >  2:4867805/1‑93 (MQ=255)
                        aTACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGGGCGCCCGCTTCGAGCGCGGTAAAgg                                                                 >  1:3830279/1‑93 (MQ=255)
                        aTACCGCTTTTAACGCTTTGATGCATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAgg                                                                 >  1:1285089/1‑93 (MQ=255)
                         tACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGc                                                                <  2:2558173/93‑1 (MQ=255)
                           ccGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCtt                                                              >  1:5161038/1‑93 (MQ=255)
                            cGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTc                                                             <  2:1148893/93‑1 (MQ=255)
                            cGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTc                                                             >  1:2611581/1‑93 (MQ=255)
                             gCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCg                                                            >  1:193677/1‑93 (MQ=255)
                              cTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGGCCGCTTCGAGCGCGGTAAAGGCTTCgg                                                           >  1:3167868/1‑93 (MQ=255)
                              cTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCgg                                                           >  1:3771889/1‑93 (MQ=255)
                              cTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCgg                                                           >  1:3618350/1‑93 (MQ=255)
                                tttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTc                                                         <  2:2707555/93‑1 (MQ=255)
                                tttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTc                                                         <  2:3979055/93‑1 (MQ=255)
                                 ttAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCg                                                        <  1:2498520/93‑1 (MQ=255)
                                  tAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGc                                                       <  2:3574762/93‑1 (MQ=255)
                                   aaCGCTTTGATGTATTGCGGACCAAAAGCCGATGAGGGAAATATTTTCAGCTCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGcc                                                      >  1:4009211/1‑93 (MQ=255)
                                    aCGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCg                                                     >  1:1780694/1‑93 (MQ=255)
                                     cGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGt                                                    <  2:2018017/93‑1 (MQ=255)
                                     cGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGt                                                    <  1:5681825/93‑1 (MQ=255)
                                      gCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTc                                                   <  2:5191134/93‑1 (MQ=255)
                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  1:879863/1‑93 (MQ=255)
                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  1:2197627/1‑93 (MQ=255)
                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  2:6125822/1‑93 (MQ=255)
                                        tttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgc                                                 >  2:785222/1‑93 (MQ=255)
                                         ttGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTcgcg                                                >  2:3624368/1‑93 (MQ=255)
                                           gATGTATTGCGGACAAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCa                                              <  1:1220349/93‑1 (MQ=255)
                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            >  1:1714518/1‑93 (MQ=255)
                                             tGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGc                                            >  2:3470436/1‑93 (MQ=255)
                                              gTATTGCTGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGcc                                           <  1:5327298/93‑1 (MQ=255)
                                              gTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGcc                                           >  1:4708055/1‑93 (MQ=255)
                                               tATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCg                                          <  1:5404696/93‑1 (MQ=255)
                                                aTTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCgg                                         <  2:2332636/93‑1 (MQ=255)
                                                   gCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGc                                      <  2:1444156/93‑1 (MQ=255)
                                                     ggACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAg                                    >  1:2035776/1‑93 (MQ=255)
                                                     ggACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAg                                    >  1:2745471/1‑93 (MQ=255)
                                                      gACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAgc                                   >  2:4481174/1‑92 (MQ=255)
                                                         cAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACgg                                >  1:958440/1‑93 (MQ=255)
                                                         cAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACgg                                >  1:3157520/1‑93 (MQ=255)
                                                         cAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACgg                                >  2:4218765/1‑93 (MQ=255)
                                                          aaaaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGAAGAgggc                               >  1:4991232/1‑92 (MQ=255)
                                                           aaaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTc                              >  1:107289/1‑93 (MQ=255)
                                                           aaaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTc                              <  1:4242520/93‑1 (MQ=255)
                                                           aaaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTc                              >  2:4375452/1‑93 (MQ=255)
                                                             aGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCAt                            <  1:6074295/93‑1 (MQ=255)
                                                             aGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCAt                            >  1:2474024/1‑93 (MQ=255)
                                                              gCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATg                           >  1:1897367/1‑93 (MQ=255)
                                                               gcGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGc                          <  2:3111915/92‑1 (MQ=255)
                                                               ccGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGTGCAGCCGGGGCAGACGGTCATGc                          >  1:4272128/1‑93 (MQ=255)
                                                               ccGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGc                          <  2:2883480/93‑1 (MQ=255)
                                                               ccGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGc                          >  1:2063268/1‑93 (MQ=255)
                                                                cGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGcc                         <  2:3218041/93‑1 (MQ=255)
                                                                   tGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      >  1:3087116/1‑93 (MQ=255)
                                                                   tGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      <  2:3767567/93‑1 (MQ=255)
                                                                   tGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTa                      <  2:106937/93‑1 (MQ=255)
                                                                    gACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTGCTGCCGTaa                     >  2:4538272/1‑92 (MQ=255)
                                                                     aCGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGc                    >  1:372855/1‑93 (MQ=255)
                                                                       ggAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGccc                  >  1:447856/1‑93 (MQ=255)
                                                                       ggAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGccc                  >  2:2079151/1‑93 (MQ=255)
                                                                        gAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCa                 >  1:3258317/1‑93 (MQ=255)
                                                                        gAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCa                 <  2:370852/93‑1 (MQ=255)
                                                                        gAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCa                 <  2:5750208/93‑1 (MQ=255)
                                                                        gAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCa                 >  1:1432197/1‑93 (MQ=255)
                                                                            tatTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGc             >  2:6064154/1‑93 (MQ=255)
                                                                              ttttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAc           >  2:470165/1‑93 (MQ=255)
                                                                               tttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcg          >  1:5605269/1‑93 (MQ=255)
                                                                               tttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcg          <  2:4049879/93‑1 (MQ=255)
                                                                               tttCAGCGCCTGCGCCCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcg          <  2:2648803/93‑1 (MQ=255)
                                                                                ttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcgg         >  2:2470702/1‑93 (MQ=255)
                                                                                ttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcgg         <  1:3020979/93‑1 (MQ=255)
                                                                                 tCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcggc        <  2:2139472/93‑1 (MQ=255)
                                                                                   aGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcggcgg      <  2:1775259/93‑1 (MQ=255)
                                                                                   aGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCAcggcgg      <  1:4365320/93‑1 (MQ=255)
                                                                                    gcgcCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGa     <  2:3352178/93‑1 (MQ=255)
                                                                                     cgcCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAt    <  2:4276125/93‑1 (MQ=255)
                                                                                      gcCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtc   >  2:2782945/1‑93 (MQ=255)
                                                                                       ccTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtca  <  2:1490852/93‑1 (MQ=255)
                                                                                       ccTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGAtca  <  2:1723640/93‑1 (MQ=255)
                                                                                            |                                                                                        
CAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCA  >  NC_000913/3870783‑3870963

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