Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,557,200 T→C V137A (GTG→GCG)  b2442 → putative prophage integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,557,2000TC100.0% 801.8 / NA 251V137A (GTG→GCG) b2442putative prophage integrase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (114/137);  total (114/137)

Only 100 of 274 total aligned reads displayed.

CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289
                                                                                          |                                                                                         
cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         >  1:7040419/1‑93 (MQ=255)
cccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCgg                                                                                         <  2:6721529/93‑1 (MQ=255)
 ccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAggcggc                                                                                        >  1:8309795/1‑93 (MQ=255)
 ccTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAggcggc                                                                                        >  1:4377895/1‑93 (MQ=255)
 ccTTATCGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAggcggc                                                                                        <  2:2142382/93‑1 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       <  2:6110128/93‑1 (MQ=255)
  cTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGgcggcg                                                                                       >  1:8214696/1‑93 (MQ=255)
     aTTGCCAATGGAATTAACCCAGTGGAACAAGCTCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGAATTTAAGAAGGCGGCGaaa                                                                                    <  1:8444603/93‑1 (MQ=255)
     aTTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGaa                                                                                    <  2:2027029/93‑1 (MQ=255)
       tGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGa                                                                                  <  1:1328567/93‑1 (MQ=255)
        gCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGAc                                                                                 <  1:9486394/93‑1 (MQ=255)
         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                <  1:1949882/93‑1 (MQ=255)
         ccAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACt                                                                                <  1:9242818/93‑1 (MQ=255)
          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               >  2:7390035/1‑93 (MQ=255)
          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               >  2:4557578/1‑93 (MQ=255)
          cAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTg                                                                               <  2:7858996/93‑1 (MQ=255)
           aaTGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTgg                                                                              <  2:1444901/93‑1 (MQ=255)
             tGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGtt                                                                            >  1:4337161/1‑93 (MQ=255)
               gAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGGTTg                                                                          >  2:8910688/1‑93 (MQ=255)
                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         <  2:3844338/93‑1 (MQ=255)
                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         <  1:2915639/93‑1 (MQ=255)
                aaTTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGc                                                                         >  1:8412488/1‑93 (MQ=255)
                  ttAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTa                                                                       <  2:3290495/93‑1 (MQ=255)
                  ttAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTa                                                                       >  2:5357332/1‑93 (MQ=255)
                     aCCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCa                                                                    <  1:10023769/93‑1 (MQ=255)
                      cccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAg                                                                   <  1:4970396/93‑1 (MQ=255)
                       ccAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAgg                                                                  <  2:9770912/93‑1 (MQ=255)
                         aGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGt                                                                >  1:2391803/1‑93 (MQ=255)
                          gTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGtt                                                               >  1:7742485/1‑93 (MQ=255)
                           tGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTg                                                              >  2:56010/1‑93 (MQ=255)
                            ggAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgg                                                             <  2:5614008/93‑1 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTggg                                                            >  2:8540975/1‑92 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            <  2:9501304/93‑1 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            <  2:8594115/93‑1 (MQ=255)
                             gAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGa                                                            <  2:1709823/93‑1 (MQ=255)
                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAg                                                           <  1:1326175/93‑1 (MQ=255)
                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAg                                                           <  2:195063/93‑1 (MQ=255)
                              aaCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAg                                                           <  1:7973379/93‑1 (MQ=255)
                               aCAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAgg                                                          >  1:2784404/1‑93 (MQ=255)
                                cAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTgggggg                                                         >  1:1161192/1‑93 (MQ=255)
                                cAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAggg                                                         <  2:8773316/93‑1 (MQ=255)
                                cAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAggg                                                         <  1:5165778/93‑1 (MQ=255)
                                cAAGCCAGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCTGCAGGATTTAAGAAGGCGGCGGAAGACTCGTTTGCTACCAGGGTTGGAggg                                                         <  2:3648074/93‑1 (MQ=255)
                                 aaGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAgggg                                                        >  2:7349431/1‑92 (MQ=255)
                                  aGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGGGggtg                                                       >  1:2303892/1‑93 (MQ=255)
                                  aGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAggggg                                                       >  1:3135364/1‑91 (MQ=255)
                                     ccGCAAGGAAAAAGCCATTGATGCGCACAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGGTGGAGGGTGGTc                                                    >  2:75869/1‑93 (MQ=255)
                                     ccGCAAGGAAAAAGACATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTc                                                    <  1:2120887/93‑1 (MQ=255)
                                      cGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCg                                                   <  2:881668/93‑1 (MQ=255)
                                        cAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGGGGAAGACTGGTTTGCTACCAGGGTTGGAGGGGGGTCGGa                                                 >  1:7234224/1‑93 (MQ=255)
                                        cAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGa                                                 <  1:9771202/93‑1 (MQ=255)
                                        cAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGa                                                 <  2:4375438/93‑1 (MQ=255)
                                        cAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGa                                                 >  1:3936000/1‑93 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  1:9127517/1‑93 (MQ=255)
                                          aGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAAt                                               >  1:4860170/1‑93 (MQ=255)
                                           ggAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATc                                              <  1:4529472/93‑1 (MQ=255)
                                           ggAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATc                                              >  2:6592329/1‑93 (MQ=255)
                                             aaaaaGCCATTGATGCGCTCAACGTGGCGGCAGGGTTTAAGAAGGCGGCGGAAGACTGGCTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCt                                            >  2:6933551/1‑93 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           <  2:3617289/93‑1 (MQ=255)
                                              aaaaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTa                                           >  1:641011/1‑93 (MQ=255)
                                                aaGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATg                                         <  2:8296556/93‑1 (MQ=255)
                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGc                                        >  2:6327888/1‑93 (MQ=255)
                                                 aGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGc                                        <  2:6028707/93‑1 (MQ=255)
                                                   ccATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaa                                      <  1:3308061/93‑1 (MQ=255)
                                                     aTTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCaaaa                                    >  1:5589299/1‑93 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTGGGAATCCTATGCAAAAc                                   >  1:3030091/1‑93 (MQ=255)
                                                      ttGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAAc                                   <  1:1524619/93‑1 (MQ=255)
                                                       tGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACa                                  <  2:4167805/93‑1 (MQ=255)
                                                       tGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACa                                  <  2:1599309/93‑1 (MQ=255)
                                                        gATGCGCTCAACATGGCGGCAGGATTTAAGAGGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTTGGAATCCTATGCAAAACAg                                 <  1:7462060/93‑1 (MQ=255)
                                                        gATGCGCTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAg                                 >  1:1344453/1‑93 (MQ=255)
                                                             gcTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAg                            <  1:6701648/93‑1 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAgt                           >  2:3566833/1‑92 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           <  1:8829052/93‑1 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           <  1:691203/93‑1 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           >  2:588903/1‑93 (MQ=255)
                                                              cTCAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGa                           <  1:3415932/93‑1 (MQ=255)
                                                               ccAACATGGCGGCAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGAt                          <  1:5524424/92‑1 (MQ=255)
                                                                      ggcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACt                   <  1:1893138/93‑1 (MQ=255)
                                                                       gcggcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  >  2:7122613/1‑93 (MQ=255)
                                                                       gcggcAGGATTGAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGTTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTa                  <  2:4290956/93‑1 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  2:9047973/1‑93 (MQ=255)
                                                                          gcAGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaa               >  2:1088877/1‑93 (MQ=255)
                                                                            aGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaa             >  2:7205099/1‑93 (MQ=255)
                                                                            aGGATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaa             >  1:9198312/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  2:2787815/1‑93 (MQ=255)
                                                                             ggATTTAAGAAGGCGGCGGAAGACTGGGTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGaaaaa            >  2:3193647/1‑93 (MQ=255)
                                                                              gATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAg           >  2:9142951/1‑93 (MQ=255)
                                                                              gATTTAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAg           >  1:412948/1‑93 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         >  1:1049825/1‑93 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         >  2:5812328/1‑93 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  2:4843943/93‑1 (MQ=255)
                                                                                tttAAGAAGGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGAc         <  1:9886684/93‑1 (MQ=255)
                                                                                   aagaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTc      <  2:6421465/93‑1 (MQ=255)
                                                                                   aagaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTc      <  1:6083991/93‑1 (MQ=255)
                                                                                    agaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCt     >  2:3806113/1‑93 (MQ=255)
                                                                                     gaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTa    <  1:5687807/93‑1 (MQ=255)
                                                                                     gaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTa    <  1:6572187/93‑1 (MQ=255)
                                                                                     gaagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTa    <  1:7625347/93‑1 (MQ=255)
                                                                                      aagGCGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTAt   <  1:4713435/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCCTTATTGCCAATGGAATTAACCCAGTGGAACAAGCCCGCAAGGAAAAAGCCATTGATGCGCTCAACATGGCGGCAGGATTTAAGAAGGTGGCGGAAGACTGGTTTGCTACCAGGGTTGGAGGGTGGTCGGAATCCTATGCAAAACAGGTCAGATCGGCACTAGAAAAAGACGTCTATC  >  NC_000913/2557110‑2557289

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