Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,564,525 G→A 100% M1M (ATG→ATA) † yffS → CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,564,5250GA99.5% 656.2 / ‑5.0 191M1M (ATG→ATA) †yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (97/93);  total (97/94)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 191 total aligned reads displayed.

GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620
                                                                                                 |                                                                                               
gtgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAg                                                                                               <  1:2266405/100‑1 (MQ=255)
  gATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGtt                                                                                             >  1:2355799/1‑100 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCATCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         >  1:1483945/1‑100 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         <  2:2568565/100‑1 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         <  2:2163033/100‑1 (MQ=255)
       gCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaa                                                                                        >  2:556972/1‑100 (MQ=255)
        cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                       <  2:1882316/100‑1 (MQ=255)
         gTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGCTTCTATGAGATGGGCGATAAAAATCATAAGTTACg                                                                                          >  1:1858075/1‑96 (MQ=255)
         gTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAAt                                                                                      <  2:927395/100‑1 (MQ=255)
          tCACTATTCCCTTCCGTTACGCTTACCCTCCCCCACTACGCGTGTCAACCTCCCAATGGATTCTATGAGACGGGCGATAACAATCATAAGTTACGAAATc                                                                                     >  1:2096035/1‑100 (MQ=255)
          tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  2:2234780/1‑100 (MQ=255)
          tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     >  1:4084507/1‑100 (MQ=255)
          tCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATc                                                                                     <  1:4237985/100‑1 (MQ=255)
           cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    >  2:2401901/1‑100 (MQ=255)
           cAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCa                                                                                    <  1:3378112/100‑1 (MQ=255)
            aaTATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCaa                                                                                   >  2:3377390/1‑100 (MQ=255)
               atTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAt                                                                                <  1:3674215/100‑1 (MQ=255)
                ttCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAAtt                                                                               <  2:4178526/100‑1 (MQ=255)
                   ccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGc                                                                            <  1:702635/100‑1 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:86714/100‑1 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  1:1880244/100‑1 (MQ=255)
                     ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          >  2:1929090/1‑100 (MQ=255)
                     ttCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGa                                                                          <  2:2541533/100‑1 (MQ=255)
                      tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                         <  1:53200/100‑1 (MQ=255)
                      tCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGAc                                                                         >  1:357369/1‑100 (MQ=255)
                       ccGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACa                                                                        >  2:181916/1‑100 (MQ=255)
                         gTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATc                                                                      >  2:1751506/1‑100 (MQ=255)
                           tACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCt                                                                    <  1:1186590/100‑1 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTACCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  2:3446194/100‑1 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  2:205405/100‑1 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  >  2:2984329/1‑100 (MQ=255)
                              gCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaa                                                                 >  1:3009068/1‑100 (MQ=255)
                              gCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaa                                                                 >  1:960559/1‑100 (MQ=255)
                                 tACCCTCCACCAATACGCGTGTTAACCTCCCAATGGCTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAgg                                                              >  1:670522/1‑100 (MQ=255)
                                 tACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAg                                                               <  1:3345081/99‑1 (MQ=255)
                                  aCCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgc                                                             >  2:2520101/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  2:1488906/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCATCaaa                                                                >  1:1455449/1‑96 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgn                                                           <  2:2102682/100‑2 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcgc                                                           >  1:216489/1‑100 (MQ=255)
                                       ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                        >  1:261777/1‑100 (MQ=255)
                                       ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                        <  2:576379/100‑1 (MQ=255)
                                          ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     >  1:3950173/1‑100 (MQ=255)
                                          ccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGAt                                                     <  2:3286643/100‑1 (MQ=255)
                                            aaTACGGGTGCTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCTAAGTTTGCGACATCCTCAAAGGCGCTGCGATgg                                                   <  1:2809454/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                >  2:1433026/1‑100 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                <  2:3687478/100‑1 (MQ=255)
                                                cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               <  2:3809266/100‑1 (MQ=255)
                                                cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               <  2:3986036/100‑1 (MQ=255)
                                                cgcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAgg                                               >  2:3635358/1‑100 (MQ=255)
                                                  cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  1:2555680/100‑1 (MQ=255)
                                                  cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             >  1:1396218/1‑100 (MQ=255)
                                                  cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  1:401079/100‑1 (MQ=255)
                                                   gtgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACa                                            >  2:895784/1‑100 (MQ=255)
                                                      ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         <  2:808110/100‑1 (MQ=255)
                                                       tAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAatat                                        <  2:1369787/100‑1 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       <  1:1117477/100‑1 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       >  1:1517404/1‑100 (MQ=255)
                                                         aCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATaa                                      <  2:3055290/100‑1 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     >  1:3303871/1‑100 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     <  1:510296/100‑1 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     >  2:1368106/1‑100 (MQ=255)
                                                            tCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAggg                                   <  1:3824304/100‑1 (MQ=255)
                                                             cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  >  2:435027/1‑100 (MQ=255)
                                                               cAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcgc                                >  1:1676977/1‑100 (MQ=255)
                                                                aaTGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTCGCGACATCCGCAAAGACGCTGCGACGGAAGGACAATATAAGGgcgc                                >  1:2886558/1‑99 (MQ=255)
                                                                   ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            >  1:3476749/1‑100 (MQ=255)
                                                                   ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            >  1:402530/1‑100 (MQ=255)
                                                                   ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            >  2:2653760/1‑100 (MQ=255)
                                                                    gATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAAc                           >  2:3115140/1‑100 (MQ=255)
                                                                     aTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACg                          >  1:2185733/1‑100 (MQ=255)
                                                                     aTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACg                          <  1:232403/100‑1 (MQ=255)
                                                                     aTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACg                          >  1:473913/1‑100 (MQ=255)
                                                                     aTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACg                          >  1:3886740/1‑100 (MQ=255)
                                                                      ttCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCATCAAAGGCGCTGCGATGGAAGGACAATATAAGGACGCTCAACGt                         >  1:1699241/1‑100 (MQ=255)
                                                                       tCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTg                        <  2:2729498/100‑1 (MQ=255)
                                                                       tCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTg                        >  2:1201644/1‑100 (MQ=255)
                                                                       tCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTg                        <  2:342925/100‑1 (MQ=255)
                                                                        cTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgg                       >  2:707290/1‑100 (MQ=255)
                                                                          aTGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTgggg                     >  1:1331466/1‑100 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  2:3794205/100‑1 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  1:2520094/100‑1 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  2:2125821/100‑1 (MQ=255)
                                                                               aTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaa                 >  1:473448/1‑99 (MQ=255)
                                                                                  ggAGATAAAAAGCATAAGTTACGAAATCGAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGc             <  1:854363/100‑1 (MQ=255)
                                                                                     gaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  1:3662918/100‑1 (MQ=255)
                                                                                      aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  2:3353942/100‑1 (MQ=255)
                                                                                      aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  1:3852587/100‑1 (MQ=255)
                                                                                      aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  1:1021247/100‑1 (MQ=255)
                                                                                      aTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATAATCAAAGGCGCTGCGATGGAAGGACAATATAAGTGCGATCAACGTGGGGCCAAATGc             >  1:2668525/1‑96 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGGCAAATGCgag          >  1:3529078/1‑98 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:3774278/1‑100 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        <  1:982232/100‑1 (MQ=255)
                                                                                         aaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGgaga      <  2:1321762/100‑1 (MQ=255)
                                                                                          aaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGAt     <  1:1360575/100‑1 (MQ=255)
                                                                                           aaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATc    <  2:3144929/100‑1 (MQ=255)
                                                                                            aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   <  1:2071523/100‑1 (MQ=255)
                                                                                            aTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   >  2:3206433/1‑100 (MQ=255)
                                                                                             tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  2:1581068/100‑1 (MQ=255)
                                                                                              cATAAGTTACGAAATCAAAATTTGCGACATGCTCAAAGGCGCTGCGATGGAAGTACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCg   <  1:4168701/98‑1 (MQ=255)
                                                                                                 |                                                                                               
GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620

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