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

GTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGCC  >  NC_000913/1045585‑1045879
                                                                                                                                                   |                                                                                                                                                    
gTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAAc                                                                                                                                                   <  1:503911/150‑1 (MQ=255)
gTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGATCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGaaa                                                                                                                                                   >  2:577596/1‑149 (MQ=255)
    gCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACGATTTTGTCCATGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTCCTGGTCAGAATTTATCGACCAGTAACAAgg                                                                                                                                               <  2:305309/150‑1 (MQ=38)
     ccTGCGACATAATCGGTTTGGCGCGCTACTGATCAATCAATACTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCTATGACATCCTTGTTCCTGGTCAGAATTTATCGACCAGTAACAAgga                                                                                                                                              <  2:373322/150‑2 (MQ=37)
     ccTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAggg                                                                                                                                              <  1:773235/150‑1 (MQ=255)
     ccTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGGCAGGTAGTTCAATGCCCTCCTTGTTACTGGTCGATAAATTCTGCACAGGAACAAgtg                                                                                                                                              >  2:649545/1‑148 (MQ=255)
      cTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGt                                                                                                                                             >  2:580738/1‑150 (MQ=255)
      cTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGt                                                                                                                                             <  2:1271950/150‑1 (MQ=255)
         cGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACa                                                                                                                                          <  2:415000/150‑1 (MQ=255)
         cGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACa                                                                                                                                          <  2:1067216/150‑1 (MQ=255)
           aaaaaaTCTGTTTGGCGCGCAAAGGATCAAACCTTCCTGTCGTGACCGGATAGATATTGACCATTTTGTCCAGGTCATACTTATTGTCAGGTAGATCAATGACATCCTTTTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATg                                                                                                                                        <  1:463710/146‑1 (MQ=255)
            cATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGc                                                                                                                                       <  1:1314890/150‑1 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATTAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGAAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTg                                                                                                                                     >  2:835974/1‑150 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATTAATCAATCCTGGCGTGAAGGGATAGATATTGACCAGTTTGTCCAGGTCATACTGATTGTCAGGGAGTTCAATGCCATCCTTGTTACTGGTCGATAAATTCTGCCCAGGAACAAGGGTACATgagg                                                                                                                                     >  2:831805/1‑147 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTg                                                                                                                                     >  1:298321/1‑150 (MQ=255)
              tAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTg                                                                                                                                     <  1:533142/150‑1 (MQ=255)
                 tCGCTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCa                                                                                                                                  <  1:714358/150‑1 (MQ=255)
                   ggTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCaaa                                                                                                                                <  2:470509/150‑1 (MQ=255)
                           gcgcAAGTGATCAATCTATGCTGGCGGGACGGGATAGTAAGTGACCATTTTGTCAAGGTCATAGTGATTGTCAGGTAGTGGAATCACATCCATGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAACCACAGAAc                                                                                                                        <  2:1173647/150‑1 (MQ=255)
                            cgcAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGAACATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACa                                                                                                                       >  2:165247/1‑150 (MQ=255)
                             gcAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAg                                                                                                                      >  2:749055/1‑150 (MQ=255)
                              cAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGc                                                                                                                     <  1:412028/150‑1 (MQ=255)
                                aCTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAATGGTACATGCTGTCAAAGACAGAACAGCCa                                                                                                                   >  2:830558/1‑150 (MQ=255)
                                 cTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGTGTACATGCTGTCAAAGACAGAACAGCCaa                                                                                                                  >  1:6821/1‑150 (MQ=255)
                                  tGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                                                                                                                 <  1:215121/150‑1 (MQ=255)
                                  tGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAt                                                                                                                 >  2:586983/1‑150 (MQ=255)
                                   gATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTAAGGTAGTTCAATGACATCCTTGATAATGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAAtt                                                                                                                >  1:972026/1‑149 (MQ=255)
                                          caatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAc                                                                                                         <  2:777724/150‑1 (MQ=255)
                                          caatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAc                                                                                                         <  2:639568/150‑1 (MQ=255)
                                           aatcCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTAcc                                                                                                        <  1:1282399/150‑1 (MQ=255)
                                                      gACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaa                                                                                             >  1:391900/1‑150 (MQ=255)
                                                       aCGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaa                                                                                            >  2:371666/1‑150 (MQ=255)
                                                       aCGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaa                                                                                            <  1:758678/150‑1 (MQ=255)
                                                          ggATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATg                                                                                         <  1:377291/150‑1 (MQ=255)
                                                          ggATAGATATTGACCATTTTGTCCAGGTCATAATGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATg                                                                                         <  2:84845/150‑1 (MQ=255)
                                                           tatagataTTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTaaaatgtt                                                                                        <  2:250918/149‑5 (MQ=255)
                                                             tagataTTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTc                                                                                      >  2:189253/1‑150 (MQ=255)
                                                                 tatTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGGCGATAAATTCTGACCAGGAACAAGGGTACATGCTGGCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCtttt                                                                                  <  1:649545/150‑1 (MQ=255)
                                                                  atTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCttttt                                                                                 <  1:453935/150‑1 (MQ=255)
                                                                     gACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCATTGTCAGTACCGAACACATAAAATTGTTCTTTTTCAt                                                                              >  2:104636/1‑150 (MQ=255)
                                                                       ccATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTg                                                                            >  2:1081403/1‑150 (MQ=255)
                                                                            ttGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTta                                                                       <  1:310904/150‑1 (MQ=255)
                                                                               tCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATGGGTTTatat                                                                    >  2:661012/1‑150 (MQ=255)
                                                                                ccAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGCACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTAtttt                                                                   >  2:622094/1‑150 (MQ=255)
                                                                                 cAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATAttt                                                                  >  2:1267240/1‑150 (MQ=255)
                                                                                    gTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATc                                                               <  2:26961/150‑1 (MQ=255)
                                                                                      cATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAg                                                             <  1:606578/150‑1 (MQ=255)
                                                                                       aTACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGt                                                            >  2:46802/1‑150 (MQ=255)
                                                                                          cTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaa                                                         >  2:227855/1‑150 (MQ=255)
                                                                                           tGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaat                                                        <  1:189253/150‑1 (MQ=255)
                                                                                            gATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATTCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaatg                                                       <  1:369799/150‑1 (MQ=255)
                                                                                            gATTGTCAGGTAGTTCAATGAAATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCaatg                                                       <  1:979838/150‑1 (MQ=255)
                                                                                                     gTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTAt                                              >  2:1205895/1‑150 (MQ=255)
                                                                                                          tCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAAtt                                         >  1:1017907/1‑150 (MQ=255)
                                                                                                          tCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAAtt                                         >  1:1220384/1‑150 (MQ=255)
                                                                                                           cAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTa                                        >  2:981885/1‑150 (MQ=255)
                                                                                                                  aTCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAgg                                 >  1:350601/1‑150 (MQ=255)
                                                                                                                   tacTTGTTAATGGTCGATAAATTCTGACCAGGAACAAGGTTACATGCTGTCAAAGACAGAACAGACAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                <  1:69832/148‑1 (MQ=255)
                                                                                                                   tCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                <  1:605903/150‑1 (MQ=255)
                                                                                                                   cccTTGTTCCTGGTCAGAATTTATCGACCAGTAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGt                                >  2:1199721/2‑150 (MQ=38)
                                                                                                                    ccTTGTTCCTGGTCAGAATTTATCGACAAGTAACTAGGGTACATGCTGACAAAGCAAGATCAGCCAATGTCAGTACCGAAAACTT‑‑‑AATGTTCTTATACAATCGTAAATATGTATCAGTCAATCTAATGTTATGAAATACGTAAGGTAAtt                            >  2:564184/1‑150 (MQ=255)
                                                                                                                    ccTTGTTACTGGTCGATAAATTCTGACCAGGAACTAGGGTACATGCTGTCAAAGACAGAACAGCCAAAGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGtt                               <  2:549432/150‑1 (MQ=255)
                                                                                                                       tGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTAtt                            >  2:2528/1‑150 (MQ=255)
                                                                                                                        gTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTACATAGGTTTATATTTATCAGTCAATGTAATTTAATGAAATACTTCAGGTTAttt                           >  2:477216/1‑149 (MQ=255)
                                                                                                                        gTTACTGGTCGATAAATTCTGAACAGGAACAAGTGTACATTCTGTCAAAGACAGAAAAGACAATTTCAGTACCTAAAACTTAAAAATGTTCTTTTTCATTGTTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTg                           <  1:24648/150‑1 (MQ=255)
                                                                                                                        gTTACTGGGCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATATATCAGTCAATGTAATGTTATGTATTACGTCAGGTTATTg                           <  2:870338/150‑1 (MQ=255)
                                                                                                                          tACTTGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACATAACATCAAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTc                         <  1:320746/150‑1 (MQ=255)
                                                                                                                            cTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCa                       >  1:1270539/1‑150 (MQ=255)
                                                                                                                            cTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCa                       >  2:1001006/1‑150 (MQ=255)
                                                                                                                              ggTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAc                     >  2:367834/1‑150 (MQ=255)
                                                                                                                              ggTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTATATATTTATCACTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCTAc                     >  2:222548/1‑150 (MQ=255)
                                                                                                                               gTCTATAAATTCTGACCATGAACAAGGGTACATGCTGTCAAAGAAAGAACAGCCAATGTCAGTACCCGAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAcc                    <  2:112200/150‑1 (MQ=255)
                                                                                                                               gTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAAcc                    <  1:371666/150‑1 (MQ=255)
                                                                                                                                 cGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCAAACCgg                  >  1:98263/1‑150 (MQ=255)
                                                                                                                                   aTAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTAc               <  2:965098/150‑1 (MQ=255)
                                                                                                                                    tAAATTCTGACCAGTTTCAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGACAGGTTATTGTCCAACCGGTAc               <  1:475900/150‑1 (MQ=255)
                                                                                                                                    tAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTAc               <  2:742004/150‑1 (MQ=255)
                                                                                                                                      aaTTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACtt             >  1:538387/1‑150 (MQ=255)
                                                                                                                                      aaTTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACAGGTACtt             >  1:1097818/1‑150 (MQ=255)
                                                                                                                                       aTTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTtc            >  2:132118/1‑150 (MQ=255)
                                                                                                                                       aTTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTtc            >  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)
                                                                                                                                                   |                                                                                                                                                    
GTTAGCCTGCGACATAATCGGTTTGGCGCGCAACTGATCAATCAATCCTGGCGTGACGGGATAGATATTGACCATTTTGTCCAGGTCATACTGATTGTCAGGTAGTTCAATGACATCCTTGTTACTGGTCGATAAATTCTGACCAGGAACAAGGGTACATGCTGTCAAAGACAGAACAGCCAATGTCAGTACCGAAAACTTAAAAATGTTCTTTTTCATTGGTTTATATTTATCAGTCAATGTAATGTTATGAATTACGTCAGGTTATTGTCCAACCGGTACTTCTCGCTCTGCC  >  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: