Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 736,480 A→G T52A (ACG→GCG)  ybfL → putative receptor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913736,4800AG86.7% 463.3 / 75.5 227T52A (ACG→GCG) ybfLputative receptor protein
Reads supporting (aligned to +/- strand):  ref base A (18/12);  new base G (84/112);  total (102/125)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 229 total aligned reads displayed.

CTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGT  >  NC_000913/736392‑736572
                                                                                        |                                                                                            
ctctTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGa                                                                                          >  1:2311372/1‑93 (MQ=18)
 tctTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGAt                                                                                         <  1:6165986/93‑1 (MQ=18)
  ctTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATa                                                                                        >  2:2917596/1‑93 (MQ=21)
   ttATGACAAGAGTCGCCGCAGGGGAGCTATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATaa                                                                                       <  2:3272298/93‑1 (MQ=255)
   ttATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCACGGCGGATaa                                                                                       >  1:2655244/1‑93 (MQ=255)
     aTGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGa                                                                                     >  1:3046267/1‑93 (MQ=37)
      tGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGaa                                                                                    >  1:5766533/1‑93 (MQ=37)
      tGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGaa                                                                                    >  1:3301691/1‑93 (MQ=37)
       gACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGaaa                                                                                   <  1:5426164/93‑1 (MQ=255)
       gACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGaaa                                                                                   <  1:1881184/93‑1 (MQ=255)
       gACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGaaa                                                                                   >  1:5616807/1‑93 (MQ=37)
             agTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAAt                                                                             >  1:3639044/1‑93 (MQ=255)
              gTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATg                                                                            >  2:5305093/1‑93 (MQ=255)
              gTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATg                                                                            >  2:6340723/1‑93 (MQ=255)
              gTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATg                                                                            >  1:794496/1‑93 (MQ=255)
              gTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATg                                                                            <  2:6479865/93‑1 (MQ=255)
               tCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATga                                                                           >  1:6214256/1‑93 (MQ=255)
                cgccgcAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATgag                                                                          <  2:1379404/93‑1 (MQ=255)
                    gcAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATca                                                                      <  1:4810014/93‑1 (MQ=255)
                    gcAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATca                                                                      >  1:1018921/1‑93 (MQ=255)
                    gcAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATca                                                                      >  1:4862703/1‑93 (MQ=38)
                     cAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATcac                                                                     <  2:5793480/93‑1 (MQ=255)
                     cAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATcac                                                                     <  2:4415986/93‑1 (MQ=255)
                      aGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATcaca                                                                    <  2:2653333/93‑1 (MQ=255)
                      aGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATcaca                                                                    >  2:1934598/1‑93 (MQ=25)
                      aGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATcaca                                                                    >  1:396198/1‑93 (MQ=35)
                       ggggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAg                                                                   <  1:401885/93‑1 (MQ=21)
                       ggggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAg                                                                   >  2:3063565/1‑93 (MQ=21)
                        gggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGc                                                                  >  1:3612341/1‑93 (MQ=18)
                        gggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGc                                                                  <  1:3512417/93‑1 (MQ=18)
                        gggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGc                                                                  <  1:339180/93‑1 (MQ=18)
                         ggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCt                                                                 >  2:4544411/1‑93 (MQ=18)
                         ggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCt                                                                 >  1:453432/1‑93 (MQ=18)
                         ggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCt                                                                 <  2:4986988/93‑1 (MQ=18)
                         ggAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCt                                                                 >  2:3414562/1‑93 (MQ=18)
                            gCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATc                                                              <  2:5392547/93‑1 (MQ=18)
                             cGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATcc                                                             <  2:5700113/93‑1 (MQ=21)
                              gATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATccc                                                            >  2:1422324/1‑93 (MQ=18)
                              gATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATccc                                                            <  2:3805900/93‑1 (MQ=21)
                               aTTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCt                                                           >  2:3372437/1‑93 (MQ=18)
                               aTTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCt                                                           <  1:6620234/93‑1 (MQ=18)
                               aTTCATGTCATTAGTGCGTGCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCt                                                           <  2:3774324/93‑1 (MQ=16)
                                ttCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTg                                                          >  1:5084363/1‑93 (MQ=18)
                                 tCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCTCTGa                                                         >  1:6415149/1‑93 (MQ=16)
                                      tCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAActtc                                                    >  2:3727313/1‑93 (MQ=21)
                                      tCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAActtc                                                    <  1:2481777/93‑1 (MQ=21)
                                      tCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAActtc                                                    <  2:1327004/93‑1 (MQ=21)
                                        aTTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAActtctt                                                  <  1:138364/93‑1 (MQ=34)
                                          tAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTaa                                                <  2:2382537/93‑1 (MQ=21)
                                           aGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCTGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAAc                                               <  1:5543575/93‑1 (MQ=16)
                                           aGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAAc                                               <  1:4243923/93‑1 (MQ=21)
                                            gTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACa                                              >  2:6110513/1‑93 (MQ=21)
                                             tGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACAt                                             >  2:914494/1‑93 (MQ=18)
                                                gTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCt                                          >  1:5874643/1‑93 (MQ=21)
                                                 ttCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTg                                         >  2:4744624/1‑93 (MQ=18)
                                                 ttCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTg                                         <  2:4907811/93‑1 (MQ=21)
                                                  tctcAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTgg                                        <  2:3177125/93‑1 (MQ=18)
                                                    tcAACAATGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGat                                      >  1:5844790/1‑93 (MQ=21)
                                                      aacaaTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGatat                                    <  1:1515159/93‑1 (MQ=18)
                                                       acaaTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATAtt                                   >  2:5529832/1‑93 (MQ=18)
                                                       acaaTGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATAtt                                   >  2:4003057/1‑93 (MQ=33)
                                                        caaTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTa                                  >  2:1254197/1‑93 (MQ=21)
                                                        caaTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTa                                  >  2:3355507/1‑93 (MQ=21)
                                                         aaTGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTaa                                 <  2:5607448/93‑1 (MQ=34)
                                                          aTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTaaa                                >  2:4825605/1‑93 (MQ=21)
                                                          aTGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTaaa                                <  1:6269648/93‑1 (MQ=18)
                                                           tGCACAGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAg                               <  2:2459263/93‑1 (MQ=18)
                                                             cacaGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAGCATGCTGGATATTAAAGGa                             <  1:6445103/93‑1 (MQ=255)
                                                              acaGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaa                            <  2:5029358/93‑1 (MQ=18)
                                                                aGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaaaa                          <  1:3544608/93‑1 (MQ=18)
                                                                aGTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaaaa                          <  1:5352794/93‑1 (MQ=18)
                                                                 gTCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGaaaaa                         <  1:1195391/93‑1 (MQ=18)
                                                                  tCTGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAat                        <  1:4151597/93‑1 (MQ=18)
                                                                   cTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAatc                       <  1:6573162/93‑1 (MQ=33)
                                                                    tGGTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAatca                      <  1:2153868/93‑1 (MQ=18)
                                                                     ggTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAatcat                     <  2:2809656/93‑1 (MQ=18)
                                                                      gTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAatcatc                    >  1:6038137/1‑93 (MQ=21)
                                                                      gTCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAatcatc                    <  2:2411973/93‑1 (MQ=18)
                                                                       tCATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAAtcatca                   <  1:1746176/93‑1 (MQ=18)
                                                                       tCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAAtcatca                   <  1:4065286/93‑1 (MQ=33)
                                                                        cATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCGCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATcac                  >  1:1250277/1‑93 (MQ=255)
                                                                        cATCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATcac                  <  1:4794445/93‑1 (MQ=21)
                                                                          tCGGGCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACaa                >  1:4998391/1‑93 (MQ=255)
                                                                             ggCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTg             <  2:1734109/93‑1 (MQ=255)
                                                                             ggCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTg             >  2:1870906/1‑93 (MQ=255)
                                                                             ggCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACGTGCTGGATATTAAAGGAAAAATCATCACAACTg             >  2:1777334/1‑93 (MQ=255)
                                                                              gCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGa            <  2:4425567/93‑1 (MQ=255)
                                                                              gCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGa            >  2:2162631/1‑93 (MQ=255)
                                                                              gCAGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTACCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGa            <  2:3127879/93‑1 (MQ=255)
                                                                               cAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGAt           <  2:2876293/93‑1 (MQ=37)
                                                                                aGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATg          >  2:6646134/1‑93 (MQ=37)
                                                                                aGATCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATg          <  2:1474264/93‑1 (MQ=255)
                                                                                  aTCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCg        >  2:3866814/1‑93 (MQ=37)
                                                                                   tCAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGa       <  2:4710686/93‑1 (MQ=37)
                                                                                    cAAGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGAt      >  1:5118160/1‑93 (MQ=37)
                                                                                     aaGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATg     >  2:5152515/1‑93 (MQ=37)
                                                                                     aaGGCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATg     <  2:855235/93‑1 (MQ=37)
                                                                                     aaGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATg     <  1:6234149/93‑1 (MQ=37)
                                                                                      aGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATgg    >  1:872786/1‑93 (MQ=37)
                                                                                        aCGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGt  <  1:434529/93‑1 (MQ=38)
                                                                                        |                                                                                            
CTCTTATGACAAGAGTCGCCGCAGGGGAGCGATTCATGTCATTAGTGCGTTCTCAACAATGCACAGTCTGGTCATCGGGCAGATCAAGACGGATAAGAAATCCAATGAGATCACAGCTATCCCTGAACTTCTTAACATGCTGGATATTAAAGGAAAAATCATCACAACTGATGCGATGGGT  >  NC_000913/736392‑736572

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