Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,764,382 C→A 6.5% D258Y (GAC→TAC)  yfjL ← CP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,764,3820CA 6.5% 264.6 / 2.5 78D258Y (GAC→TAC) yfjLCP4‑57 putative defective prophage, DUF4297/DUF1837 polymorphic toxin family protein
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base A (2/4);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTC  >  NC_000913/2764259‑2764521
                                                                                                                           |                                                                                                                                           
tATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGAtt                                                                                                                             >  2:190651/1‑140 (MQ=255)
tATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGAtt                                                                                                                             >  2:1539012/1‑140 (MQ=255)
 aTCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGAttt                                                                                                                            >  2:1351346/1‑140 (MQ=255)
  tCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGAtttt                                                                                                                           >  1:404731/1‑140 (MQ=255)
        aTGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGcc                                                                                                                     >  2:1055652/1‑140 (MQ=255)
        aTGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTACCGTTTGACGTAGATTTTGCGGcc                                                                                                                     >  2:1078353/1‑140 (MQ=255)
        aTGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTACCGTTTGACGTAGATTTTGCGGcc                                                                                                                     >  2:1534660/1‑140 (MQ=255)
          gtgCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAAGGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTg                                                                                                                   <  1:1083630/140‑1 (MQ=255)
            gCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGAc                                                                                                                 <  2:907665/140‑1 (MQ=255)
               tGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCga                                                                                                              <  1:1359760/140‑1 (MQ=255)
               tGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCga                                                                                                              <  1:721333/140‑1 (MQ=255)
                     ggTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTgg                                                                                                        >  1:1491373/1‑140 (MQ=255)
                      gTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGt                                                                                                       <  2:711812/140‑1 (MQ=255)
                       gctttCTTAAGCGAGGTCTGCAGATGTCGACCCTGGCCTCTAACCTTCGCGCAGACCGAAGCCAAAGGCGCTCTACAAGGCGTCCGGCTGAGGGGTAAGTCCCGAGTGACGTAGCTTTTGCGGCCTGACCGAGATTGGTc                                                                                                      <  2:720369/138‑1 (MQ=255)
                        aTTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCa                                                                                                     >  2:1453840/1‑140 (MQ=255)
                         tttCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCAt                                                                                                    <  2:1526740/140‑1 (MQ=255)
                                 gTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGaca                                                                                            >  2:1572715/1‑140 (MQ=255)
                                 gTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGaca                                                                                            >  2:200530/1‑140 (MQ=255)
                                   gAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTATATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACg                                                                                          <  1:83466/140‑1 (MQ=255)
                                   gAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACg                                                                                          <  2:769491/140‑1 (MQ=255)
                                   gAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACg                                                                                          <  2:918886/140‑1 (MQ=255)
                                   gAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACg                                                                                          <  1:1460338/140‑1 (MQ=255)
                                    aGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACgc                                                                                         >  1:353622/1‑140 (MQ=255)
                                      tACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACgctg                                                                                       >  1:360989/1‑140 (MQ=255)
                                      tACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACgctg                                                                                       >  1:360939/1‑140 (MQ=255)
                                              aTGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAAcc                                                                               >  1:1233023/1‑140 (MQ=255)
                                                gTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAACTGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCAc                                                                             <  1:1430247/140‑1 (MQ=255)
                                                gTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCAc                                                                             <  2:404731/140‑1 (MQ=255)
                                                gTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCAc                                                                             <  2:1357625/140‑1 (MQ=255)
                                                   gAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCa                                                                          >  1:400826/1‑140 (MQ=255)
                                                   gAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCa                                                                          >  1:403535/1‑140 (MQ=255)
                                                     aCCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAg                                                                        >  2:235176/1‑140 (MQ=255)
                                                      ccTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGt                                                                       <  1:1678881/140‑1 (MQ=255)
                                                      ccTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGt                                                                       <  1:1218678/140‑1 (MQ=255)
                                                             tttACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTg                                                                <  1:503602/140‑1 (MQ=255)
                                                               tACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAATTGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGt                                                              <  2:1367414/140‑1 (MQ=255)
                                                               tACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGt                                                              >  2:565649/1‑140 (MQ=255)
                                                                aCGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAACTGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGAGTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTc                                                             <  1:717409/140‑1 (MQ=255)
                                                                aCGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTc                                                             <  1:1453840/140‑1 (MQ=255)
                                                                aCGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTc                                                             <  1:1021423/140‑1 (MQ=255)
                                                                aCGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTc                                                             <  1:1500124/140‑1 (MQ=255)
                                                                      aGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCACGTCGGTCTGGTCGTAAGt                                                       >  2:1527967/1‑140 (MQ=255)
                                                                      aGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGt                                                       <  2:299680/140‑1 (MQ=255)
                                                                        cgcAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTACCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAg                                                     <  1:1534660/140‑1 (MQ=255)
                                                                        cgcAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTACCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAg                                                     <  1:1078353/140‑1 (MQ=255)
                                                                           aGACCGAAGCCATAGGCGCTCTATATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGAc                                                  <  1:904860/140‑1 (MQ=255)
                                                                           aGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGAc                                                  <  1:664990/140‑1 (MQ=255)
                                                                             aCCGAAGCCATAGTCGCTCTATATGGCGTTAGGATGTTGGGTATGTCGCGTTTGAAGTTGATTTTGCGGCCTGACAGTGATTGTTCATCTTCGTCTCGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGTCtt                                                <  2:717847/140‑1 (MQ=255)
                                                                             aCCGAAGCCATAGGCGCTCTACATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGAACGAGATTGGTCATCTTCGACACGCTGCTTGAACCCCTCAAGTAGGTCTGGTCGTAAGTAGGACtt                                                <  1:762985/140‑1 (MQ=255)
                                                                             aCCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACtt                                                <  2:785471/140‑1 (MQ=255)
                                                                             aCCGAAGCCATAGGAGCTCTAACAGGCGTAAGGATGTGCGGTAAGTACCGATTGACGTAGATTTTGCGGCCTGACCGAGATTGGACATCTTCGACACGCTGCTTGAACCACTCCAGTAGGTCTGGTCGTAAGTAGGACtt                                                <  2:1676188/140‑1 (MQ=255)
                                                                               cGAAGACATAGGCGCTCTACATGGCTTAAGGATGTGGGGTAAGTACCGTTTGTAGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGt                                              <  1:909637/140‑1 (MQ=255)
                                                                                 aaGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGtca                                            >  1:1318723/1‑140 (MQ=255)
                                                                                  aGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGtcat                                           <  1:216943/140‑1 (MQ=255)
                                                                                  aGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGtcat                                           <  1:1156783/140‑1 (MQ=255)
                                                                                  aGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGtcat                                           <  1:218845/140‑1 (MQ=255)
                                                                                  aGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGATTGACGTAGATTTTGCGGCCTGACCGAGATAGGTCATCATCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGACGTAAGTAGGACTTGtcat                                           >  1:52335/1‑140 (MQ=255)
                                                                                       tAGGCGCTCTCCATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCCCGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCa                                      <  2:825405/140‑1 (MQ=255)
                                                                                         ggCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACt                                    >  2:376465/1‑140 (MQ=255)
                                                                                         ggCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACt                                    >  2:1066018/1‑140 (MQ=255)
                                                                                         ggCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACCCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACt                                    >  2:376466/1‑140 (MQ=255)
                                                                                            gcTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACtgtg                                 <  1:127793/140‑1 (MQ=255)
                                                                                            gcTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACtgtg                                 <  1:127611/140‑1 (MQ=255)
                                                                                              tctAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGaa                               >  2:302717/1‑140 (MQ=255)
                                                                                              tctAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGaa                               >  2:302721/1‑140 (MQ=255)
                                                                                              tctAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGaa                               >  1:472733/1‑140 (MQ=255)
                                                                                               ctAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCAGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAAt                              <  1:887766/140‑1 (MQ=255)
                                                                                                        gTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTAGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCgaga                     >  2:1531120/1‑140 (MQ=255)
                                                                                                            ggATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGc                 >  2:1194511/1‑140 (MQ=255)
                                                                                                                 tgGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGcc            >  1:415232/1‑140 (MQ=255)
                                                                                                                 tgGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGcc            >  2:811604/1‑140 (MQ=255)
                                                                                                                 tgGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGcc            >  2:274893/1‑140 (MQ=255)
                                                                                                                       caGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGGCATCGTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGCGAATGCTCGCCCTTgcg      <  2:1082795/139‑1 (MQ=255)
                                                                                                                       aaGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTgcg      >  2:504162/1‑140 (MQ=255)
                                                                                                                       aaGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTgcg      <  2:1087563/140‑1 (MQ=255)
                                                                                                                       aaGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATGGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTAGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCATTgcg      <  2:1082794/140‑1 (MQ=255)
                                                                                                                       aaGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGAGTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTgcg      <  2:1082750/140‑1 (MQ=255)
                                                                                                                        aGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTgcgc     >  1:707009/1‑140 (MQ=255)
                                                                                                                           cccGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTc  <  1:1017661/140‑1 (MQ=255)
                                                                                                                           |                                                                                                                                           
TATCCGCGATGTGCTTGTAGCGGTATTTCTTAAGTGAGTACTGCTGATGTAGAACCTGGCCTTTACGCTTAGCGCAGACCGAAGCCATAGGCGCTCTAAATGGCGTAAGGATGTGGGGTAAGTCCCGTTTGACGTAGATTTTGCGGCCTGACCGAGATTGGTCATCTTCGACACGCTGCTTGAACCACTCAAGTAGGTCTGGTCGTAAGTAGGACTTGTCATCAGCACTGTGAATGCGTCGAGAATGCTCGCCCTTGCGCGTC  >  NC_000913/2764259‑2764521

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