Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 850,995 T→G 6.4% intergenic (+30/+19) ompX → / ← opgE outer membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913850,9950TG 6.4% 303.8 / 11.2 111intergenic (+30/+19)ompX/opgEouter membrane protein X/OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base T (39/65);  new base G (4/3);  total (43/68)
Fisher's exact test for biased strand distribution p-value = 4.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.08e-01

Only 100 of 112 total aligned reads displayed.

TGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTC  >  NC_000913/850848‑851143
                                                                                                                                                   |                                                                                                                                                    
tGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCg                                                                                                                                                    >  2:20546/1‑150 (MQ=255)
 gCGTGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgg                                                                                                                                                   <  1:512832/150‑1 (MQ=255)
   ggggCTGCAGTTCAACCCGATGTAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgggg                                                                                                                                                 <  1:1089090/150‑1 (MQ=255)
   gggTCTGCAGATCAACCCGATGGAAAACGTTGCGCTGGACTTCTCTTACGAGCAGCGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCgggg                                                                                                                                                 <  2:1181576/150‑1 (MQ=255)
     gTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCg                                                                                                                                               <  2:841307/150‑1 (MQ=255)
        tGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAt                                                                                                                                            >  2:801247/1‑150 (MQ=255)
        tGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAt                                                                                                                                            >  2:553648/1‑150 (MQ=255)
         gCATTTCAACCCGATGGAAAACGTTGCTCTTGACTTCTCTTACTAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGTTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtt                                                                                                                                           <  2:796870/150‑1 (MQ=255)
         gCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACATGGATTGACGGTGTTGGTTACCGCTTATAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtt                                                                                                                                           >  2:1091116/1‑150 (MQ=255)
          cAGTTCAACCCGATGGAAAACGTGGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAttt                                                                                                                                          <  1:1528659/150‑1 (MQ=255)
           aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtttt                                                                                                                                         <  1:820881/150‑1 (MQ=255)
           aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGAGAGGCGGAtttt                                                                                                                                         <  1:1057846/150‑1 (MQ=255)
           aGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGGAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAtttt                                                                                                                                         <  2:1434148/150‑1 (MQ=255)
            gtgcAACCCGATGGACACCGTTGCTCGGGACTTCTCTTACGAGCCGAGCCGTATTCGTAGCGTTGACGTAGGCACCGGGATTGCCGGCGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCCCGCGAGGCGGAttttt                                                                                                                                        <  1:1595350/147‑1 (MQ=255)
            gTTCAACCCGATGGAACACGTTGCTCTGGACTGCTCTTACGAGGAGAGCCGTATTCGTCGCGCTGACGTAGGCTCCTGGATGGCCGGTGTTGGGTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                        <  2:1530356/150‑1 (MQ=255)
            gTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGAttttt                                                                                                                                        >  1:196200/1‑150 (MQ=255)
             ttCCACCCGCTGGAAAACGTTGCGCTGGTCTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                                                                       <  2:134773/150‑1 (MQ=255)
             gtCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTg                                                                                                                                       <  1:349271/149‑1 (MQ=255)
                   ccGATGGAAAACTTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTa                                                                                                                                 <  2:164045/150‑1 (MQ=255)
                   ccGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTa                                                                                                                                 >  1:1323814/1‑150 (MQ=255)
                    ggATGGAATACGTTGCTCTGGACTTCTCTTAATATCATAGTCGTATTCGTAGCGTTGACGTAGTCACCTGGATTGACGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                                <  2:255553/149‑1 (MQ=255)
                    cGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTTTTTTTaa                                                                                                                                >  1:1822718/1‑150 (MQ=255)
                    cGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTaa                                                                                                                                <  1:1407255/150‑1 (MQ=255)
                     gATGGAAAACGTTGCTCTGGCCTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                                                                               <  1:502688/150‑1 (MQ=255)
                     gATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAg                                                                                                                               <  1:434781/150‑1 (MQ=255)
                         gAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGttt                                                                                                                           >  1:1000580/1‑150 (MQ=255)
                             ccGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCggg                                                                                                                       <  1:973410/149‑1 (MQ=255)
                             aCGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCggg                                                                                                                       >  2:1229562/1‑150 (MQ=255)
                              cGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGt                                                                                                                      <  2:1757049/150‑1 (MQ=255)
                               gTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTc                                                                                                                     <  1:1025920/150‑1 (MQ=255)
                                ttGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCg                                                                                                                    <  1:472168/150‑1 (MQ=255)
                                   ctctGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGaaa                                                                                                                 >  2:792517/1‑150 (MQ=255)
                                        gACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATc                                                                                                            >  2:1797395/1‑150 (MQ=255)
                                          cTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCgg                                                                                                          <  1:1284251/150‑1 (MQ=255)
                                           ttCTCTTACGAGCAGCGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGt                                                                                                         >  1:251543/1‑150 (MQ=255)
                                             ctctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTc                                                                                                       <  1:1748196/150‑1 (MQ=255)
                                              tctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTcc                                                                                                      <  1:525883/150‑1 (MQ=255)
                                               ctTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCg                                                                                                     >  1:1290482/1‑150 (MQ=255)
                                                      gCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGt                                                                                              <  2:14233/150‑1 (MQ=255)
                                                        atagCCGTATTCGTAGCGTTGACGTAGGCACCTGGAATGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCGGAGATGGTTg                                                                                            <  2:169681/148‑1 (MQ=255)
                                                         gagCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTATGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGt                                                                                           <  2:252112/150‑1 (MQ=255)
                                                           gCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATAGGTTCCGAGATGGTTGTaa                                                                                         >  2:277468/1‑150 (MQ=255)
                                                             cGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATc                                                                                       <  1:1224644/150‑1 (MQ=255)
                                                             cGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGATTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGGCGAAAATATCGGTTCCGAGATGGTTGTAATc                                                                                       <  2:344059/150‑1 (MQ=255)
                                                              gTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATcc                                                                                      >  1:789619/1‑150 (MQ=255)
                                                                aTTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTTATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCAc                                                                                    <  1:20796/150‑1 (MQ=255)
                                                                    gTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTc                                                                                <  1:1512505/150‑1 (MQ=255)
                                                                      aGCGTTGACGTAGGCACCTGGATCGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTg                                                                              >  2:878841/1‑150 (MQ=255)
                                                                       gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                             <  1:1032255/150‑1 (MQ=255)
                                                                       gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                             <  1:407950/150‑1 (MQ=255)
                                                                       gCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCGCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGt                                                                             <  1:932472/150‑1 (MQ=255)
                                                                        cggtGACGTAGCCACCAGGAGTGCCGGTTGTTGTTACCGCTTCTAATCACTTTGGGGATATAAAAAATCCGCATCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGGCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTa                                                                            <  2:158912/147‑1 (MQ=255)
                                                                        ccttGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTa                                                                            <  1:1522611/148‑1 (MQ=255)
                                                                        cGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTa                                                                            <  2:424687/150‑1 (MQ=255)
                                                                          ttGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAAc                                                                          <  1:1242155/150‑1 (MQ=255)
                                                                                  ggCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGtt                                                                  >  1:146152/1‑150 (MQ=255)
                                                                                    cACCTTGATTGCCGGTGTTGGTTACCGCTTCTAATCCCTTTCGCGAAATAAAAACTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGGTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                                <  2:297165/150‑1 (MQ=255)
                                                                                    cACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATTGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                                >  2:501762/1‑150 (MQ=255)
                                                                                    cACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                                <  1:1022664/150‑1 (MQ=255)
                                                                                    cACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATACTCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTGCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTgg                                                                <  1:687785/150‑1 (MQ=255)
                                                                                        tGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGAt                                                            <  2:333371/150‑1 (MQ=255)
                                                                                          gATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGt                                                          <  2:1399016/150‑1 (MQ=255)
                                                                                           aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta                                                         >  1:1074101/1‑150 (MQ=255)
                                                                                           aTTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGta                                                         <  1:582790/150‑1 (MQ=255)
                                                                                                   tgtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGcc                                                 >  2:4300/1‑150 (MQ=255)
                                                                                                    gtTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGccc                                                <  1:1811596/150‑1 (MQ=255)
                                                                                                      tGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCgg                                              >  2:1020086/1‑150 (MQ=255)
                                                                                                         gttcggCTTTTAATCTCTTTGGTGATATAAATAATACGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATGGGCCCGGCtt                                           <  2:339004/145‑1 (MQ=255)
                                                                                                           aCCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttc                                         <  1:324486/150‑1 (MQ=255)
                                                                                                            ccGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGCAAGTTGGTGATGTATATCGGCCCGGCTttct                                        >  2:1357606/1‑150 (MQ=255)
                                                                                                              gCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTATGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGAAATCCACATTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTttcttc                                      >  1:1449413/1‑150 (MQ=255)
                                                                                                              gCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCAACCTTCTGTAAATGTAAGTTGGTGATGTATATCGGCCCGGCTttctta                                      >  1:1442203/1‑149 (MQ=255)
                                                                                                               cTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTtcttct                                     >  1:504955/1‑150 (MQ=255)
                                                                                                                ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTTTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTtctc                                    <  2:1605301/150‑1 (MQ=255)
                                                                                                                ttCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTtctc                                    <  1:1797395/150‑1 (MQ=255)
                                                                                                                  cTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAg                                  <  2:732481/150‑1 (MQ=255)
                                                                                                                  cTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAg                                  <  2:546951/150‑1 (MQ=255)
                                                                                                                  cTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATAGGTTCCTAGATGGTTGTAATCAACCTTCTGTAACTGGAAGTTGGTGATGTATATCTCCACGGCTTTCTTCTCaa                                  >  2:113710/1‑149 (MQ=255)
                                                                                                                   tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                                 >  2:1716533/1‑150 (MQ=255)
                                                                                                                   tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                                 <  1:1371466/150‑1 (MQ=255)
                                                                                                                   tAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCAAACTTCTGTAACTGGAAGTAGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                                 >  1:1219895/1‑150 (MQ=255)
                                                                                                                     aTCACTTTGGTTATATAAAAAATCCGCCTCTCGGGGAGTATTTTTGTTTTTAAGGTTTCGGGGCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGata                               <  2:268436/150‑1 (MQ=255)
                                                                                                                           ttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGATATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACg                         <  1:1091116/150‑1 (MQ=255)
                                                                                                                           ttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACg                         <  1:876095/150‑1 (MQ=255)
                                                                                                                           ttGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAAACCACCTTCTGTAAATGGAAGTTGTTGATTTATATCTGCCCGGATTTCTTCTCAGATATAAACg                         >  2:1094592/1‑150 (MQ=255)
                                                                                                                             ggTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACggg                       <  1:20546/150‑1 (MQ=255)
                                                                                                                             ggTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACggg                       <  2:1712482/150‑1 (MQ=255)
                                                                                                                               tGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTa                     >  1:1077213/1‑150 (MQ=255)
                                                                                                                                gATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTAt                    >  2:875162/1‑150 (MQ=255)
                                                                                                                                   ataAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                 >  2:692222/1‑150 (MQ=255)
                                                                                                                                   aacaaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                 >  2:363342/4‑150 (MQ=255)
                                                                                                                                    ttataaaTGCGCCTCTAGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATTTCGGTGCCGAGATGGTTGTAATCCACCTTCTGTATCTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGa                <  1:305064/146‑1 (MQ=255)
                                                                                                                                     caaaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGTTTTCGGTGCGAAAATATCGGTGCCGAGATGGTTGTAATCCACCTGGGGTAACTGGAAGTTGTTTAGTTATATCGGCCCGGATTTATTGTCAGAATTAAACGGGTTTTTGAt               >  2:658516/2‑150 (MQ=255)
                                                                                                                                     aaaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATATACGGGTATTTGAt               >  2:1192561/1‑150 (MQ=255)
                                                                                                                                      aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg              <  2:1000580/150‑1 (MQ=255)
                                                                                                                                      aaaaaTCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg              >  1:1784520/1‑150 (MQ=255)
                                                                                                                                       aaaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGt             >  1:1330093/1‑150 (MQ=255)
                                                                                                                                        aaaTCCGCCCCGAGAGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGtt            >  2:925108/1‑150 (MQ=255)
                                                                                                                                           tCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAt         <  2:533152/150‑1 (MQ=255)
                                                                                                                                                  ctcGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTc  >  2:1810721/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGCGGGTCTGCAGTTCAACCCGATGGAAAACGTTGCTCTGGACTTCTCTTACGAGCAGAGCCGTATTCGTAGCGTTGACGTAGGCACCTGGATTGCCGGTGTTGGTTACCGCTTCTAATCACTTTGGTGATATAAAAAATCCGCCTCTCGGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTC  >  NC_000913/850848‑851143

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