Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,356 C→G 5.4% P16P (CCC→CCG ybgE → putative inner membrane protein in cydABX‑ybgE operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,3560CG 5.4% 868.6 / 33.0 239P16P (CCC→CCGybgEputative inner membrane protein in cydABX‑ybgE operon
Reads supporting (aligned to +/- strand):  ref base C (121/105);  new base G (7/6);  total (128/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 246 total aligned reads displayed.

GCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGT  >  NC_000913/774208‑774505
                                                                                                                                                    |                                                                                                                                                     
gCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGccccc                                                                                                                                                      >  1:1469311/1‑150 (MQ=255)
      aTTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCggg                                                                                                                                                >  1:318877/1‑150 (MQ=255)
      aTTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCggg                                                                                                                                                <  2:673284/150‑1 (MQ=255)
       ttCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGgc                                                                                                                                               <  1:2721997/150‑1 (MQ=255)
            ggAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCttt                                                                                                                                          <  2:2778534/150‑1 (MQ=255)
             gAACGCTTCTTGCCTGTTCGTTTGGGTTAATCCCCGCGCTGGCGCTTTAACACGTCTAATCTGGCAATGCGGGTAAAGAATACATCTGATTTTTAATTTTATCGCGACTTTGTATGCGTTAATGTACAAGCGCCCGCAGGGGGCGCttgt                                                                                                                                         <  1:3449679/150‑3 (MQ=255)
             gAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCGCAGGGGGCGCttgt                                                                                                                                         <  2:787475/150‑3 (MQ=255)
               aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCt                                                                                                                                       >  2:456007/1‑150 (MQ=255)
               aCGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTGTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCt                                                                                                                                       >  2:2919891/1‑150 (MQ=255)
                  cttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCg                                                                                                                                    <  2:3347911/150‑1 (MQ=255)
                   ttcttGCCTGTTCGTTTGGGTTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                   <  1:3128963/150‑1 (MQ=255)
                   ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                   <  2:2696871/150‑1 (MQ=255)
                   ttcttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGt                                                                                                                                   <  2:758729/150‑1 (MQ=255)
                     cttGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGa                                                                                                                                 <  2:2644810/150‑1 (MQ=255)
                           tGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGgcgc                                                                                                                           >  2:2892373/1‑150 (MQ=255)
                           tGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACAAGTCGAATCAGGCAAACCCGGTCAAGAAGACATCTGATGTGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGtcgc                                                                                                                           >  2:2546929/1‑150 (MQ=255)
                            gTTCGTTTGGGGAAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCt                                                                                                                          >  2:3205997/1‑150 (MQ=255)
                               cGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCt                                                                                                                       <  1:2295703/150‑1 (MQ=255)
                               cGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCt                                                                                                                       <  1:1218715/150‑1 (MQ=255)
                                gTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTg                                                                                                                      >  1:3422625/1‑150 (MQ=255)
                                 tttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGTGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGt                                                                                                                     >  2:402754/1‑150 (MQ=255)
                                 tttGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGt                                                                                                                     <  2:3020437/150‑1 (MQ=255)
                                   tGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTa                                                                                                                   <  1:1255376/150‑1 (MQ=255)
                                     gggTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                 <  1:3663173/150‑1 (MQ=255)
                                     gggTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGc                                                                                                                 <  1:891136/150‑1 (MQ=255)
                                       gTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAg                                                                                                               <  2:2165188/150‑1 (MQ=255)
                                         aaTCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGa                                                                                                             <  1:2935900/150‑1 (MQ=255)
                                          aTCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                            <  2:463558/150‑1 (MQ=255)
                                          aTCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                            >  2:1333616/1‑150 (MQ=255)
                                          aTCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTACTTCGTGATGGCGCTTCTGTTAGCAGGat                                                                                                            >  2:3461097/1‑150 (MQ=255)
                                              ccGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtat                                                                                                        >  1:2811729/1‑150 (MQ=255)
                                              ccGCGCTGGCGCGTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACGGTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCGTTCCTGCGTGCTGGCGCTTCGGTTAGCAGGatgtat                                                                                                        >  1:2010807/1‑150 (MQ=255)
                                                gcgcTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGatgtatgt                                                                                                      >  1:803449/1‑150 (MQ=255)
                                                   cTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGtttt                                                                                                   <  1:4070906/150‑1 (MQ=255)
                                                    tGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTg                                                                                                  <  1:1818142/150‑1 (MQ=255)
                                                      gcgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTggg                                                                                                >  2:2091901/1‑150 (MQ=255)
                                                      gcgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTggg                                                                                                <  1:3027669/150‑1 (MQ=255)
                                                      gcgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTggg                                                                                                <  1:1233830/150‑1 (MQ=255)
                                                       cgcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGa                                                                                               <  1:2563686/150‑1 (MQ=255)
                                                        gcTTGAACACGTCGAATCAGGCAAAGCGGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTGTTGGGAc                                                                                              >  1:3217276/1‑150 (MQ=255)
                                                        gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                              >  1:3608539/1‑150 (MQ=255)
                                                        gcTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAc                                                                                              <  2:125681/150‑1 (MQ=255)
                                                         cTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAcc                                                                                             >  1:2254/1‑150 (MQ=255)
                                                            gAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGAcc                                                                                             <  1:2891637/147‑1 (MQ=255)
                                                            gAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCAt                                                                                          <  1:665047/150‑1 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCATGATGTATGTTTTGGGACCCATCTCGtt                                                                                    <  2:3123465/150‑1 (MQ=255)
                                                                  gTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtt                                                                                    <  1:3684895/150‑1 (MQ=255)
                                                                   tCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCCTCTCGttt                                                                                   >  1:538056/1‑150 (MQ=255)
                                                                   tCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGttt                                                                                   >  1:92689/1‑150 (MQ=255)
                                                                   tCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGttt                                                                                   >  2:698658/1‑150 (MQ=255)
                                                                    cGAATCAGGCAAAGCCGGTCAAGAAGACATCTGAAGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCTCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGtttt                                                                                  <  2:1968646/150‑1 (MQ=255)
                                                                      aaTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGc                                                                                >  2:1229745/1‑150 (MQ=255)
                                                                      aaTCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGc                                                                                >  2:968335/1‑150 (MQ=255)
                                                                        tCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcg                                                                              <  1:376579/150‑1 (MQ=255)
                                                                           ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGa                                                                           >  1:1011162/1‑150 (MQ=255)
                                                                           ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGa                                                                           >  2:3993198/1‑150 (MQ=255)
                                                                           ggCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGa                                                                           <  2:2811729/150‑1 (MQ=255)
                                                                            gCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                          >  1:2116436/1‑150 (MQ=255)
                                                                            gCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGaa                                                                          >  1:3994239/1‑150 (MQ=255)
                                                                            gCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCATTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATATCGTTTTGCAGCGaa                                                                          >  1:3817176/1‑150 (MQ=255)
                                                                                    ggTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcg                                                                              >  2:424789/1‑138 (MQ=255)
                                                                                    ggTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCcg                                                                              <  1:424789/138‑1 (MQ=255)
                                                                                        aagaagACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACt                                                              >  2:3552820/1‑150 (MQ=255)
                                                                                          gaagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTgg                                                            <  2:1463867/150‑1 (MQ=255)
                                                                                           aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGa                                                           <  2:932691/150‑1 (MQ=255)
                                                                                           aagaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGa                                                           <  1:3244608/150‑1 (MQ=255)
                                                                                             gaCATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGaaa                                                         <  2:2391805/150‑1 (MQ=255)
                                                                                                aTCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCt                                                      <  2:1133932/150‑1 (MQ=255)
                                                                                                 tCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTg                                                     <  1:3830026/150‑1 (MQ=255)
                                                                                                  cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                    >  2:2122144/1‑150 (MQ=255)
                                                                                                  cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                    >  2:1647170/1‑150 (MQ=255)
                                                                                                  cTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTgg                                                    <  1:3554284/150‑1 (MQ=255)
                                                                                                      tgaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATg                                                <  1:2126829/150‑1 (MQ=255)
                                                                                                       gaGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATgg                                               <  1:1627164/150‑1 (MQ=255)
                                                                                                          tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGttt                                            >  2:3774612/1‑150 (MQ=255)
                                                                                                          tAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTACATCGTGATGGCGCAACTGTGAGCAGGATGTATGTTTTTGGACCCAACGCGTTTTTCCGCGAAGACCACTGAACTGGAAATCTGGCATGGgtt                                            >  1:1311791/1‑150 (MQ=255)
                                                                                                             gATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTTGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCAGGCATGGTTTAtt                                         >  1:240115/1‑150 (MQ=255)
                                                                                                             gATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTAtt                                         <  1:1119220/150‑1 (MQ=255)
                                                                                                                 aTCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTg                                     <  2:1716534/150‑1 (MQ=255)
                                                                                                                   cgcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAt                                   <  1:2197612/150‑1 (MQ=255)
                                                                                                                    gcgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCAAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGAtg                                  >  2:1943978/1‑150 (MQ=255)
                                                                                                                     cgACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCGGAtgt                                 >  2:3677267/1‑150 (MQ=255)
                                                                                                                      gACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTTGGTCCCAGCTCGTTTTTCCGCGAAGACCAGTGTACTGGAAATCTGGGATGGTTTATTTGTGAtgtg                                >  1:1750519/1‑150 (MQ=255)
                                                                                                                       aCTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTTGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGCAATCTGGCATGGTTTATTGCTGATGTgg                               >  1:2515098/1‑150 (MQ=255)
                                                                                                                       aCTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                               <  1:2714333/150‑1 (MQ=255)
                                                                                                                       aCTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTgg                               <  2:2553178/150‑1 (MQ=255)
                                                                                                                                gCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgt                      <  2:532902/150‑1 (MQ=255)
                                                                                                                                 cGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTTGGACCCAACTCGTTTTTCCGCGAAGACCAGTGAACTGGAAAAATGGCATGGTTTATTGCTGATGTGGGCCGTCtgtg                     >  1:2994589/1‑150 (MQ=255)
                                                                                                                                 cGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCtgtg                     >  2:3870231/1‑150 (MQ=255)
                                                                                                                                         gACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGtgtg             >  2:3730976/1‑150 (MQ=255)
                                                                                                                                          gaaaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGa            >  2:1840182/3‑150 (MQ=255)
                                                                                                                                           aaaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAt           >  2:228225/2‑150 (MQ=255)
                                                                                                                                           aaaGCGCCCGCAGGGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAt           >  2:1185834/2‑150 (MQ=255)
                                                                                                                                            aaGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGAtt          <  2:4052542/150‑1 (MQ=255)
                                                                                                                                             aGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTc         >  2:1468687/1‑150 (MQ=255)
                                                                                                                                               cgcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCAc       >  2:2887/1‑150 (MQ=255)
                                                                                                                                                gcCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACg      >  2:1495421/1‑150 (MQ=255)
                                                                                                                                                 cccccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACgg     <  1:1009645/150‑1 (MQ=255)
                                                                                                                                                  ccccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGc    >  2:3903521/1‑150 (MQ=255)
                                                                                                                                                   cccTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCg   >  2:639454/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCATGGATTCTGGGAACGCTTCTTGCCTGTTCGTTTGGGGTAATCACCGCGCTGGCGCTTGAACACGTCGAATCAGGCAAAGCCGGTCAAGAAGACATCTGATGAGTAAGATTATCGCGACTTTGTATGCGGTAATGGACAAGCGCCCCCTGCGGGCGCTTTCCTTCGTGATGGCGCTTCTGTTAGCAGGATGTATGTTTTGGGACCCATCTCGTTTTGCCGCGAAGACCAGTGAACTGGAAATCTGGCATGGTTTATTGCTGATGTGGGCCGTCTGTGCTGGTGTGATTCACGGCGT  >  NC_000913/774208‑774505

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: