Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,564,525 G→A 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,5250GA94.3% 272.4 / 5.2 88M1M (ATG→ATA) †yffSCPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (42/41);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 3.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620
                                                                                                 |                                                                                               
gtgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAg                                                                                               <  2:1260198/100‑1 (MQ=255)
 tgATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGt                                                                                              >  1:1516376/1‑100 (MQ=255)
   aTCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGCTTCTATGAGATGGGAGATAAAAATCATAAGTTa                                                                                            >  2:1473835/1‑100 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         >  1:392254/1‑100 (MQ=255)
      tGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGa                                                                                         >  1:1400526/1‑100 (MQ=255)
        cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                       >  1:764929/1‑100 (MQ=255)
        cGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGaaa                                                                                       >  1:1188387/1‑100 (MQ=255)
         gTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAAt                                                                                      <  2:96307/100‑1 (MQ=255)
                  cccTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTg                                                                             >  2:1018197/1‑100 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           >  2:1398957/1‑100 (MQ=255)
                    cTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCg                                                                           <  2:389771/100‑1 (MQ=255)
                        cGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACAt                                                                       >  2:694709/1‑100 (MQ=255)
                         gTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATc                                                                      >  1:863461/1‑100 (MQ=255)
                         gTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATc                                                                      >  1:794412/1‑100 (MQ=255)
                            aCGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTc                                                                   <  2:534268/100‑1 (MQ=255)
                             ngCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGCGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  >  2:636393/2‑100 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  2:691932/100‑1 (MQ=255)
                             cGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCa                                                                  <  2:152309/100‑1 (MQ=255)
                              gCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaa                                                                 <  1:1398967/100‑1 (MQ=255)
                               cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                <  1:544358/100‑1 (MQ=255)
                               cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                <  2:1307083/100‑1 (MQ=255)
                               cTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCaaa                                                                >  1:498873/1‑100 (MQ=255)
                                ttACCCTCCATCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATACGTTACGAAATCAAAATTTGCGACATCCTCAAAg                                                               <  2:1035283/100‑1 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGACATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  2:1374236/1‑100 (MQ=255)
                                   cccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgcg                                                            >  1:109852/1‑100 (MQ=255)
                                    ccTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGgtgc                                                           >  2:735453/1‑100 (MQ=255)
                                     cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          <  2:392253/100‑1 (MQ=255)
                                     cTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          >  2:570502/1‑100 (MQ=255)
                                     cTCCACCAATACGCGTGTTAACCTCCCAATGGATGCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCt                                                          >  1:1326720/1‑100 (MQ=255)
                                      tCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTg                                                         <  1:639739/100‑1 (MQ=255)
                                      tCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTg                                                         <  2:794407/100‑1 (MQ=255)
                                       ccaccaATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGc                                                        >  2:14470/1‑100 (MQ=255)
                                            aaTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATgg                                                   <  2:1536601/100‑1 (MQ=255)
                                             aTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa                                                  <  2:326805/100‑1 (MQ=255)
                                             aTACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGa                                                  <  1:768396/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                <  2:1403199/100‑1 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                >  2:990896/1‑100 (MQ=255)
                                               aCGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAg                                                <  1:1003569/100‑1 (MQ=255)
                                                 gcgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGa                                              <  2:481851/100‑1 (MQ=255)
                                                  cgTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGAc                                             <  1:990904/100‑1 (MQ=255)
                                                   gtgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACa                                            <  2:363987/100‑1 (MQ=255)
                                                    tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  2:814575/1‑100 (MQ=255)
                                                    tgtTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACaa                                           >  1:282399/1‑100 (MQ=255)
                                                      ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         >  2:276940/1‑100 (MQ=255)
                                                      ttAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAata                                         <  1:1441733/100‑1 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       >  2:43521/1‑100 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAtata                                       <  2:830085/100‑1 (MQ=255)
                                                        aaCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATAAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAAcata                                       >  2:396381/1‑100 (MQ=255)
                                                         aCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATaa                                      >  1:489936/1‑100 (MQ=255)
                                                         aCCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATaa                                      >  1:947860/1‑100 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     <  2:1448132/100‑1 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     >  2:1466153/1‑100 (MQ=255)
                                                          ccTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAg                                     <  1:1018204/100‑1 (MQ=255)
                                                            tCCCAATTGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCTACAGACGCAAAGGCGCTGCGAGGGAAGGACAATATAAggg                                   <  2:833160/100‑1 (MQ=255)
                                                             cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  >  2:5327/1‑100 (MQ=255)
                                                             cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  <  1:955851/100‑1 (MQ=255)
                                                             cccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgc                                  >  2:1259074/1‑100 (MQ=255)
                                                              ccAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcg                                 >  1:758496/1‑100 (MQ=255)
                                                               cAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgccc                                >  1:1101204/1‑98 (MQ=255)
                                                               cAATGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGgcgc                                >  1:605259/1‑100 (MQ=255)
                                                                 aTGGATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTc                              <  1:496356/100‑1 (MQ=255)
                                                                   ggATTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCaa                            <  2:271597/100‑1 (MQ=255)
                                                                     aTTCTATGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACg                          >  2:112110/1‑100 (MQ=255)
                                                                           tgggaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  2:334808/97‑1 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  1:145361/100‑1 (MQ=255)
                                                                           tGAGATGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGc                    <  1:67135/100‑1 (MQ=255)
                                                                            gagaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGcc                   >  2:1521967/1‑100 (MQ=255)
                                                                            gagaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGcc                   >  1:1549341/1‑100 (MQ=255)
                                                                             agaTGGGAGATAAAACTCATAAGTTACGATATCAAAATTTGCGACAGCCTCAAAGGCGCTGCGATGGAAGTACAATATAAGGGCGCTCAACGTGGGGCCa                  <  2:542883/100‑1 (MQ=255)
                                                                             agaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCa                  >  2:1455264/1‑100 (MQ=255)
                                                                              gaTGGGAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCaa                 >  2:231798/1‑100 (MQ=255)
                                                                                 gggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATg              >  2:557931/1‑100 (MQ=255)
                                                                                  ggAGATAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGc             >  1:877096/1‑100 (MQ=255)
                                                                                   gagaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCg            <  2:1516366/100‑1 (MQ=255)
                                                                                    agaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCga           <  2:605256/100‑1 (MQ=255)
                                                                                    agaTAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCga           <  1:1521977/100‑1 (MQ=255)
                                                                                     cttAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCGCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          <  1:687356/98‑1 (MQ=255)
                                                                                      aTAAAAATCATAAGTCACGAAATCAAAATCTGCGACATCCCCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagg         <  2:1400516/100‑1 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  1:373492/1‑100 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgagga        >  2:780958/1‑100 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          >  2:687353/1‑98 (MQ=255)
                                                                                       tAAAAATCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCgag          >  2:1401620/1‑98 (MQ=255)
                                                                                        aaaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGTGCTCAACGTGGGGCCAAATGCgaggag       <  1:1242729/100‑1 (MQ=255)
                                                                                         aaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGgaga      >  1:61748/1‑100 (MQ=255)
                                                                                         aaaaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGgaga      >  2:202412/1‑100 (MQ=255)
                                                                                           aaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATc    <  1:117623/100‑1 (MQ=255)
                                                                                           aaTCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATc    <  1:711019/100‑1 (MQ=255)
                                                                                             tCATAAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGc  <  2:1188380/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTGATCTGCGTCAATATTCCCTTCCGTTACGCTTACCCTCCACCAATACGCGTGTTAACCTCCCAATGGATTCTATGAGATGGGAGATAAAAATCATGAGTTACGAAATCAAAATTTGCGACATCCTCAAAGGCGCTGCGATGGAAGGACAATATAAGGGCGCTCAACGTGGGGCCAAATGCGAGGAGATCGC  >  NC_000913/2564428‑2564620

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