Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,427,075 T→C V68A (GTC→GCC)  b3279 → putative transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,427,0750TC96.7% 718.6 / 8.6 243V68A (GTC→GCC) b3279putative transferase
Reads supporting (aligned to +/- strand):  ref base T (3/5);  new base C (104/130);  total (108/135)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

Only 100 of 244 total aligned reads displayed.

GATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGCTGC  >  NC_000913/3426986‑3427166
                                                                                         |                                                                                           
 aTCTGGCGGCTCGTTGTGATTCGGGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACt                                                                                         >  2:139684/1‑93 (MQ=255)
 aTCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGTAGTATGTTGCATGCCACt                                                                                         >  1:3392647/1‑93 (MQ=255)
 aTCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACt                                                                                         >  2:3632013/1‑93 (MQ=255)
    tGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCAt                                                                                      <  2:4285198/93‑1 (MQ=255)
     ggCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATa                                                                                     >  2:2306414/1‑93 (MQ=255)
         gCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATc                                                                                 <  2:1483450/93‑1 (MQ=255)
         gCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATc                                                                                 <  1:466091/93‑1 (MQ=255)
          cTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATcc                                                                                <  2:2624016/93‑1 (MQ=255)
           tCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCt                                                                               >  1:4503830/1‑93 (MQ=255)
           tCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCt                                                                               <  1:3002450/93‑1 (MQ=255)
           tCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCt                                                                               >  1:1739992/1‑93 (MQ=255)
            cGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTc                                                                              <  1:1629569/93‑1 (MQ=255)
            cGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTc                                                                              <  2:3791524/93‑1 (MQ=255)
                gtgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTAc                                                                          >  2:1423665/1‑93 (MQ=255)
                gtgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTAc                                                                          <  2:403284/93‑1 (MQ=255)
                gtgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTAc                                                                          >  1:2491749/1‑93 (MQ=255)
                 tgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACa                                                                         >  2:3719798/1‑93 (MQ=255)
                 tgATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACa                                                                         <  1:510084/93‑1 (MQ=255)
                  gATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACaa                                                                        >  1:2851270/1‑93 (MQ=255)
                   aTTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAAc                                                                       >  1:2715896/1‑93 (MQ=255)
                    ttCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAAcc                                                                      >  1:1018241/1‑93 (MQ=255)
                    ttCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAAcc                                                                      >  1:1847671/1‑93 (MQ=255)
                    ttCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAAcc                                                                      <  2:2825938/93‑1 (MQ=255)
                       gTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAg                                                                   >  1:1596846/1‑93 (MQ=255)
                       gTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAg                                                                   <  2:3198561/93‑1 (MQ=255)
                        tGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGa                                                                  >  2:2970427/1‑93 (MQ=255)
                        tGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGa                                                                  <  2:4495020/93‑1 (MQ=255)
                         ggAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGAt                                                                 >  1:3041493/1‑93 (MQ=255)
                         ggAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGAt                                                                 >  1:1938431/1‑93 (MQ=255)
                         ggAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGAt                                                                 <  1:96578/93‑1 (MQ=255)
                         ggAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTCGCATGCCACTCATAAATCCTCGTACAACCCAGAt                                                                 <  2:2976671/93‑1 (MQ=255)
                          gagaTGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATg                                                                >  2:4353647/1‑93 (MQ=255)
                          gagaTGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATg                                                                >  1:3191130/1‑93 (MQ=255)
                             aTGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCa                                                             >  1:1196679/1‑93 (MQ=255)
                              tGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCaa                                                            >  1:4140238/1‑93 (MQ=255)
                              tGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCaa                                                            >  1:4352850/1‑93 (MQ=255)
                               gTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAAc                                                           <  2:726573/93‑1 (MQ=255)
                                tACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAAcc                                                          <  2:2175740/93‑1 (MQ=255)
                                 aCATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAAccc                                                         <  2:4044036/93‑1 (MQ=255)
                                  cATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCa                                                        <  2:576274/93‑1 (MQ=255)
                                   aTTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCAt                                                       >  1:3163732/1‑93 (MQ=255)
                                    ttATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCAtt                                                      <  1:736457/93‑1 (MQ=255)
                                    ttATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCAtt                                                      <  1:1803403/93‑1 (MQ=255)
                                      aTGTACAGATCGGAGCGCGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTaa                                                    <  1:43039/93‑1 (MQ=255)
                                      aTGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTaa                                                    >  2:3262650/1‑93 (MQ=255)
                                       tGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAAc                                                   >  1:2550017/1‑93 (MQ=255)
                                       tGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAAc                                                   <  1:2251868/93‑1 (MQ=255)
                                       tGTACAGATCGGAGCACGCACCAAAATCCCGGATGGCAGGATGTTGCATGCCACTCATAAATCCTCGTACAAACCAGATGGCAACCCATTaaa                                                   >  2:1664724/1‑92 (MQ=255)
                                        gTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAAcc                                                  <  2:1947254/93‑1 (MQ=255)
                                        gTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAAcc                                                  <  1:2382109/93‑1 (MQ=255)
                                         tACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCa                                                 <  2:235804/93‑1 (MQ=255)
                                         tACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCa                                                 <  1:673135/93‑1 (MQ=255)
                                         tACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCa                                                 >  2:126176/1‑93 (MQ=255)
                                         tACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCa                                                 <  2:386309/93‑1 (MQ=255)
                                           cAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCACCCCATTAACCAtt                                               <  1:2002802/93‑1 (MQ=255)
                                            aGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTg                                              <  1:933901/93‑1 (MQ=255)
                                              aTCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGc                                            >  1:879140/1‑93 (MQ=255)
                                                 ggAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGaa                                         <  2:5399/93‑1 (MQ=255)
                                                 ggAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGaa                                         <  1:595999/93‑1 (MQ=255)
                                                  gAGGACGCACCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAg                                        >  2:3961485/1‑93 (MQ=255)
                                                   aggacggacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAGCCCATTAACCATTGGCGAAGa                                       <  2:4070495/90‑1 (MQ=255)
                                                    gcacgcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGAt                                      >  1:2671324/1‑93 (MQ=255)
                                                     cacgcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATg                                     >  2:3502552/1‑93 (MQ=255)
                                                       cgcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTc                                   <  1:2043983/93‑1 (MQ=255)
                                                        gcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCa                                  >  2:795544/1‑93 (MQ=255)
                                                        gcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCa                                  <  2:1570200/93‑1 (MQ=255)
                                                        gcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCa                                  >  1:2939808/1‑93 (MQ=255)
                                                        gcacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGAAGGCAACCCATTAACCATTGGCGAAGATGTCa                                  <  1:3300463/93‑1 (MQ=255)
                                                         cacCAATATCCAGGATGGCTGTATGTTGCATGCCACTCATACATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCAc                                 <  1:4205057/93‑1 (MQ=255)
                                                         cacCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCAc                                 >  1:2327312/1‑93 (MQ=255)
                                                          acCAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACt                                <  1:2061235/93‑1 (MQ=255)
                                                           ccAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTg                               >  1:1515376/1‑93 (MQ=255)
                                                            cAATATCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGt                              >  2:1129708/1‑93 (MQ=255)
                                                              atatCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTg                            >  1:3580983/1‑93 (MQ=255)
                                                              atatCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTg                            <  2:221701/93‑1 (MQ=255)
                                                                atCCAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGt                          <  1:281558/93‑1 (MQ=255)
                                                                  ccAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTca                        <  1:46885/93‑1 (MQ=255)
                                                                   cAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTcac                       <  1:2086891/93‑1 (MQ=255)
                                                                   cAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTcac                       <  2:3647125/93‑1 (MQ=255)
                                                                   cAGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTAGCGAAGATGTCACTGTTGGTcac                       >  1:358729/1‑93 (MQ=255)
                                                                    aGGATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTcaca                      >  2:4215590/1‑93 (MQ=255)
                                                                     ggATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACaa                     >  1:765506/1‑93 (MQ=255)
                                                                     ggATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACaa                     <  1:4635922/93‑1 (MQ=255)
                                                                     ggATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACaa                     <  2:2141824/93‑1 (MQ=255)
                                                                      gATGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAg                    >  1:3799531/1‑93 (MQ=255)
                                                                       aTGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAgg                   <  2:1303126/93‑1 (MQ=255)
                                                                        tGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAggg                  >  2:3855347/1‑92 (MQ=255)
                                                                        tGGCAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGt                  <  1:2278073/93‑1 (MQ=255)
                                                                           cAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGAt               >  2:4105214/1‑93 (MQ=255)
                                                                           cAGTATGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGAt               <  2:4663710/93‑1 (MQ=255)
                                                                               aTGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTc           >  2:3856559/1‑93 (MQ=255)
                                                                               aTGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTc           >  2:416602/1‑93 (MQ=255)
                                                                                tGTTTCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTcc          <  1:3919264/93‑1 (MQ=255)
                                                                                tGTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTcc          <  1:4077675/93‑1 (MQ=255)
                                                                                 gTTGCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCa         <  1:644507/93‑1 (MQ=255)
                                                                                    gCATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACgg      <  2:1781799/93‑1 (MQ=255)
                                                                                     cATGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGc     <  1:4425668/93‑1 (MQ=255)
                                                                                      aTGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGCt    <  2:230701/93‑1 (MQ=255)
                                                                                      aTGCCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGCt    <  2:2889098/93‑1 (MQ=255)
                                                                                        gTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGCTGc  <  2:2011609/93‑1 (MQ=255)
                                                                                         |                                                                                           
GATCTGGCCGCTCGTTGTGATTCGTGGAGATGTACATTATGTACAGATCGGAGCACGCACCAATATCCAGGATGGCAGTATGTTGCATGTCACTCATAAATCCTCGTACAACCCAGATGGCAACCCATTAACCATTGGCGAAGATGTCACTGTTGGTCACAAGGTGATGCTCCACGGCTGC  >  NC_000913/3426986‑3427166

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