Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,562,886 A→G V299A (GTT→GCT)  b1489 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,562,8860AG100.0% 779.8 / NA 238V299A (GTT→GCT) b1489hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (99/139);  total (99/139)

Only 100 of 244 total aligned reads displayed.

GCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGA  >  NC_000913/1562796‑1562975
                                                                                          |                                                                                         
gCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCa                                                                                         <  1:4378253/93‑1 (MQ=255)
gCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCa                                                                                         >  2:5930250/1‑93 (MQ=255)
 cTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCAt                                                                                        <  1:3503672/93‑1 (MQ=255)
 cTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCAt                                                                                        <  2:1014965/93‑1 (MQ=255)
  tGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATg                                                                                       >  2:4180778/1‑93 (MQ=255)
   gACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGc                                                                                      >  1:4692754/1‑93 (MQ=255)
   gACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGc                                                                                      >  2:5243693/1‑93 (MQ=255)
   gACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGc                                                                                      >  2:6398912/1‑93 (MQ=255)
    aCCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCg                                                                                     <  1:1308217/93‑1 (MQ=255)
    aCCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCg                                                                                     <  1:1584105/93‑1 (MQ=255)
     ccAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGa                                                                                    >  2:4609641/1‑93 (MQ=255)
     ccAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGa                                                                                    <  1:5096988/93‑1 (MQ=255)
      cAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGAt                                                                                   <  1:5298968/93‑1 (MQ=255)
      cAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGAt                                                                                   >  2:2403644/1‑93 (MQ=255)
        gCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTc                                                                                 >  1:2539199/1‑93 (MQ=255)
          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               <  2:3501897/93‑1 (MQ=255)
          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               <  2:3494204/93‑1 (MQ=255)
          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               <  1:3582584/93‑1 (MQ=255)
          tttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGt                                                                               >  2:5061629/1‑93 (MQ=255)
           ttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGtt                                                                              <  2:5335859/93‑1 (MQ=255)
           ttGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGtt                                                                              >  1:123513/1‑93 (MQ=255)
             gcgcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACGGCGAAGCACGCGATTCGttta                                                                            >  1:185883/1‑91 (MQ=255)
              cgcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAg                                                                           <  2:489031/93‑1 (MQ=255)
               gcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGt                                                                          <  2:4376396/93‑1 (MQ=255)
               gcATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGt                                                                          <  1:5367606/93‑1 (MQ=255)
                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       <  1:294815/93‑1 (MQ=255)
                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       <  2:4767186/93‑1 (MQ=255)
                  ttttGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACa                                                                       <  1:739163/93‑1 (MQ=255)
                      gAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGAtt                                                                   >  2:1800296/1‑93 (MQ=255)
                      gAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGAtt                                                                   <  1:4509232/93‑1 (MQ=255)
                      gAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGAtt                                                                   <  2:5689961/93‑1 (MQ=255)
                          ttCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGAt                                                               <  2:3157460/93‑1 (MQ=255)
                           tCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATg                                                              <  1:463676/93‑1 (MQ=255)
                           tCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATg                                                              <  2:5815310/93‑1 (MQ=255)
                           tCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATg                                                              <  2:2929071/93‑1 (MQ=255)
                            cTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGt                                                             <  2:1168823/93‑1 (MQ=255)
                             tGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGtt                                                            >  2:1045684/1‑93 (MQ=255)
                               cACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTAc                                                          >  1:3309704/1‑93 (MQ=255)
                               cACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTAc                                                          <  2:1939795/93‑1 (MQ=255)
                                cccGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACg                                                         <  1:3897942/92‑1 (MQ=255)
                                 ccGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGa                                                        <  2:5976556/93‑1 (MQ=255)
                                 ccGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGa                                                        >  1:1745099/1‑93 (MQ=255)
                                   gTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACa                                                      >  1:5012605/1‑93 (MQ=255)
                                    tGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaa                                                     >  2:4539649/1‑93 (MQ=255)
                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                    >  2:4468788/1‑93 (MQ=255)
                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                    <  1:2864165/93‑1 (MQ=255)
                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACaaa                                                    <  2:223018/93‑1 (MQ=255)
                                     gggAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTGACGACaaa                                                    <  1:6561439/93‑1 (MQ=255)
                                       gAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATg                                                  <  2:1279993/93‑1 (MQ=255)
                                        aaGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGa                                                 <  1:3613575/93‑1 (MQ=255)
                                        aaGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGa                                                 >  2:6615761/1‑93 (MQ=255)
                                         aGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAt                                                >  2:2830116/1‑93 (MQ=255)
                                         aGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAt                                                >  1:2083264/1‑93 (MQ=255)
                                          gATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGAtt                                               >  1:6232249/1‑93 (MQ=255)
                                             gACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTcc                                            <  1:2408333/93‑1 (MQ=255)
                                             gACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTcc                                            <  1:6383305/93‑1 (MQ=255)
                                              aCGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTccc                                           >  2:5880707/1‑93 (MQ=255)
                                               cGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTcccc                                          <  2:6338839/93‑1 (MQ=255)
                                                gCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTccccc                                         <  2:5810283/93‑1 (MQ=255)
                                                 cccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCa                                        <  1:4907829/93‑1 (MQ=255)
                                                  ccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCAt                                       >  2:1592594/1‑93 (MQ=255)
                                                  ccAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCAt                                       >  1:2134579/1‑93 (MQ=255)
                                                   cAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCAtt                                      >  2:3748183/1‑93 (MQ=255)
                                                     gtgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTc                                    >  1:2800644/1‑93 (MQ=255)
                                                     gtgtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTc                                    >  1:4046732/1‑93 (MQ=255)
                                                       gtATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCAt                                  <  2:2577575/93‑1 (MQ=255)
                                                        tATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATg                                 <  2:5468951/93‑1 (MQ=255)
                                                        tATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATg                                 <  1:5825046/93‑1 (MQ=255)
                                                         aTCGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGa                                <  1:1672718/93‑1 (MQ=255)
                                                           cGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGaaa                              <  1:2397429/93‑1 (MQ=255)
                                                           cGGCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGaaa                              <  1:6355316/93‑1 (MQ=255)
                                                             gCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATg                            >  2:2973562/1‑93 (MQ=255)
                                                             gCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATg                            <  2:807807/93‑1 (MQ=255)
                                                             gCATCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATg                            >  1:4935329/1‑93 (MQ=255)
                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         <  2:313986/93‑1 (MQ=255)
                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         >  2:678309/1‑93 (MQ=255)
                                                                tCCCGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         >  2:2387098/1‑93 (MQ=255)
                                                                acccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCg                         <  2:3657054/92‑1 (MQ=255)
                                                                 cccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCgg                        >  2:5638248/1‑93 (MQ=255)
                                                                 cccGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCgg                        <  2:5431834/93‑1 (MQ=255)
                                                                   cGTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgc                      <  1:6672164/93‑1 (MQ=255)
                                                                    gTTGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTTCCCCATTTCATGAAATGGCGGgct                     >  1:2169068/1‑93 (MQ=255)
                                                                      tGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgctgc                   <  1:2237868/93‑1 (MQ=255)
                                                                      tGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgctgc                   <  2:4048764/93‑1 (MQ=255)
                                                                      tGCGCAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGgctgc                   <  2:4006286/93‑1 (MQ=255)
                                                                         gcAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCtgctgc                <  1:6057971/93‑1 (MQ=255)
                                                                         gcAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCtgctgc                >  2:3689362/1‑93 (MQ=255)
                                                                         gcAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCtgctgc                <  2:4644695/93‑1 (MQ=255)
                                                                          cAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGTTGca               >  2:211462/1‑93 (MQ=255)
                                                                          cAGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGca               >  2:3326751/1‑93 (MQ=255)
                                                                           aGTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGcac              >  1:4825105/1‑93 (MQ=255)
                                                                            gTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGgaca             >  2:4314940/1‑93 (MQ=255)
                                                                            gTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGcaca             <  1:2250794/93‑1 (MQ=255)
                                                                            gTGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGcaca             <  1:1932530/93‑1 (MQ=255)
                                                                             tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGGTGCTGCACAt            >  1:4731130/1‑93 (MQ=255)
                                                                             tGCGAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACAt            <  1:5665364/93‑1 (MQ=255)
                                                                                gAACAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggc         >  2:721299/1‑93 (MQ=255)
                                                                                  aCAGCGAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATggcgg       <  2:2687360/93‑1 (MQ=255)
                                                                                       gAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGGTGCTGCACATGGCGGCgaga  >  1:3386619/1‑93 (MQ=255)
                                                                                       gAAGCATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCgaga  >  2:3311043/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTGACCAGCTTTGCGCATTTTGAATTTCTGCACCGTGGGAAGATGACGCCCAGTGTATCGGCATCCCGTTGCGCAGTGCGAACAGCGAAACATGCGATTCGTTGAGTACAGATTCGATGTTACGACAAATGATTTCCCCCATTTCATGAAATGGCGGGCTGCTGCACATGGCGGCGAGA  >  NC_000913/1562796‑1562975

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: