Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,045,731 G→T 5.4% P25T (CCT→ACT)  gfcE ← putative O‑antigen capsule outer membrane auxillary protein export channel

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,045,7310GT 5.4% 224.6 / 2.5 83P25T (CCT→ACT) gfcEputative O‑antigen capsule outer membrane auxillary protein export channel
Reads supporting (aligned to +/- strand):  ref base G (40/38);  new base T (2/3);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

GTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGCC  >  NC_000913/1045585‑1045879
                                                                                                                                                   |                                                                                                                                                    
gTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAAc                                                                                                                                                   <  1:503911/150‑1 (MQ=255)
gTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGATCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGaaa                                                                                                                                                   >  2:577596/1‑149 (MQ=255)
    gCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACGATTTTGTCCATGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTCCTGGTC‑AGAATTTATCGACCAGTAACAAgg                                                                                                                                               <  2:305309/150‑1 (MQ=38)
     ccTGCGACATAATCGGTTTGGCGCGCTACTGATCAATCAATACTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCTATGACATCCTTGTTCCTGGTC‑AGAATTTATCGACCAGTAACAAgga                                                                                                                                              <  2:373322/150‑2 (MQ=37)
     ccTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAggg                                                                                                                                              <  1:773235/150‑1 (MQ=255)
     ccTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGGCAGGTAGTTCAATGCCCTCCTTGTTACTGGTCGATAAATTCT‑GCACAGGAACAAgtg                                                                                                                                              >  2:649545/1‑148 (MQ=255)
      cTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGt                                                                                                                                             >  2:580738/1‑150 (MQ=255)
      cTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGt                                                                                                                                             <  2:1271950/150‑1 (MQ=255)
         cGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACa                                                                                                                                          <  2:415000/150‑1 (MQ=255)
         cGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACa                                                                                                                                          <  2:1067216/150‑1 (MQ=255)
           aaaaaaTCTGTTTGGCGCGCAAAGGATCAAACCTTCCTGTCGTGACCGGATAGATATTGACCATTTTGTCCAGGTCATACTTATTGTCAGGTAGATCAATGACATCCTTTTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATg                                                                                                                                        <  1:463710/146‑1 (MQ=255)
            cATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGc                                                                                                                                       <  1:1314890/150‑1 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATTAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGAAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTg                                                                                                                                     >  2:835974/1‑150 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATTAATCAATCCTGGCGTGAAGGGATAGATATTGACCAGTTTGTCCAGGTCATACTGATTGTCAGGGAGTTCAATGCCATCCTTGTTACTGGTCGATAAATTCT‑GCCCAGGAACAAGGGTACATgagg                                                                                                                                     >  2:831805/1‑147 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTg                                                                                                                                     >  1:298321/1‑150 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTg                                                                                                                                     <  1:533142/150‑1 (MQ=255)
                 tCGCTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCa                                                                                                                                  <  1:714358/150‑1 (MQ=255)
                   ggTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCaaa                                                                                                                                <  2:470509/150‑1 (MQ=255)
                           gcgcAAGTGATCAATCTATGCTGGCGGGACGGGATAGTAAGTGACCATTTTGTCAAGGTCATAGTGATTGTCAGGTAGTGGAATCACATCCATGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAACCACAGAAc                                                                                                                        <  2:1173647/150‑1 (MQ=255)
                            cgcAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGAACATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACa                                                                                                                       >  2:165247/1‑150 (MQ=255)
                             gcAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAg                                                                                                                      >  2:749055/1‑150 (MQ=255)
                              cAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGc                                                                                                                     <  1:412028/150‑1 (MQ=255)
                                aCTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAATGGTACATGCTGTCAAAGACAGAACAGCCa                                                                                                                   >  2:830558/1‑150 (MQ=255)
                                 cTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGTGTACATGCTGTCAAAGACAGAACAGCCaa                                                                                                                  >  1:6821/1‑150 (MQ=255)
                                  tGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                                                                                                                 <  1:215121/150‑1 (MQ=255)
                                  tGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                                                                                                                 >  2:586983/1‑150 (MQ=255)
                                   gATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTAAGGTAGTTCAATGACATCCTTGATAATGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAtt                                                                                                                >  1:972026/1‑149 (MQ=255)
                                          caatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAc                                                                                                         <  2:777724/150‑1 (MQ=255)
                                          caatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAc                                                                                                         <  2:639568/150‑1 (MQ=255)
                                           aatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAcc                                                                                                        <  1:1282399/150‑1 (MQ=255)
                                                      gACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaa                                                                                             >  1:391900/1‑150 (MQ=255)
                                                       aCGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaa                                                                                            >  2:371666/1‑150 (MQ=255)
                                                       aCGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaa                                                                                            <  1:758678/150‑1 (MQ=255)
                                                          ggATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATg                                                                                         <  1:377291/150‑1 (MQ=255)
                                                          ggATAGATATTGACCATTTTGTCCAGGTCATAATGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATg                                                                                         <  2:84845/150‑1 (MQ=255)
                                                           tatagataTTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaatgtt                                                                                        <  2:250918/149‑5 (MQ=255)
                                                             tagataTTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTc                                                                                      >  2:189253/1‑150 (MQ=255)
                                                                 tatTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGGCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGGCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCtttt                                                                                  <  1:649545/150‑1 (MQ=255)
                                                                  atTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCttttt                                                                                 <  1:453935/150‑1 (MQ=255)
                                                                     gACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCATTGTCAGTACCGAACACATAAAATTGTTCTTTTTCAt                                                                              >  2:104636/1‑150 (MQ=255)
                                                                       ccATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTg                                                                            >  2:1081403/1‑150 (MQ=255)
                                                                            ttGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTta                                                                       <  1:310904/150‑1 (MQ=255)
                                                                               tCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATGGGTTTatat                                                                    >  2:661012/1‑150 (MQ=255)
                                                                                ccAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGCACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTAtttt                                                                   >  2:622094/1‑150 (MQ=255)
                                                                                 cAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATAttt                                                                  >  2:1267240/1‑150 (MQ=255)
                                                                                    gTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATc                                                               <  2:26961/150‑1 (MQ=255)
                                                                                      cATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAg                                                             <  1:606578/150‑1 (MQ=255)
                                                                                       aTACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGt                                                            >  2:46802/1‑150 (MQ=255)
                                                                                          cTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaa                                                         >  2:227855/1‑150 (MQ=255)
                                                                                           tGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaat                                                        <  1:189253/150‑1 (MQ=255)
                                                                                            gATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATTCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaatg                                                       <  1:369799/150‑1 (MQ=255)
                                                                                            gATTGTCAGGTAGTTCAATGAAATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaatg                                                       <  1:979838/150‑1 (MQ=255)
                                                                                                     gTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTAt                                              >  2:1205895/1‑150 (MQ=255)
                                                                                                          tCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAAtt                                         >  1:1017907/1‑150 (MQ=255)
                                                                                                          tCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAAtt                                         >  1:1220384/1‑150 (MQ=255)
                                                                                                           cAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTa                                        >  2:981885/1‑150 (MQ=255)
                                                                                                                  aTCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAgg                                 >  1:350601/1‑150 (MQ=255)
                                                                                                                   tacTTGTTAATGGTCGATAAATTCT‑GACCAGGAACAAGGTTACATGCTGTCAAAGACAGAACAGACAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                <  1:69832/148‑1 (MQ=255)
                                                                                                                   tCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                <  1:605903/150‑1 (MQ=255)
                                                                                                                   cccTTGTTCCTGGTC‑AGAATTTATCGACCAGTAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                >  2:1199721/2‑150 (MQ=38)
                                                                                                                    ccTTGTTCCTGGTC‑AGAATTTATCGACAAGTAACTAGGGTACATGCTGACAAAGCAAGATCAGCCAATGTCAGTACCGAAAACTT‑‑‑AATGTTCTTATACAATCGTAAATATGTATCAGTCAATCTAATGTTATGAAATACGTAAGGTAAtt                            >  2:564184/1‑150 (MQ=255)
                                                                                                                    ccTTGTTACTGGTCGATAAATTCT‑GACCAGGAACTAGGGTACATGCTGTCAAAGACAGAACAGCCAAAGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGtt                               <  2:549432/150‑1 (MQ=255)
                                                                                                                       tGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTAtt                            >  2:2528/1‑150 (MQ=255)
                                                                                                                        gTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTACATAGGTTTATATTTATCAGTCAATGTAATTTAATGAAATACTTCAGGTTAttt                           >  2:477216/1‑149 (MQ=255)
                                                                                                                        gTTACTGGTCGATAAATTCT‑GAACAGGAACAAGTGTACATTCTGTCAAAGACAGAAAAGACAATTTCAGTACCTAAAACTTAAAAATGTTCTTTTTCATTGTTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTg                           <  1:24648/150‑1 (MQ=255)
                                                                                                                        gTTACTGGGCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATATATCAGTCAATGTAATGTTATGTATTACGTCAGGTTATTg                           <  2:870338/150‑1 (MQ=255)
                                                                                                                          tACTTGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACATAACATCAAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTc                         <  1:320746/150‑1 (MQ=255)
                                                                                                                            cTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCa                       >  1:1270539/1‑150 (MQ=255)
                                                                                                                            cTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCa                       >  2:1001006/1‑150 (MQ=255)
                                                                                                                              ggTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAc                     >  2:367834/1‑150 (MQ=255)
                                                                                                                              ggTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTATATATTTATCACTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCTAc                     >  2:222548/1‑150 (MQ=255)
                                                                                                                               gTCTATAAATTCT‑GACCATGAACAAGGGTACATGCTGTCAAAGAAAGAACAGCCAATGTCAGTACCCGAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAcc                    <  2:112200/150‑1 (MQ=255)
                                                                                                                               gTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAcc                    <  1:371666/150‑1 (MQ=255)
                                                                                                                                 cGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCAAACCgg                  >  1:98263/1‑150 (MQ=255)
                                                                                                                                   aTAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATA‑TTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTAc               <  2:965098/150‑1 (MQ=255)
                                                                                                                                    tAAATTCT‑GACCAGTTTCAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGACAGGTTATTGTCCAACCGGTAc               <  1:475900/150‑1 (MQ=255)
                                                                                                                                    tAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTAc               <  2:742004/150‑1 (MQ=255)
                                                                                                                                      aaTTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACtt             >  1:538387/1‑150 (MQ=255)
                                                                                                                                      aaTTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACAGGTACtt             >  1:1097818/1‑150 (MQ=255)
                                                                                                                                       aTTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTtc            >  2:132118/1‑150 (MQ=255)
                                                                                                                                       aTTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTtc            >  2:132115/1‑150 (MQ=255)
                                                                                                                                             gACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGct       >  1:868727/1‑150 (MQ=255)
                                                                                                                                               ccAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGctct     >  1:884414/1‑150 (MQ=255)
                                                                                                                                                  ggAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGcc  >  1:309437/1‑150 (MQ=255)
                                                                                                                                                  ggAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGcc  >  1:1254682/1‑150 (MQ=255)
                                                                                                                                                  ggAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGcc  >  2:951683/1‑150 (MQ=255)
                                                                                                                                                  ggAACAAGGGAACATGCTCTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGcc  <  2:472170/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCT‑GACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGCC  >  NC_000913/1045585‑1045879

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: ‑