Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 573,471 A→G Y34C (TAC→TGC)  rusA → DLP12 prophage; crossover junction endodeoxyribonuclease RusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913573,4710AG100.0% 1489.9 / NA 356Y34C (TAC→TGC) rusADLP12 prophage; crossover junction endodeoxyribonuclease RusA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (195/161);  total (195/161)

Only 100 of 356 total aligned reads displayed.

GAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGC  >  NC_000913/573373‑573568
                                                                                                  |                                                                                                 
gAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGcc                                                                                                 >  2:1118859/1‑101 (MQ=255)
  aTACCTACAGCATGACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCcgc                                                                                               >  2:1118863/1‑101 (MQ=255)
  aTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCcgc                                                                                               <  1:1118122/101‑1 (MQ=255)
  aTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCcgc                                                                                               >  1:1118866/1‑101 (MQ=255)
   tACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCcgcg                                                                                              >  1:1118872/1‑101 (MQ=255)
   tACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCcgcg                                                                                              >  2:1118873/1‑101 (MQ=255)
     ccTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGAt                                                                                            <  1:1118434/101‑1 (MQ=255)
      cTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATa                                                                                           <  1:1118211/101‑1 (MQ=255)
       tACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATaa                                                                                          >  1:1118881/1‑101 (MQ=255)
       tACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATaa                                                                                          >  2:1118880/1‑101 (MQ=255)
        aCAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAAc                                                                                         >  1:1118888/1‑101 (MQ=255)
        aCAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAAc                                                                                         <  1:1118275/101‑1 (MQ=255)
        aCAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAAc                                                                                         <  1:1118697/101‑1 (MQ=255)
        aCAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAAc                                                                                         <  2:1118412/101‑1 (MQ=255)
        aCAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAAc                                                                                         >  2:1118886/1‑101 (MQ=255)
              tCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGcc                                                                                   >  1:1118895/1‑101 (MQ=255)
               cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGccc                                                                                  >  2:1118898/1‑101 (MQ=255)
                   ttccccTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATa                                                                                           <  2:1118527/85‑1 (MQ=255)
                   ttACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAAt                                                                              <  1:1117834/101‑1 (MQ=255)
                   ttACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAAt                                                                              >  2:1118902/1‑101 (MQ=255)
                     aCCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCATGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCa                                                                            >  1:1118904/1‑101 (MQ=255)
                      cccTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCAt                                                                           <  2:1118375/101‑1 (MQ=255)
                       ccTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAAt                                                                              <  1:1118451/97‑1 (MQ=255)
                         tGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaa                                                                        >  2:1118910/1‑101 (MQ=255)
                         tGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaa                                                                        >  1:1118912/1‑101 (MQ=255)
                           gCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaaaa                                                                      >  2:1118918/1‑101 (MQ=255)
                           gCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaaaa                                                                      >  1:1118917/1‑101 (MQ=255)
                            ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaaaaa                                                                     <  1:1118343/101‑1 (MQ=255)
                              tCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                   >  2:1118926/1‑101 (MQ=255)
                              tCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                   >  1:1118927/1‑101 (MQ=255)
                                 gAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCaa                                                                >  1:1118936/1‑101 (MQ=255)
                                 gAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCaa                                                                >  2:1118932/1‑101 (MQ=255)
                                  aGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAAt                                                               <  1:1118508/101‑1 (MQ=255)
                                   gCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATg                                                              >  2:1118939/1‑101 (MQ=255)
                                     aataatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                            <  2:1118170/101‑1 (MQ=255)
                                      ataatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTg                                                           >  2:1118944/1‑101 (MQ=255)
                                      ataatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTg                                                           <  2:1118273/101‑1 (MQ=255)
                                      ataatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTg                                                           >  1:1118946/1‑101 (MQ=255)
                                       taatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                                          >  2:1118950/1‑101 (MQ=255)
                                       taatCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                                          >  1:1118951/1‑101 (MQ=255)
                                          tCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGata                                                       >  1:1118958/1‑101 (MQ=255)
                                           cGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGatat                                                      >  2:1118960/1‑101 (MQ=255)
                                            gCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATc                                                     <  1:1118517/101‑1 (MQ=255)
                                             cTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCg                                                    >  2:1118964/1‑101 (MQ=255)
                                             cTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCg                                                    >  1:1118965/1‑101 (MQ=255)
                                              tATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCgg                                                   <  1:1118615/101‑1 (MQ=255)
                                                ttACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGAt                                                                                            <  2:1118481/58‑1 (MQ=255)
                                                ttACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGcc                                                 <  2:1118386/101‑1 (MQ=255)
                                                 tACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCt                                                <  1:1118480/101‑1 (MQ=255)
                                                  aCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACg                                                                                        >  2:1118973/1‑60 (MQ=255)
                                                  aCCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTg                                               >  1:1118972/1‑101 (MQ=255)
                                                   ccGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGa                                                                                <  2:1118356/67‑1 (MQ=255)
                                                   ccGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAAt                                                                              <  2:1118329/69‑1 (MQ=255)
                                                   ccGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCAt                                                                           <  1:1118308/72‑1 (MQ=255)
                                                   ccGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                  <  2:1118220/97‑1 (MQ=255)
                                                   ccGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCt                                                <  1:1118491/99‑1 (MQ=255)
                                                     gCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                            >  2:1118979/1‑101 (MQ=255)
                                                     gCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCt                                            >  1:1118977/1‑101 (MQ=255)
                                                       cATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTc                                                                                      >  1:1118985/1‑57 (MQ=255)
                                                       cATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                            <  2:1118399/83‑1 (MQ=255)
                                                       cATAATCGCGGGCGCACGCACGTCAGCCCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTAt                                          <  1:1118266/101‑1 (MQ=255)
                                                        aTAATCGCGGGCGCAGGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATg                                         >  2:1118988/1‑101 (MQ=255)
                                                         tAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATc                                                     <  2:1118580/88‑1 (MQ=255)
                                                         tAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGc                                        <  1:1118417/101‑1 (MQ=255)
                                                          aaTCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGcc                                       >  1:1118993/1‑101 (MQ=255)
                                                           aTCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCt                                      <  1:1118251/101‑1 (MQ=255)
                                                           aTCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCt                                      <  2:1118253/101‑1 (MQ=255)
                                                            tCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtg                                     >  1:1118996/1‑101 (MQ=255)
                                                             cgcgGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgt                                    <  2:1118518/101‑1 (MQ=255)
                                                             cgcgGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgt                                    >  2:1118999/1‑101 (MQ=255)
                                                             cgcgGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgt                                    <  1:1118318/101‑1 (MQ=255)
                                                              gcgGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgtg                                   >  1:1119005/1‑101 (MQ=255)
                                                              gcgGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCtgtg                                   >  2:1119003/1‑101 (MQ=255)
                                                                gggCGCACGCACGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaaaa                                                                      >  1:1119010/1‑64 (MQ=255)
                                                                    gcacgcacGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATc                             >  1:1119015/1‑101 (MQ=255)
                                                                    gcacgcacGTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATc                             <  1:1118372/101‑1 (MQ=255)
                                                                     cacgcacgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATcc                            >  1:1119020/1‑101 (MQ=255)
                                                                     cacgcacgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATcc                            >  2:1119018/1‑101 (MQ=255)
                                                                       cgcacgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTaaaa                                                                      <  1:1118637/57‑1 (MQ=255)
                                                                        gcacgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCAGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCa                         >  1:1119025/1‑101 (MQ=255)
                                                                          acgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCAtt                       <  1:1118400/101‑1 (MQ=255)
                                                                          acgTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCAtt                       >  2:1119028/1‑101 (MQ=255)
                                                                            gTCAGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATcc                            <  2:1118291/94‑1 (MQ=255)
                                                                               aGCGCAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTg                  >  2:1119033/1‑101 (MQ=255)
                                                                                  gcAGAGGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGc                                             <  2:1118093/71‑1 (MQ=255)
                                                                                    agagGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcgca             <  1:1118814/101‑1 (MQ=255)
                                                                                    agagGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcaca             >  2:1119038/1‑101 (MQ=255)
                                                                                    agagGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcaca             >  1:1119036/1‑101 (MQ=255)
                                                                                     gagGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaa                                 <  2:1118446/80‑1 (MQ=255)
                                                                                     gagGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACAt            <  1:1118426/101‑1 (MQ=255)
                                                                                      agGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATg           >  1:1119044/1‑101 (MQ=255)
                                                                                      agGGGCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATg           >  1:1119040/1‑101 (MQ=255)
                                                                                       ggggCAGGCATGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGc          <  1:1118385/101‑1 (MQ=255)
                                                                                         ggcaggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCg        <  2:1118573/101‑1 (MQ=255)
                                                                                         ggcaggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCg        <  1:1118327/101‑1 (MQ=255)
                                                                                          gcaggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCgg       >  2:1119050/1‑101 (MQ=255)
                                                                                           caggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGa      >  1:1119051/1‑101 (MQ=255)
                                                                                             ggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtc    <  2:1118437/101‑1 (MQ=255)
                                                                                             ggcaTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtc    <  1:1118569/101‑1 (MQ=255)
                                                                                               caTGCCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtcgc  >  1:1119058/1‑101 (MQ=255)
                                                                                                  |                                                                                                 
GAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCGC  >  NC_000913/573373‑573568

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑