Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,259,162 A→G intergenic (+277/‑86) yjbM → / → yjbN hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,259,1620AG99.7% 1203.3 / ‑5.5 374intergenic (+277/‑86)yjbM/yjbNhypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (187/186);  total (188/186)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 406 total aligned reads displayed.

TTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGC  >  NC_000913/4259072‑4259251
                                                                                          |                                                                                         
ttattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTatggtatggt                                                                                         <  2:344829/93‑1 (MQ=255)
ttattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTatggtatggt                                                                                         >  2:7443303/1‑93 (MQ=255)
 tattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCATTGAGTCAGGCGCTTCCTGTTGTTAtggtatggta                                                                                        >  2:3143767/1‑93 (MQ=255)
 tattTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCCTCCTGTTGTTAtggtatggta                                                                                        >  1:1609606/1‑93 (MQ=255)
  attTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTAc                                                                                       >  1:7463046/1‑93 (MQ=255)
   tttCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTAcc                                                                                      >  2:6940479/1‑93 (MQ=255)
         tCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGTGCTTCCTGTTGTTATGGTATGGTACCCCGCTa                                                                                <  1:6473966/93‑1 (MQ=255)
         tCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTa                                                                                <  2:4465579/93‑1 (MQ=255)
         tCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTa                                                                                >  2:1982724/1‑93 (MQ=255)
          cACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAt                                                                               <  1:5899158/93‑1 (MQ=255)
          cACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAt                                                                               <  1:4184109/93‑1 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              <  2:3630610/93‑1 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              >  1:7143617/1‑93 (MQ=255)
           aCTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTAtt                                                                              >  2:3903795/1‑93 (MQ=255)
            cTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTg                                                                             <  2:1516148/93‑1 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            >  1:4948629/1‑93 (MQ=255)
             ttCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGa                                                                            <  1:3554937/93‑1 (MQ=255)
             ttCGGGAGAGCTTCATGCCTTGAGAGGAACGCAAATTTCTCATGCAATGAGACAGGCGCTGCCCGGTGTTAGGGTATGGTACCCCGCGACTGa                                                                            >  1:7009027/1‑93 (MQ=255)
                gTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGTTATGGTACCCCGCTATTGAGcc                                                                         >  1:5058046/1‑93 (MQ=255)
                 tGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcct                                                                        >  1:4868872/1‑93 (MQ=255)
                 tGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcct                                                                        >  2:872514/1‑93 (MQ=255)
                  gagaGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcctc                                                                       >  1:6093165/1‑93 (MQ=255)
                  gagaGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcctc                                                                       >  2:7559526/1‑93 (MQ=255)
                    gagCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcctcct                                                                     >  2:731013/1‑93 (MQ=255)
                    gagCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGcctcct                                                                     >  1:3379165/1‑93 (MQ=255)
                      gCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGa                                                                   >  1:4208000/1‑93 (MQ=255)
                      gCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGa                                                                   >  1:3076702/1‑93 (MQ=255)
                       cTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGaa                                                                  >  2:6063905/1‑93 (MQ=255)
                       cTTCATGCCTTGACAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGaa                                                                  <  2:5721684/93‑1 (MQ=255)
                           aTGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGt                                                              <  1:4447319/93‑1 (MQ=255)
                           aTGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGt                                                              >  1:7325697/1‑93 (MQ=255)
                             gCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGa                                                            >  2:6898783/1‑93 (MQ=255)
                             gCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGa                                                            <  1:3174083/93‑1 (MQ=255)
                              ccTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACGGGCGCTTCCTGTTGTTATGGTATGGTACCACGCTTTTGAGCCTCCTGAATAGTGAt                                                           <  2:7027882/93‑1 (MQ=255)
                              ccTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGAt                                                           <  1:1002826/93‑1 (MQ=255)
                              ccTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGAt                                                           <  1:7899972/93‑1 (MQ=255)
                              ccTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGAt                                                           <  2:4935638/93‑1 (MQ=255)
                                  gagaGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTg                                                       <  1:6008593/93‑1 (MQ=255)
                                    gagGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGaa                                                     >  2:3786752/1‑93 (MQ=255)
                                     agGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAat                                                    <  2:2167639/93‑1 (MQ=255)
                                     agGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAat                                                    >  2:4356069/1‑93 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTCCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                <  1:260702/93‑1 (MQ=255)
                                         tctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataaca                                                >  2:3298732/1‑93 (MQ=255)
                                          ctcAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacat                                               >  2:1399834/1‑93 (MQ=255)
                                            cAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataa                                             >  2:6004888/1‑93 (MQ=255)
                                             aaTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataac                                            <  2:1134943/93‑1 (MQ=255)
                                             aaTTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAataacataac                                            <  1:5789870/93‑1 (MQ=255)
                                                tttCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAAcccc                                         >  1:1170336/1‑93 (MQ=255)
                                                 ttCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCa                                        >  1:6609492/1‑93 (MQ=255)
                                                 ttCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCa                                        <  1:2449766/93‑1 (MQ=255)
                                                 ttCTTTTGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCa                                        >  1:2120119/1‑93 (MQ=255)
                                                  tCTTTTGCAATGAGACAGGCGCCTCCTGCTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCAt                                       <  1:1872600/93‑1 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     <  1:6212907/93‑1 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     >  2:65796/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     >  2:6492442/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     >  1:2764064/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     >  2:534332/1‑93 (MQ=255)
                                                    ttttGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGa                                     <  2:4658015/93‑1 (MQ=255)
                                                     tttGCAATGAGACAGGCGCCTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGat                                    <  1:1231302/93‑1 (MQ=255)
                                                      ttGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGata                                   <  2:3670294/93‑1 (MQ=255)
                                                       tGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatat                                  <  1:7232054/93‑1 (MQ=255)
                                                        gCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata                                 >  1:1287629/1‑93 (MQ=255)
                                                        gCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata                                 <  1:413431/93‑1 (MQ=255)
                                                        gCAATGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGatata                                 <  2:2779371/93‑1 (MQ=255)
                                                          aaTGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATc                               <  1:4046079/93‑1 (MQ=255)
                                                           aTGAGCCGGGCGCTTCTTGTTGTAATGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGa                             <  1:3920003/93‑1 (MQ=255)
                                                            tGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTTAATAACATAACCCCATGATATATCGa                             >  1:845204/1‑93 (MQ=255)
                                                            tGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGa                             <  1:5295494/93‑1 (MQ=255)
                                                            tGAGACAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGa                             >  2:6762251/1‑93 (MQ=255)
                                                             gagaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGAt                            >  2:6327421/1‑93 (MQ=255)
                                                             gagaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGAt                            <  2:5120115/93‑1 (MQ=255)
                                                              agaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATa                           <  2:7093398/93‑1 (MQ=255)
                                                               gaCAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaa                          >  2:4073720/1‑93 (MQ=255)
                                                                 cAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaaaa                        <  1:3781822/93‑1 (MQ=255)
                                                                 cAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaaaa                        <  1:7489882/93‑1 (MQ=255)
                                                                 cAGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATaaaa                        <  1:2209443/93‑1 (MQ=255)
                                                                  aGGCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaat                       <  1:3857276/93‑1 (MQ=255)
                                                                   ggCGCTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGAGATGCTGAATAACATAACCCCATGATATATCGATAAaata                      >  1:5578969/1‑93 (MQ=255)
                                                                    gcgcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaataa                     <  2:7825272/93‑1 (MQ=255)
                                                                    gcgcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaataa                     <  2:7187954/93‑1 (MQ=255)
                                                                     cgcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaataat                    <  2:4879476/93‑1 (MQ=255)
                                                                     cgcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaataat                    <  2:254928/93‑1 (MQ=255)
                                                                     cgcTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAaataat                    <  2:1826942/93‑1 (MQ=255)
                                                                       cTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtct                  >  1:7757543/1‑93 (MQ=255)
                                                                       cTTCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtct                  >  1:1771193/1‑93 (MQ=255)
                                                                        ttCCTGTTGTTATGGTATGGTCCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtctc                 <  1:5450355/93‑1 (MQ=255)
                                                                        ttCCTGTTGTTATGGTATGGTACCCCGCTATTGCGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAAtctc                 >  1:6935265/1‑93 (MQ=255)
                                                                         tCCTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATctct                >  2:305041/1‑93 (MQ=255)
                                                                          tcTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTa               <  1:6855366/92‑1 (MQ=255)
                                                                          ccTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTa               <  1:6653099/93‑1 (MQ=255)
                                                                          ccTGTTGTTATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTa               >  1:4535977/1‑93 (MQ=255)
                                                                             gttgttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACAt            >  2:2960826/1‑93 (MQ=255)
                                                                              ttgttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACAtt           >  1:6987960/1‑93 (MQ=255)
                                                                               tgttATGGTATTGTACCCCGGTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACAttt          <  1:1405854/93‑1 (MQ=255)
                                                                                 ttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGa        <  1:4698313/93‑1 (MQ=255)
                                                                                 ttATGGTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGa        <  2:1328433/93‑1 (MQ=255)
                                                                                   atggtatggtACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaaa      <  2:550604/93‑1 (MQ=255)
                                                                                    tggtatggtaCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaaaa     <  2:6716382/93‑1 (MQ=255)
                                                                                    tggtatggtaCCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGaaaa     <  2:2389731/93‑1 (MQ=255)
                                                                                      gTATGGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATg   >  2:1700431/1‑93 (MQ=255)
                                                                                          |                                                                                         
TTATTTCCATCACTTCGTGAGAGCTTCATGCCTTGAGAGGATCTCAATTTTCTTTTGCAATGAGACAGGCGCTTCCTGTTGTTATGGTATAGTACCCCGCTATTGAGCCTCCTGAATAGTGATGCTGAATAACATAACCCCATGATATATCGATAAAATAATCTCTACATTTGAAAATGC  >  NC_000913/4259072‑4259251

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