Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,005,381 A→G S174G (AGC→GGC)  pldA → outer membrane phospholipase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,005,3810AG99.3% 516.2 / ‑4.1 145S174G (AGC→GGC) pldAouter membrane phospholipase A
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (71/73);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

Only 100 of 145 total aligned reads displayed.

AGAGTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCAGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTATCAC  >  NC_000913/4005233‑4005527
                                                                                                                                                    |                                                                                                                                                  
 gagTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCt                                                                                                                                                 <  4:1262536/151‑1 (MQ=255)
 gagTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCt                                                                                                                                                 <  4:423007/151‑1 (MQ=255)
 gagTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCt                                                                                                                                                 <  1:803007/151‑1 (MQ=255)
  agTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTg                                                                                                                                                <  1:1508442/151‑1 (MQ=255)
       tCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGaa                                                                                                                                             <  1:101429/149‑1 (MQ=255)
       tCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGaa                                                                                                                                             <  4:1433897/149‑1 (MQ=255)
       tCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAAcc                                                                                                                                           >  1:1099887/1‑151 (MQ=255)
        cACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCg                                                                                                                                          <  3:416166/151‑1 (MQ=255)
           cGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCt                                                                                                                                       <  4:1405106/151‑1 (MQ=255)
           cGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCt                                                                                                                                       <  4:524396/151‑1 (MQ=255)
              ttCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCttt                                                                                                                                     >  4:1082965/1‑150 (MQ=255)
                 gTGAAACCAACTACGACCCGCAATTGTTCCTCGGTTTTGCCACCGATTACCTTTTTGCGGGTTGGACGCTGCGCGATGTGGAGAGGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTta                                                                                                                                    <  1:1047565/148‑1 (MQ=255)
                      aCCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGcc                                                                                                                            >  4:747615/1‑151 (MQ=255)
                        cAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTg                                                                                                                          >  2:346005/1‑151 (MQ=255)
                         aaCTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTg                                                                                                                          >  1:326813/1‑150 (MQ=255)
                              cGATCCGCAATTGTTCTGCGGTTTTCCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGTTATAACCACGACTCTAACGGGCGTTCCGACCCGCCCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATg                                                                                                                       <  4:326784/148‑1 (MQ=255)
                                acccGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTCCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGa                                                                                                                  <  3:247141/149‑1 (MQ=255)
                                aaCCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGGGATGTGGAGATGGGGTATAACCCCGACTCTAACGGGCGTTCCGACCCGACCCCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGa                                                                                                                  <  2:126937/151‑1 (MQ=255)
                                aaCCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGc                                                                                                                     >  1:1310062/1‑148 (MQ=255)
                                aaCCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGc                                                                                                                     >  4:383930/1‑148 (MQ=255)
                                       aTTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGt                                                                                                           >  1:920368/1‑151 (MQ=255)
                                           ttCCTCGGTTTTGCTACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAAct                                                                                                       >  4:250548/1‑151 (MQ=255)
                                             ccTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAActg                                                                                                      >  4:1407586/1‑150 (MQ=255)
                                              cTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAActggc                                                                                                    >  1:60054/1‑151 (MQ=255)
                                                cGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGCTGTGGAGATGGGGTATAACCACGACTCTAACGCGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAActggctg                                                                                                  >  4:920377/1‑151 (MQ=255)
                                                 ggTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAActggctgg                                                                                                 <  4:1290626/151‑1 (MQ=255)
                                                   ttttGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTCACGGGCGTTCCGACCAGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCGGATGGCAGAAAACGGTAACTGGCTGGt                                                                                                >  1:1477958/1‑150 (MQ=255)
                                                     ttGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATACCCCCGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAg                                                                                              <  4:756348/150‑1 (MQ=255)
                                                      tGCCACCGATTCCCGTTTGGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGa                                                                                             <  4:1724179/150‑1 (MQ=255)
                                                      tGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGa                                                                                             >  4:1780100/1‑150 (MQ=255)
                                                       gCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATACCCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGaa                                                                                            <  4:1109401/150‑1 (MQ=255)
                                                        ccACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAg                                                                                           >  4:1562658/1‑150 (MQ=255)
                                                        ccACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGt                                                                                          >  1:264934/1‑151 (MQ=255)
                                                           ccGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGa                                                                                        >  4:206220/1‑150 (MQ=255)
                                                             gATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGc                                                                                     <  2:249618/151‑1 (MQ=255)
                                                              aTTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGc                                                                                     <  4:1191014/150‑1 (MQ=255)
                                                               ttACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATACCCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGcc                                                                                    <  1:1459853/150‑1 (MQ=255)
                                                               ttACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCg                                                                                   >  4:1399754/1‑151 (MQ=255)
                                                               ttACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCg                                                                                   >  3:407971/1‑151 (MQ=255)
                                                               ttACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCAGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCg                                                                                   >  4:674813/1‑151 (MQ=255)
                                                                    gTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGt                                                                               >  1:1610287/1‑150 (MQ=255)
                                                                    gTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGt                                                                               >  3:4802/1‑150 (MQ=255)
                                                                      tttGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtg                                                                            >  3:288017/1‑151 (MQ=255)
                                                                       ttGCAGGTTGGGCGCTGCGCGATGTGGAGATGGGGTATACCCAGGACTCTAACGGGCGTTCCGCCCCGACCTCCCGGGGCGGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtg                                                                            <  4:1356360/150‑1 (MQ=255)
                                                                       ttGCAGGTTGGACGCTGCGGGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtg                                                                            <  4:243836/150‑1 (MQ=255)
                                                                       ttGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtgt                                                                           >  2:425639/1‑151 (MQ=255)
                                                                        tGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtgtg                                                                          >  3:528867/1‑151 (MQ=255)
                                                                         gcgggTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTAtgtg                                                                          <  4:990491/147‑1 (MQ=255)
                                                                            ggTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATgtggtg                                                                       <  2:49244/150‑1 (MQ=255)
                                                                              ttGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGt                                                                    >  2:77698/1‑151 (MQ=255)
                                                                               tGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGt                                                                    >  2:132116/1‑150 (MQ=255)
                                                                                ggACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTaa                                                                  <  2:218126/151‑1 (MQ=255)
                                                                                ggACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTa                                                                   <  4:1828763/150‑1 (MQ=255)
                                                                                ggACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTa                                                                   <  2:187416/150‑1 (MQ=255)
                                                                                   cGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATAc                                                               <  2:98519/151‑1 (MQ=255)
                                                                                   cGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATAc                                                               <  1:230785/151‑1 (MQ=255)
                                                                                      tGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTg                                                             >  2:382511/1‑150 (MQ=255)
                                                                                       gcgcTATGTGGAGAGGGGGTATAACCGAGACTCTAAGGGGCGTCCCGACCCGACCTCCCGCGGCTGGACCCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTg                                                             <  4:381649/149‑1 (MQ=255)
                                                                                         gcgATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACg                                                          <  4:1737130/150‑1 (MQ=255)
                                                                                          cgATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGAt                                                        >  1:1117850/1‑151 (MQ=255)
                                                                                           gATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGAt                                                        <  2:19009/150‑1 (MQ=255)
                                                                                            aTGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGAt                                                        <  1:1082499/149‑1 (MQ=255)
                                                                                            aTGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATa                                                       >  2:319616/1‑150 (MQ=255)
                                                                                             tgtgGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAAc                                                     >  2:201360/1‑151 (MQ=255)
                                                                                                ggAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCg                                                  >  4:517508/1‑151 (MQ=255)
                                                                                                 gagaTGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCgg                                                 >  1:843768/1‑151 (MQ=255)
                                                                                                      ggggTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGatat                                             <  1:1349445/150‑1 (MQ=255)
                                                                                                      ggggTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGata                                              <  1:436996/149‑1 (MQ=255)
                                                                                                       gggTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGATGATAACCCGGATATCa                                           <  4:784118/151‑1 (MQ=255)
                                                                                                        ggTATAACCACGACTCTAACGGGCGTTCCGACCCGCCCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCAc                                          <  1:20689/151‑1 (MQ=255)
                                                                                                        ggTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCAc                                          <  4:1613372/151‑1 (MQ=255)
                                                                                                         gTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCAcc                                         >  1:122333/1‑151 (MQ=255)
                                                                                                           ataACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCa                                        >  1:898017/1‑150 (MQ=255)
                                                                                                             aaCCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCaaa                                      >  1:434954/1‑150 (MQ=255)
                                                                                                             aaCCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAat                                     >  4:535266/1‑151 (MQ=255)
                                                                                                                cACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAatata                                  <  1:487490/151‑1 (MQ=255)
                                                                                                                     ctctAACGGGCGTGCCGGCCCGCCCTCCCGCGGTGGAAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATgg                               <  4:451828/149‑1 (MQ=255)
                                                                                                                      tctAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGtt                            >  1:1362136/1‑151 (MQ=255)
                                                                                                                        tAACGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGt                             <  4:1138014/148‑1 (MQ=255)
                                                                                                                         aaCGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTtact                         >  1:80809/1‑151 (MQ=255)
                                                                                                                         aaCGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTtac                          <  2:511604/150‑1 (MQ=255)
                                                                                                                           cGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTtactac                       >  1:107622/1‑151 (MQ=255)
                                                                                                                           cGGGCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTtacta                        <  4:823857/150‑1 (MQ=255)
                                                                                                                            gggcgaagagacccgaccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGCTGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTAcc                      >  1:811780/10‑151 (MQ=255)
                                                                                                                            gggCTTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTtacta                        <  4:1376148/149‑1 (MQ=255)
                                                                                                                            gggCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTAcc                      <  4:677833/151‑1 (MQ=255)
                                                                                                                              gCGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAg                    <  1:1668363/151‑1 (MQ=255)
                                                                                                                               cGTTCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTCCTACCAGc                   >  3:73327/1‑151 (MQ=255)
                                                                                                                                  tCCGACCCGACCTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCtt                 >  1:339813/1‑150 (MQ=255)
                                                                                                                                   ccgccccgccCTCCCGCGTCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTa                <  4:1036096/150‑1 (MQ=255)
                                                                                                                                     gacccgaccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTaaa              >  1:1140697/1‑150 (MQ=255)
                                                                                                                                       cccgaccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAAt            <  1:1114175/150‑1 (MQ=255)
                                                                                                                                       cccgaccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATc           <  4:1485142/151‑1 (MQ=255)
                                                                                                                                         cgaccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCgg         >  1:19928/1‑151 (MQ=255)
                                                                                                                                            ccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGc        <  3:185562/149‑1 (MQ=255)
                                                                                                                                            ccTCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTa      <  1:581931/151‑1 (MQ=255)
                                                                                                                                              tCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTATc    >  1:836616/1‑151 (MQ=255)
                                                                                                                                              tCCCGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTATc    >  1:1663712/1‑151 (MQ=255)
                                                                                                                                               cccGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTAt     <  4:1324644/149‑1 (MQ=255)
                                                                                                                                                ccGCGGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTATCAc  >  4:1910938/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGAGTCTTCACCGTTTCGTGAAACCAACTACGAACCGCAATTGTTCCTCGGTTTTGCCACCGATTACCGTTTTGCAGGTTGGACGCTGCGCGATGTGGAGATGGGGTATAACCACGACTCTAACGGGCGTTCCGACCCGACCTCCCGCAGCTGGAACCGCCTTTATACTCGCCTGATGGCAGAAAACGGTAACTGGCTGGTAGAAGTGAAGCCGTGGTATGTGGTGGGTAATACTGACGATAACCCGGATATCACCAAATATATGGGTTACTACCAGCTTAAAATCGGCTATCAC  >  NC_000913/4005233‑4005527

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