Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,770,319 C→G 5.8% P206R (CCG→CGG)  yfjR → CP4‑57 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,770,3190CG 5.8% 494.2 / 17.9 137P206R (CCG→CGG) yfjRCP4‑57 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (65/64);  new base G (4/4);  total (69/68)
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.77e-01

Only 100 of 143 total aligned reads displayed.

TACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAT  >  NC_000913/2770175‑2770464
                                                                                                                                                |                                                                                                                                                 
tACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGt                                                                                                                                                >  2:443649/1‑148 (MQ=255)
 aCTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTc                                                                                                                                               <  3:793898/148‑1 (MQ=255)
  cTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCt                                                                                                                                              <  2:592149/148‑1 (MQ=255)
    cATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttc                                                                                                                                            >  2:863318/1‑148 (MQ=255)
        caccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttcattc                                                                                                                                        <  1:443649/148‑1 (MQ=255)
          ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACTCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTc                                                                                                                                               >  2:985110/1‑139 (MQ=255)
          ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACTCTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTc                                                                                                                                               <  1:985110/139‑1 (MQ=255)
          ccaGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAg                                                                                                                                      >  3:313405/1‑148 (MQ=255)
              cagACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGTAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCg                                                                                                                                                  <  2:963774/132‑1 (MQ=255)
                    tGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                                                                             >  1:1020671/1‑131 (MQ=255)
                       tATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCt                                                                                                                                              <  1:1064688/127‑1 (MQ=255)
                       tATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCt                                                                                                                                              >  2:1064688/1‑127 (MQ=255)
                          cTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttc                                                                                                                                            <  3:953167/126‑1 (MQ=255)
                          cTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttc                                                                                                                                            >  4:953167/1‑126 (MQ=255)
                            ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGt                                                                                                                                                >  3:457741/1‑120 (MQ=255)
                            ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGt                                                                                                                                                <  4:457741/120‑1 (MQ=255)
                            ggtggCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                                                                    >  3:1133858/1‑148 (MQ=255)
                                   gAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                                   >  1:559831/1‑126 (MQ=255)
                                   gAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGcc                                                                                                                                   <  2:559831/126‑1 (MQ=255)
                                   gAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                             >  3:206590/1‑148 (MQ=255)
                                    aaCACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCa                                                                                                            <  4:313405/148‑1 (MQ=255)
                                        cGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCttt                                                                                                                 <  3:166057/139‑1 (MQ=255)
                                        cGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCttt                                                                                                                 >  4:166057/1‑139 (MQ=255)
                                          aaGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTCTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCt                                                                                                                   >  4:779103/1‑135 (MQ=255)
                                            gggCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCAc                                                                                                    <  4:206590/148‑1 (MQ=255)
                                                 tatCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCtt                                                                                               >  3:359797/1‑148 (MQ=255)
                                                  atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTa                                                                                                                                <  4:97117/114‑1 (MQ=255)
                                                  atCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATGCAGGCCTTa                                                                                                                                >  3:97117/1‑114 (MQ=255)
                                                   tCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGc                                                                                                          >  4:399504/1‑135 (MQ=255)
                                                   tCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGc                                                                                                          <  3:399504/135‑1 (MQ=255)
                                                           ttCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                                                                     <  4:1133858/148‑1 (MQ=255)
                                                            tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                                                                <  4:434138/120‑1 (MQ=255)
                                                            tCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTa                                                                                                                >  3:434138/1‑120 (MQ=255)
                                                               cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCa                                                                                                     >  2:251074/1‑128 (MQ=255)
                                                               cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCa                                                                                                     <  1:251074/128‑1 (MQ=255)
                                                               cacTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                                                                     <  1:279341/144‑1 (MQ=255)
                                                                 cTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                                                                             <  4:67720/86‑1 (MQ=255)
                                                                 cTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCtt                                                                                                                                             >  3:67720/1‑86 (MQ=255)
                                                                     atgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTg                                                                                                  <  4:1025138/125‑1 (MQ=255)
                                                                      tgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                             <  4:408203/113‑1 (MQ=255)
                                                                      tgaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                             >  3:408203/1‑113 (MQ=255)
                                                                       gaAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGcc                                                                         <  3:166650/148‑1 (MQ=255)
                                                                         aaTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcccc                                                                                                                            >  1:1071996/1‑95 (MQ=255)
                                                                         aaTCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTAcccc                                                                                                                            <  2:1071996/95‑1 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTa                                                                                                                                <  3:550297/89‑1 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTa                                                                                                                                >  4:550297/1‑89 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCg                                                                                                                     >  2:1064975/1‑100 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                                                                     <  4:535040/132‑1 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTc                                                                                     >  3:535040/1‑132 (MQ=255)
                                                                           tCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCAcc                                                                     >  3:547589/1‑148 (MQ=255)
                                                                             cATCTGATTCGGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCgg                                                                   >  1:593427/1‑148 (MQ=255)
                                                                              aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                                                                    <  3:500536/98‑1 (MQ=255)
                                                                              aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGc                                                                                                                    >  4:500536/1‑98 (MQ=255)
                                                                              aTCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGa                                                                  <  3:256655/148‑1 (MQ=255)
                                                                               tCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGt                                                                                       >  1:780358/1‑126 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttca                                                                                                                                           <  1:715390/73‑1 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCttca                                                                                                                                           >  2:715390/1‑73 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACg                                                                                            <  3:912590/120‑1 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                                >  3:1046330/1‑148 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                                >  1:575818/1‑148 (MQ=255)
                                                                                cTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                                >  1:817117/1‑148 (MQ=255)
                                                                                    ttCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgcg                                                            <  4:1046330/148‑1 (MQ=255)
                                                                                            ttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCt                                                                                                                                  >  3:628970/1‑70 (MQ=255)
                                                                                            ttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCt                                                                                                                                  <  4:628970/70‑1 (MQ=255)
                                                                                            ttGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAgcgc                                                           <  4:18021/141‑1 (MQ=255)
                                                                                               cAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgtt                                                 <  4:274024/148‑1 (MQ=255)
                                                                                                 ggAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTAt                                               <  1:46309/148‑1 (MQ=255)
                                                                                                   agaCTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACg                                                                                            <  4:538621/101‑1 (MQ=255)
                                                                                                        tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGttt                                                                                         <  2:229348/99‑1 (MQ=255)
                                                                                                        tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGc                                                                          <  4:720312/114‑1 (MQ=255)
                                                                                                        tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGc                                                                          >  3:720312/1‑114 (MQ=255)
                                                                                                        tttCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAAt                                                                <  1:935142/124‑1 (MQ=255)
                                                                                                                cAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATcc                                                                                                             >  1:649115/1‑71 (MQ=255)
                                                                                                                 aaTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgg                               <  3:824559/148‑1 (MQ=255)
                                                                                                                  aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgga                              <  2:593427/148‑1 (MQ=255)
                                                                                                                    gACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGttt                                                                                         <  4:876993/87‑1 (MQ=255)
                                                                                                                     aCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACACCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACAGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACAt                                 >  2:474116/1‑142 (MQ=255)
                                                                                                                      cAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGgaggag                          <  2:817117/148‑1 (MQ=255)
                                                                                                                       aGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGc                         <  2:575818/148‑1 (MQ=255)
                                                                                                                          cTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgga                              >  2:668710/1‑140 (MQ=255)
                                                                                                                                cGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgtt                                                 <  3:647597/115‑1 (MQ=255)
                                                                                                                                cGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgtt                                                 >  4:647597/1‑115 (MQ=255)
                                                                                                                                   cTGGTTGAAGACGCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACg                                                                                            >  4:1062594/1‑69 (MQ=255)
                                                                                                                                   cTGGTTGAAGACGCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACg                                                                                            <  3:1062594/69‑1 (MQ=255)
                                                                                                                                   cTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt                                                  >  3:599698/1‑111 (MQ=255)
                                                                                                                                   cTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt                                                  <  4:599698/111‑1 (MQ=255)
                                                                                                                                      aatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGa                                                                  <  4:85340/90‑1 (MQ=255)
                                                                                                                                      aatGAAGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGa                                                                  >  3:85340/3‑92 (MQ=255)
                                                                                                                                        tGAAGACCGGGTCTTCAACCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATgga                              <  1:111600/126‑1 (MQ=255)
                                                                                                                                          aaGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg               <  4:306722/139‑1 (MQ=255)
                                                                                                                                          aaGACCGGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg               >  3:306722/1‑139 (MQ=255)
                                                                                                                                           aGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCAt     >  2:438575/1‑148 (MQ=255)
                                                                                                                                            gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt                                                  >  3:930538/1‑102 (MQ=255)
                                                                                                                                            gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCgttgt                                                  <  4:930538/102‑1 (MQ=255)
                                                                                                                                              cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt  <  1:438575/148‑1 (MQ=255)
                                                                                                                                              cccGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt  >  4:916657/1‑148 (MQ=255)
                                                                                                                                                cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTa                                                >  2:101376/1‑100 (MQ=255)
                                                                                                                                                cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTa                                                <  1:101376/100‑1 (MQ=255)
                                                                                                                                                cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGtt             <  2:876492/135‑1 (MQ=255)
                                                                                                                                                cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGtt             >  1:876492/1‑135 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
TACTCATCCACCAGCAGACCTGGTATCTGGTGGCTGAACACGAAGGGCATATCGCCGTATTCACACTTGATGAAATCCATCTGATTCAGCCTTTGCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAT  >  NC_000913/2770175‑2770464

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: