Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,431,882 T→C V945A (GTG→GCG)  stfR → Rac prophage, putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,431,8820TC98.7% 260.3 / ‑6.4 74V945A (GTG→GCG) stfRRac prophage, putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (38/35);  minor base A (1/0);  total (39/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAA  >  NC_000913/1431747‑1432016
                                                                                                                                       |                                                                                                                                      
caccaGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGAt                                                                                                                                     <  1:442580/139‑1 (MQ=255)
             cccAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGt                                                                                                                        >  1:893861/1‑139 (MQ=255)
               cAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTgg                                                                                                                      <  2:893861/139‑1 (MQ=255)
                 gTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCt                                                                                                                    >  2:1733666/1‑139 (MQ=255)
                  tCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTg                                                                                                                   >  2:1199704/1‑139 (MQ=255)
                    ggCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACACATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 <  1:1281983/139‑1 (MQ=255)
                    ggCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 >  1:1226889/1‑139 (MQ=255)
                    ggCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 <  1:621210/139‑1 (MQ=255)
                    agacagccAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGTTGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 <  1:1257406/133‑1 (MQ=255)
                     gcgcACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGAc                                                                                                                >  1:526813/1‑139 (MQ=255)
                                   tGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  >  1:1576381/1‑139 (MQ=255)
                                     gCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGAAGGGGGAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGAGATCCCTGATATGCGTGGCTGGACGATTAAGGGcctgtct                                                                                                >  2:1257406/1‑133 (MQ=14)
                                                  gTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                   >  2:1160637/1‑139 (MQ=255)
                                                  gTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCggg                                                                                   <  2:905469/139‑1 (MQ=255)
                                                   ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 <  2:36306/108‑1 (MQ=255)
                                                   ttCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGa                                                                                                                 >  1:36306/1‑108 (MQ=255)
                                                       gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                              >  1:1415658/1‑139 (MQ=255)
                                                       gTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTa                                                                              <  1:1364328/139‑1 (MQ=255)
                                                           ggTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGt                                                                          <  2:1226889/139‑1 (MQ=255)
                                                             ttATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTg                                                                                                                                <  2:283238/83‑1 (MQ=255)
                                                             ttATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTg                                                                                                                                >  1:283238/1‑83 (MQ=255)
                                                              tATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                       <  1:1160637/139‑1 (MQ=255)
                                                                  cccTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGa                                                                                                                                      >  1:1081499/1‑72 (MQ=255)
                                                                    cTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                       <  1:920171/133‑1 (MQ=255)
                                                                    cTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                       >  2:920171/1‑133 (MQ=255)
                                                                      gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                       >  2:274065/1‑131 (MQ=255)
                                                                      gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGtctc                                                                       <  1:274065/131‑1 (MQ=255)
                                                                      gATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAg                                                               >  2:1593047/1‑139 (MQ=255)
                                                                       aTGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAgg                                                              >  2:542036/1‑139 (MQ=255)
                                                                        tGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGa                                                             <  2:1576381/139‑1 (MQ=255)
                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATTCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTgg                                                                                        <  1:1332736/108‑1 (MQ=255)
                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTgg                                                                                        >  2:1332736/1‑108 (MQ=255)
                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACg                                                           >  1:1178645/1‑137 (MQ=255)
                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGc                                                         >  1:1604226/1‑139 (MQ=255)
                                                                            ggggCAGGCTTTTGACAAATCTGCTTACCCGAAAATTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACg                                                           <  2:1178645/137‑1 (MQ=255)
                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  >  2:68127/1‑97 (MQ=255)
                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  <  1:68127/97‑1 (MQ=255)
                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCa                                                        >  2:1752320/1‑139 (MQ=255)
                                                                             gggCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCa                                                        >  2:1356192/1‑139 (MQ=255)
                                                                               gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAg                                                                     <  1:1245970/124‑1 (MQ=255)
                                                                               gCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAg                                                                     >  2:1245970/1‑124 (MQ=255)
                                                                                cAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTa                                                     <  1:1356192/139‑1 (MQ=255)
                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                               <  2:1062956/112‑1 (MQ=255)
                                                                                 aGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGt                                                                               >  1:1062956/1‑112 (MQ=255)
                                                                                  ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  <  1:956874/92‑1 (MQ=255)
                                                                                  ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  >  2:956874/1‑92 (MQ=255)
                                                                                  ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTaaa                                                   >  2:980238/1‑139 (MQ=255)
                                                                                  ggCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTaaa                                                   >  1:853784/1‑139 (MQ=255)
                                                                                        tGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCAt                                             >  1:1607375/1‑139 (MQ=255)
                                                                                         gACAAATCTGCTTACCCGAAACTTGCAGTCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATa                                            <  2:1714825/139‑1 (MQ=255)
                                                                                             aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGc                                                                                  <  2:1195189/97‑1 (MQ=255)
                                                                                             aaTCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGc                                                                                  >  1:1195189/1‑97 (MQ=255)
                                                                                                 tGCTTACCCGAAACTTGCGGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagc                                    >  2:324248/1‑139 (MQ=255)
                                                                                                   cTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATa                                            >  1:476186/1‑129 (MQ=255)
                                                                                                   cTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATa                                            <  2:476186/129‑1 (MQ=255)
                                                                                                   cTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgc                                  >  2:49735/1‑139 (MQ=255)
                                                                                                       cccGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagc                              <  1:980238/139‑1 (MQ=255)
                                                                                                         cGAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagcgc                            <  2:526813/139‑1 (MQ=255)
                                                                                                          gAAACTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCa                           <  1:1336563/139‑1 (MQ=255)
                                                                                                           aaaCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGc                                                                                                      >  1:360479/1‑63 (MQ=255)
                                                                                                           aaaCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGc                                                                                                      <  2:360479/63‑1 (MQ=255)
                                                                                                            aaCTTGCGGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATc                         <  1:324248/139‑1 (MQ=255)
                                                                                                            aaCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATc                         <  2:1415658/139‑1 (MQ=255)
                                                                                                            aaCTTGCAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATc                         <  2:561790/139‑1 (MQ=255)
                                                                                                                  cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagc                              >  1:1190568/1‑128 (MQ=255)
                                                                                                                  cAGCCGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCAcagcgccagc                              <  2:1190568/128‑1 (MQ=255)
                                                                                                                      cGCTTATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGAt               >  1:1553137/1‑139 (MQ=255)
                                                                                                                         ttATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  <  1:1745044/53‑1 (MQ=16)
                                                                                                                         ttATCCGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAAc                                                                                                  >  2:1745044/1‑53 (MQ=255)
                                                                                                                              cGTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGAc       <  2:1492468/139‑1 (MQ=255)
                                                                                                                               gTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCg                                                <  1:1622456/97‑1 (MQ=255)
                                                                                                                               gTCAGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCg                                                >  2:1622456/1‑97 (MQ=255)
                                                                                                                                  aGGCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGaaa   >  2:780238/1‑139 (MQ=255)
                                                                                                                                   ggCGCGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGaaaa  >  1:311757/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                      
CACCAGAAAGCTACCCAGTCGGCGCACCAATCCCGTGGCCATCAGATACCGTTCCGTCTGGTTATGCCCTGATGCAGGGGCAGGCTTTTGACAAATCTGCTTACCCGAAACTTGCAGCCGCTTATCCGTCAGGCGTGATCCCTGATATGCGTGGCTGGACGATTAAGGGCAAACCTGCCAGTGGTCGGGCCGTATTGTCTCAGGAACAGGACGGCATTAAATCGCATACCCACAGCGCCAGCGCATCCAGTACAGATTTGGGGACGAAAA  >  NC_000913/1431747‑1432016

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