Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 552,300 C→T A8A (GCG→GCA purE ← phosphoribosylaminoimidazole carboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913552,3000CT99.3% 984.6 / ‑4.0 279A8A (GCG→GCApurEphosphoribosylaminoimidazole carboxylase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (117/159);  total (118/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-01

Only 100 of 288 total aligned reads displayed.

GTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGA  >  NC_000913/552210‑552389
                                                                                          |                                                                                         
 tggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGcc                                                                                        <  1:760157/93‑1 (MQ=255)
 tggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGcc                                                                                        <  1:6072408/93‑1 (MQ=255)
  ggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCg                                                                                       <  1:7375720/93‑1 (MQ=255)
  ggtgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCg                                                                                       <  1:718594/93‑1 (MQ=255)
    tgCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGa                                                                                     <  1:843464/93‑1 (MQ=255)
       gggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGatta                                                                                  >  1:672106/1‑93 (MQ=255)
        ggACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGattat                                                                                 <  2:7632673/93‑1 (MQ=255)
          aCATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTg                                                                               <  1:1044576/93‑1 (MQ=255)
           cATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGc                                                                              >  1:7930206/1‑93 (MQ=255)
           cATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGc                                                                              >  1:4360134/1‑93 (MQ=255)
            aTTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCg                                                                             <  2:4642727/93‑1 (MQ=255)
             ttCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCgg                                                                            >  2:2495664/1‑93 (MQ=255)
             ttCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCgg                                                                            <  1:3016809/93‑1 (MQ=255)
              tCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCggg                                                                           >  2:32559/1‑93 (MQ=255)
               cAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGa                                                                          >  1:3735048/1‑93 (MQ=255)
                aGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGaa                                                                         <  2:2750752/93‑1 (MQ=255)
                 ggATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAg                                                                        >  2:3072128/1‑93 (MQ=255)
                  gATTTCGAAGATTTCGGCGGCGAACTGCATGGTAACCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGa                                                                       >  2:6131103/1‑93 (MQ=255)
                   aTTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGaa                                                                      >  1:262719/1‑92 (MQ=255)
                   aTTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGAc                                                                      >  2:7531195/1‑93 (MQ=255)
                    tttCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACa                                                                     <  1:4123891/93‑1 (MQ=255)
                     ttCGAAGATGTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACAt                                                                    <  1:4527260/93‑1 (MQ=255)
                      tCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATg                                                                   >  2:2942226/1‑93 (MQ=255)
                        gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGCCACGTGCCGGATTATTGCGGGAAGACATGCg                                                                 >  1:5892129/1‑93 (MQ=255)
                        gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCg                                                                 <  2:4114184/93‑1 (MQ=255)
                        gAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCg                                                                 >  1:100240/1‑93 (MQ=255)
                         aaGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGt                                                                <  2:4811864/93‑1 (MQ=255)
                          aGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTc                                                               <  2:110762/93‑1 (MQ=255)
                          aGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTc                                                               <  2:1204694/93‑1 (MQ=255)
                             tttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTa                                                            >  2:396220/1‑93 (MQ=255)
                             tttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTa                                                            >  1:3110245/1‑93 (MQ=255)
                              ttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaa                                                           <  1:6617690/93‑1 (MQ=255)
                                cggcggCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTaaaa                                                         >  1:3340296/1‑93 (MQ=255)
                                 ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAAc                                                        >  1:290773/1‑93 (MQ=255)
                                 ggcggcGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAAc                                                        <  2:3560270/93‑1 (MQ=255)
                                    ggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTcc                                                     >  2:4036817/1‑93 (MQ=255)
                                    ggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTcc                                                     >  2:8535443/1‑93 (MQ=255)
                                    ggcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTcc                                                     <  1:6149190/93‑1 (MQ=255)
                                     gcgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCt                                                    >  2:787013/1‑93 (MQ=255)
                                      cgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtg                                                   <  1:4072241/93‑1 (MQ=255)
                                       gAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCtgt                                                  >  2:7508663/1‑93 (MQ=255)
                                         aCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtgg                                                >  2:7148339/1‑93 (MQ=255)
                                          cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                               <  1:3628801/93‑1 (MQ=255)
                                          cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                               <  2:1882060/93‑1 (MQ=255)
                                          cTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggt                                               <  2:5728372/93‑1 (MQ=255)
                                           tGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTgtggtg                                              <  2:5328280/93‑1 (MQ=255)
                                            gCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGc                                             >  2:3786195/1‑93 (MQ=255)
                                            gCATGCTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGc                                             <  1:6653391/93‑1 (MQ=255)
                                              aTGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGcac                                           <  2:4535110/93‑1 (MQ=255)
                                                   aGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActc                                      >  2:1011311/1‑93 (MQ=255)
                                                   aGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActc                                      <  2:6291856/93‑1 (MQ=255)
                                                    gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct                                     <  2:6584192/93‑1 (MQ=255)
                                                    gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct                                     <  2:5474526/93‑1 (MQ=255)
                                                    gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct                                     >  1:1687976/1‑93 (MQ=255)
                                                    gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActct                                     <  2:6399836/93‑1 (MQ=255)
                                                    gCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAActcc                                     >  2:3811770/1‑92 (MQ=255)
                                                     cccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACtctc                                    <  1:8212157/93‑1 (MQ=255)
                                                     cccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACtctc                                    <  1:4192425/93‑1 (MQ=255)
                                                     cccACTCGCTTTTGGCCCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACtctc                                    <  2:7397623/93‑1 (MQ=255)
                                                      ccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCg                                   >  1:3440543/1‑93 (MQ=255)
                                                      ccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCg                                   >  1:3843448/1‑93 (MQ=255)
                                                      ccAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCg                                   <  2:8099780/93‑1 (MQ=255)
                                                       cAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                  >  2:1395983/1‑93 (MQ=255)
                                                       cAGTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCgg                                  >  1:4445546/1‑93 (MQ=255)
                                                         gTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTGTTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                                <  1:1406396/93‑1 (MQ=255)
                                                         gTCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCt                                >  2:3387951/1‑93 (MQ=255)
                                                          tCGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCtt                               <  1:6989791/93‑1 (MQ=255)
                                                           cGCTTTTGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCttt                              <  1:3916654/93‑1 (MQ=255)
                                                           cGCTTTTGGACCCCATCACGATGACGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCttt                              <  2:5993043/93‑1 (MQ=255)
                                                              ttttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTaga                           <  1:5064993/93‑1 (MQ=255)
                                                              ttttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTaga                           <  1:2185849/93‑1 (MQ=255)
                                                                ttGGACCCCATCGCGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                         <  2:5679669/93‑1 (MQ=255)
                                                                ttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                         <  1:5637123/93‑1 (MQ=255)
                                                                ttGGACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAgg                         <  1:594840/93‑1 (MQ=255)
                                                                  ggACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGc                       <  2:8559632/93‑1 (MQ=255)
                                                                   gACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGca                      <  2:2031582/93‑1 (MQ=255)
                                                                   gACCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGca                      <  1:8021090/93‑1 (MQ=255)
                                                                    aCCCCATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGccc                     >  1:1635692/1‑91 (MQ=255)
                                                                      cccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACag                   >  2:685661/1‑93 (MQ=255)
                                                                       ccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                  <  2:3414336/93‑1 (MQ=255)
                                                                       ccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                  >  2:160045/1‑93 (MQ=255)
                                                                       ccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                  >  2:5036837/1‑93 (MQ=255)
                                                                       ccATCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACaga                  >  1:2733410/1‑93 (MQ=255)
                                                                         aTCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAgaga                <  2:8169047/93‑1 (MQ=255)
                                                                         aTCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAgaga                >  2:1459439/1‑93 (MQ=255)
                                                                         aTCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATCCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAgaga                >  2:8781044/1‑93 (MQ=255)
                                                                          tCACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGaa               <  1:1878250/93‑1 (MQ=255)
                                                                          ccaggATGGTGACAGCTGCCGGATTATTGCGGGAAGATATGCGTCTTAAAGCTCCTGTGGGGCACAACTCTCGGCGGTCGGGGGCACAGAGaa               <  2:1339203/89‑1 (MQ=255)
                                                                           cACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAAt              <  2:425829/93‑1 (MQ=255)
                                                                           cACGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAAt              >  1:1541895/1‑93 (MQ=255)
                                                                            aCGATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATa             >  2:8578119/1‑93 (MQ=255)
                                                                              gATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGc           >  2:7185483/1‑93 (MQ=255)
                                                                              gATGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGc           <  1:6405116/93‑1 (MQ=255)
                                                                               aTGGCGACACGTGCCGGATTATTGCGGGAGGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCa          <  2:4868708/93‑1 (MQ=255)
                                                                                ttgCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc         <  1:5191699/91‑1 (MQ=255)
                                                                                tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc         <  2:2304189/93‑1 (MQ=255)
                                                                                tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc         <  2:1519407/93‑1 (MQ=255)
                                                                                tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc         <  1:7814748/93‑1 (MQ=255)
                                                                                tGGCGACACGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCAc         >  1:7545181/1‑93 (MQ=255)
                                                                                       acGTGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAaga  >  1:2525816/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTGGTGCGGGACATTCAGGATTTCGAAGATTTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGA  >  NC_000913/552210‑552389

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