Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,004,812 T→C R448R (CGT→CGC recQ → ATP‑dependent DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,004,8120TC99.6% 764.0 / NA 245R448R (CGT→CGCrecQATP‑dependent DNA helicase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (130/114);  minor base A (0/1);  total (130/115)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 254 total aligned reads displayed.

TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901
                                                                                          |                                                                                         
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         <  1:1455238/93‑1 (MQ=255)
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         >  1:2821934/1‑93 (MQ=255)
tCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCAt                                                                                         <  1:276584/93‑1 (MQ=255)
 cAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATc                                                                                        <  1:804019/93‑1 (MQ=255)
  aGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATcc                                                                                       >  2:1592167/1‑93 (MQ=255)
   gATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcg                                                                                      >  1:2881907/1‑93 (MQ=255)
     ttGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgcg                                                                                    <  2:533566/93‑1 (MQ=255)
     ttGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCcgcg                                                                                    >  1:2843686/1‑93 (MQ=255)
         ccTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTa                                                                                >  2:51219/1‑93 (MQ=255)
          cTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCCTCCGCGACTAt                                                                               >  1:2036718/1‑93 (MQ=255)
          cTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTAt                                                                               >  2:2326412/1‑93 (MQ=255)
          cTTGCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTAt                                                                               <  1:207463/93‑1 (MQ=255)
           tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATg                                                                              >  1:873394/1‑93 (MQ=255)
           tttCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATg                                                                              <  1:3482287/93‑1 (MQ=255)
               caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCa                                                                          >  1:2250776/1‑93 (MQ=255)
               caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCa                                                                          <  1:1512541/93‑1 (MQ=255)
               caccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCa                                                                          >  1:1212325/1‑93 (MQ=255)
                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         <  2:3971305/93‑1 (MQ=255)
                accaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCAt                                                                         >  2:603750/1‑93 (MQ=255)
                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  2:2153190/1‑93 (MQ=255)
                 ccaTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATg                                                                        >  2:3145416/1‑93 (MQ=255)
                  caTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGa                                                                       <  2:2088535/93‑1 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      >  2:2512785/1‑93 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      <  2:2649083/93‑1 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGAc                                                                      <  1:2372632/93‑1 (MQ=255)
                   aTTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGAGCATCCGCGACTATGGTCATGAc                                                                      >  2:2893953/1‑93 (MQ=255)
                    ttGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACa                                                                     >  1:3081329/1‑93 (MQ=255)
                     tGGTCGTGTGAATCAGCGGTTTGGGATAGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACaa                                                                    <  1:751171/93‑1 (MQ=255)
                       gTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAAc                                                                  <  2:3360591/93‑1 (MQ=255)
                        tCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACt                                                                 >  1:549710/1‑93 (MQ=255)
                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               <  2:1330665/93‑1 (MQ=255)
                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               <  1:3971629/93‑1 (MQ=255)
                          gtgtGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGa                                                               <  2:4556903/93‑1 (MQ=255)
                           tgtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaa                                                              >  1:3207243/1‑93 (MQ=255)
                           tgtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaa                                                              >  2:3436429/1‑93 (MQ=255)
                            gtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaaa                                                             <  2:3915002/93‑1 (MQ=255)
                            gtgAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGaaa                                                             <  2:671445/93‑1 (MQ=255)
                                aTCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCt                                                         <  2:932429/93‑1 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        <  1:3273781/93‑1 (MQ=255)
                                 tCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTa                                                        <  1:2525183/93‑1 (MQ=255)
                                  cAGCGGTGTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATCACAAACTGAAAGTCTAt                                                       <  1:3816557/93‑1 (MQ=255)
                                   aGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATg                                                      >  2:3572486/1‑93 (MQ=255)
                                   aGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATg                                                      <  1:4074497/93‑1 (MQ=255)
                                    gCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATgg                                                     >  1:3170674/1‑93 (MQ=255)
                                      ggTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCa                                                   >  2:3097872/1‑93 (MQ=255)
                                       gTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCAt                                                  <  1:1884191/93‑1 (MQ=255)
                                        tttGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATg                                                 <  2:1489581/93‑1 (MQ=255)
                                          tGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATggg                                               >  1:1805204/1‑93 (MQ=255)
                                            ggATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGcc                                             >  2:4656267/1‑93 (MQ=255)
                                             gATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCg                                            >  1:1694090/1‑93 (MQ=255)
                                               tGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTg                                          >  1:4322019/1‑93 (MQ=255)
                                                gggTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGa                                         <  2:1111266/93‑1 (MQ=255)
                                                 ggTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGAt                                        >  2:2360333/1‑93 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       <  2:4596820/93‑1 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       <  2:1068125/93‑1 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       >  2:1253558/1‑93 (MQ=255)
                                                  gTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATa                                       >  1:4529505/1‑93 (MQ=255)
                                                   ttATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaa                                      >  1:3648068/1‑93 (MQ=255)
                                                    tATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaa                                     <  1:1154442/93‑1 (MQ=255)
                                                     aTGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaaa                                    <  2:796274/93‑1 (MQ=255)
                                                     aTGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATaaaa                                    <  1:2499545/93‑1 (MQ=255)
                                                       gtggtgGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGc                                  >  2:1453148/1‑93 (MQ=255)
                                                       gtggtgGAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGc                                  >  2:3807126/1‑93 (MQ=255)
                                                         ggtggAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCa                                <  1:3002706/93‑1 (MQ=255)
                                                          gcggAAGTGATTGGTGCTGTTAATACCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGGATATGGCATGGGCCGTGATAAAGGCCAt                               <  2:2974481/91‑1 (MQ=255)
                                                             gAAGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGaa                            >  1:1015617/1‑93 (MQ=255)
                                                              aaGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAAc                           >  1:3524250/1‑93 (MQ=255)
                                                              aaGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAAc                           <  1:4377984/93‑1 (MQ=255)
                                                               aGTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACCGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACa                          <  1:3312332/93‑1 (MQ=255)
                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         >  1:1547985/1‑93 (MQ=255)
                                                                gTGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAt                         <  2:3281030/93‑1 (MQ=255)
                                                                 tGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAtt                        >  2:2762237/1‑93 (MQ=255)
                                                                 tGATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACAtt                        >  1:3953676/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  2:4345512/1‑93 (MQ=255)
                                                                  gATTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTg                       >  2:1561046/1‑93 (MQ=255)
                                                                   aTTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACCAACTGAAAGTCTATGGCATGGGCCGTGCTAAAAGCCATGAACATTgg                      >  2:1931694/1‑93 (MQ=255)
                                                                   aTTCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTgg                      <  1:3743445/93‑1 (MQ=255)
                                                                    ttCGTGGTGCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTggg                     <  1:1233019/93‑1 (MQ=255)
                                                                       gtggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGa                  <  1:1379831/93‑1 (MQ=255)
                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 >  1:3943988/1‑93 (MQ=255)
                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 >  1:4468747/1‑93 (MQ=255)
                                                                        tggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAg                 >  2:4245938/1‑93 (MQ=255)
                                                                         ggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGc                <  2:749828/93‑1 (MQ=255)
                                                                         ggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGc                >  1:4040203/1‑93 (MQ=255)
                                                                         ggtgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGc                >  2:3839777/1‑93 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              >  1:1203372/1‑93 (MQ=255)
                                                                           tgCTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGt              >  1:3615270/1‑93 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            >  1:4606914/1‑93 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            <  2:1550229/93‑1 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            <  2:4648136/93‑1 (MQ=255)
                                                                             cTAATAACCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGa            >  1:2933418/1‑93 (MQ=255)
                                                                              taataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGAt           <  2:4481949/93‑1 (MQ=255)
                                                                              taataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGAt           >  2:292572/1‑93 (MQ=255)
                                                                              taataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGAt           >  2:2751788/1‑93 (MQ=255)
                                                                               aataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATc          >  1:1027596/1‑93 (MQ=255)
                                                                                ataaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATcc         <  1:959140/93‑1 (MQ=255)
                                                                                 taaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGCGCGTGATCCg        >  1:933122/1‑93 (MQ=255)
                                                                                 taaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCg        >  2:43871/1‑93 (MQ=255)
                                                                                 taaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCg        <  2:170986/93‑1 (MQ=255)
                                                                                  aaCCAGCGCATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGc       <  2:4146220/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCAGATTGCCCTTTCCACCATTGGTCGTGTGAATCAGCGGTTTGGGATGGGTTATGTGGTGGAAGTGATTCGTGGTGCTAATAACCAGCGTATCCGCGACTATGGTCATGACAAACTGAAAGTCTATGGCATGGGCCGTGATAAAAGCCATGAACATTGGGTGAGCGTGATCCGCCAGCT  >  NC_000913/4004722‑4004901

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