Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,731,386 T→C V13A (GTA→GCA)  xylH → XylH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,731,3860TC99.5% 592.7 / ‑5.2 185V13A (GTA→GCA) xylHXylH
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (82/102);  total (83/102)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 191 total aligned reads displayed.

TCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGTACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTCTTT  >  NC_000913/3731296‑3731477
                                                                                          |                                                                                            
tCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCAc                                                                                            >  2:650991/1‑93 (MQ=255)
 cATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCAcc                                                                                           <  1:1091543/93‑1 (MQ=255)
    aaCCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGAc                                                                                        >  1:3629680/1‑93 (MQ=255)
     aCCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACa                                                                                       <  1:3805965/93‑1 (MQ=255)
     aCCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACa                                                                                       >  2:4294311/1‑93 (MQ=255)
       cTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATc                                                                                     <  2:2817173/93‑1 (MQ=255)
       cTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATc                                                                                     <  2:2152495/93‑1 (MQ=255)
       cTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATc                                                                                     <  1:1660007/93‑1 (MQ=255)
        cgACTCCGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCACTCCGTCTGAAGTGAAATTGGCCGCACCGACATcc                                                                                    <  2:1905552/92‑1 (MQ=255)
         gACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCg                                                                                   >  1:2496167/1‑93 (MQ=255)
         gACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCg                                                                                   <  1:827444/93‑1 (MQ=255)
          aCTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCgg                                                                                  <  2:142341/93‑1 (MQ=255)
           cTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCggg                                                                                 >  2:1199763/1‑92 (MQ=255)
           cTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGt                                                                                 <  1:1778851/93‑1 (MQ=255)
           cTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGt                                                                                 <  2:3415304/93‑1 (MQ=255)
            tCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTg                                                                                <  1:4509182/93‑1 (MQ=255)
               ggAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCt                                                                             >  2:4283573/1‑93 (MQ=255)
               ggAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGACGAAAAGCAATCCTTTTTAAGTGAATTTGGCGGAACCGCATTCCGGTGACt                                                                             >  2:3306683/1‑93 (MQ=255)
               ggAGCAGGTGATGGAAGCCGCATTGAGAAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCt                                                                             >  1:320068/1‑93 (MQ=255)
                 aGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtc                                                                           >  2:2512384/1‑93 (MQ=255)
                 aGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtc                                                                           <  1:537324/93‑1 (MQ=255)
                  gCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtct                                                                          >  2:3387007/1‑93 (MQ=255)
                  gCAGGTGATGGAAGCCGCATGGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtct                                                                          <  2:363445/93‑1 (MQ=255)
                   gaGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtctc                                                                         <  2:1364593/92‑1 (MQ=255)
                   cAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtctc                                                                         <  1:4459239/93‑1 (MQ=255)
                   cAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTtctc                                                                         <  1:1123284/93‑1 (MQ=255)
                    aGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTcc                                                                        >  2:1213177/1‑93 (MQ=255)
                    aGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTcc                                                                        <  1:393472/93‑1 (MQ=255)
                    aGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTcc                                                                        >  2:2116271/1‑93 (MQ=255)
                     ggTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGGCGCACCGACATCCGGTGGCTTCTCCg                                                                       >  1:1968592/1‑93 (MQ=255)
                      gTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCgg                                                                     <  2:3780442/93‑1 (MQ=255)
                       tGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCggg                                                                    >  1:360521/1‑93 (MQ=255)
                        gATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGc                                                                   >  1:3038134/1‑93 (MQ=255)
                        gATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGc                                                                   >  1:2422924/1‑93 (MQ=255)
                        gATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGc                                                                   <  2:3605712/93‑1 (MQ=255)
                        gATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGc                                                                   >  1:3936193/1‑93 (MQ=255)
                         aTGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGTTTCTCCGGGCt                                                                  <  1:506794/93‑1 (MQ=255)
                          tGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTg                                                                 >  2:467570/1‑93 (MQ=255)
                               gCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATc                                                            >  2:2680552/1‑93 (MQ=255)
                                  gCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACt                                                         >  2:4143122/1‑93 (MQ=255)
                                       gaggagCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAAtt                                                    >  1:4234310/1‑93 (MQ=255)
                                       gaggagCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAAtt                                                    <  1:43786/93‑1 (MQ=255)
                                         ggagCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTg                                                  <  2:2518202/93‑1 (MQ=255)
                                         ggagCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTg                                                  <  1:1726860/93‑1 (MQ=255)
                                          gagCGAACATCATGTCGAAAAGCCATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGc                                                 >  2:2434281/1‑93 (MQ=255)
                                           agCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCa                                                <  1:404575/93‑1 (MQ=255)
                                           agCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGCCTTCTCCGGGCTGAAATCACTGAATTTGCa                                                <  1:404576/93‑1 (MQ=255)
                                            gCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAg                                               >  1:905291/1‑93 (MQ=255)
                                            gCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAg                                               >  2:1049867/1‑93 (MQ=255)
                                            gCGAACATCATGTCGAAAAGCAATCCGGCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAg                                               <  1:867343/93‑1 (MQ=255)
                                             cGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAgg                                              >  2:3808678/1‑93 (MQ=255)
                                              gAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGt                                             >  2:698632/1‑93 (MQ=255)
                                              gAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGt                                             <  1:552741/93‑1 (MQ=255)
                                                aCATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCt                                           <  2:1915729/93‑1 (MQ=255)
                                                  atcatGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTc                                         >  2:282242/1‑93 (MQ=255)
                                                  atcatGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTc                                         <  2:540038/93‑1 (MQ=255)
                                                   tcatGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTTGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCg                                        >  2:3626571/1‑93 (MQ=255)
                                                   tcatGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCg                                        <  2:2448272/93‑1 (MQ=255)
                                                   tcatGTAAAAAAGCAATCCGTCTGAAGTGAAATTGCCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCg                                        <  1:880428/93‑1 (MQ=255)
                                                    catGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGt                                       <  2:1007069/93‑1 (MQ=255)
                                                     atGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGtg                                      <  1:3415012/93‑1 (MQ=255)
                                                     atGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGtg                                      >  1:3817565/1‑93 (MQ=255)
                                                      tGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGtga                                     >  2:2860264/1‑93 (MQ=255)
                                                      tGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGtga                                     <  1:2189530/93‑1 (MQ=255)
                                                       gTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGAGGCTGAATCACTGAATTTGCAGGTCTTCGtgat                                    >  2:4386855/1‑93 (MQ=255)
                                                            aaaGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTg                               >  1:846022/1‑93 (MQ=255)
                                                            aaaGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTg                               <  1:2045202/93‑1 (MQ=255)
                                                             aaGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGc                              >  2:1739697/1‑93 (MQ=255)
                                                             aaGCAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGc                              >  1:3112521/1‑93 (MQ=255)
                                                                cAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGc                           <  1:258316/93‑1 (MQ=255)
                                                                cAATCCGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGc                           >  1:2846760/1‑93 (MQ=255)
                                                                    ccGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatc                       <  2:3066401/93‑1 (MQ=255)
                                                                     cGTCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatca                      <  2:2775374/93‑1 (MQ=255)
                                                                       tCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatcatc                    <  1:214124/93‑1 (MQ=255)
                                                                       tCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatcatc                    <  2:4545550/93‑1 (MQ=255)
                                                                       tCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatcatc                    <  1:580943/93‑1 (MQ=255)
                                                                       tCTGAAGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTatcatc                    <  1:4700097/93‑1 (MQ=255)
                                                                         cgAAGTGAAATTGGCCGTATCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGc                  <  1:2689912/92‑1 (MQ=255)
                                                                           aaGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCaa                <  1:1480059/93‑1 (MQ=255)
                                                                           aaGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCaa                >  1:2051859/1‑93 (MQ=255)
                                                                            aGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCaac               >  2:3933743/1‑92 (MQ=255)
                                                                            aGTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAAt               >  1:1762262/1‑93 (MQ=255)
                                                                             gTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATc              >  1:3317173/1‑93 (MQ=255)
                                                                             gTGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATc              >  2:157992/1‑93 (MQ=255)
                                                                             gTGAAATTGGCCGCACCGACATCCGGCGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATc              <  1:3256935/93‑1 (MQ=255)
                                                                              tGAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCa             <  1:961328/93‑1 (MQ=255)
                                                                               gAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCAt            <  2:4160490/93‑1 (MQ=255)
                                                                               gAAATTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCAt            >  2:4444585/1‑93 (MQ=255)
                                                                                aaaTTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATg           <  1:133784/93‑1 (MQ=255)
                                                                                aaaTTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATg           >  2:315739/1‑93 (MQ=255)
                                                                                 aaTTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGc          <  1:4193107/93‑1 (MQ=255)
                                                                                 aaTTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGc          <  1:3089013/93‑1 (MQ=255)
                                                                                  aTTGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCt         >  1:1041222/1‑93 (MQ=255)
                                                                                   ttGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTg        >  2:2714026/1‑93 (MQ=255)
                                                                                    tGGCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGt       >  2:1226213/1‑93 (MQ=255)
                                                                                      gCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTc     >  2:1487741/1‑93 (MQ=255)
                                                                                      gCCGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTc     <  2:4727230/93‑1 (MQ=255)
                                                                                       ccGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTCt    <  1:4282765/93‑1 (MQ=255)
                                                                                       ccGCACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTCt    <  1:1205111/93‑1 (MQ=255)
                                                                                         ctACCGGCACCCGGAGGCTTCTCTGGGCTGAAATCACTGTAATTGTAGGTCTTCGTGATGATTGCTGCTATGATCGCGATCCTGCTGTTCttt  >  2:3669324/2‑93 (MQ=255)
                                                                                          |                                                                                            
TCATAACCTGACTCAGGAGCAGGTGATGGAAGCCGCATTGAGGAGCGAACATCATGTCGAAAAGCAATCCGTCTGAAGTGAAATTGGCCGTACCGACATCCGGTGGCTTCTCCGGGCTGAAATCACTGAATTTGCAGGTCTTCGTGATGATTGCAGCTATCATCGCAATCATGCTGTTCTTT  >  NC_000913/3731296‑3731477

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